ELEMENTS OF 
HEAT TRANSFER 
AND 

INSULATION 



ELEMENTS OF 

HEAT TRANSFER 

AND 

INSULATION 

MAX JAKOB 

RGsecirch Professor of Mechanical Engineering 
Illinois Institate of Technology 


GEORGE A. HAWKINS 


Associate Professor of Mechanical Engineering 
Purdue University 


NEW YORK 

JOHN WILEY & SONS, inc. 

LONDON; CHAPMAN & HALL, limited 

1942 



Copyright, 1942 

By max JAKOB and GEORGE A. HAWKINS 


All Rights Reserved 

This book or any part thereof must not 
be reproduced in any form without 
the written permission of the publisher. 


PRINTED IN THE UNITED STATES OF AMERICA 



PREFACE 


“ We are frequently told that heat transfer is a field for gn 
Study; that it is too difficult a subject in its rational form for the i 
graduate student. If we accept this view, we may well be aecu 
accepting mechanical engineering, in the majority of its practices, 
Empirical art.” 

! Approximately three years before Dean Wilkinson publishe 
ibove statement* the authors had begun to teach regular el 
jourses in heat transfer, designed for senior mechanical and ch( 
jngineering students, following the general lines of the present 
The subjects of converting and conveying heat energy are of s: 
jnportance. The first deals with thermodynamics and is indue 
llmost every engineering curriculum; the latter is now beginning 
ooked upon favorably for inclusion in the undergraduate plan of £ 
It, therefore, seems timely to present a textbook for courses in 
gansmission, not designed for graduate level nor exhaustive in s 
jut adapted to the average junior or senior engineering student. 
{)ecial course in heat transmission is offered, this text may serve 
idditional brief book in the usual thermodynamics course. Ir 
lay the student will be given a more nearly complete picture c 
fiientific fundamentals of heat engineering. 

In order to be suited for these purposes the book is restrictp(ft 
iw basic principles of heat transfer and insulation, and to t^ir ap 
ion to simple problems. The various subjects are presented 
toad manner with examples following the derivations. 

However, scientific rigor has been retained to as great an exte 
'■ssible at the intended level of presentation. This applies espei 
t the use of consistent units, which are so important in the scienci 
jplication of heat transfer. In this matter no concessions to supei 
'actices have been made. The student is induced to insure in 
ecial case that the balance of physical dimensions is correct, 
rticularly to avoid the usual confusion of mass and force units. 
Conduction, convection, and radiation are treated separately 
d then in their combinations. A whole chapter is devoted t( 
sic property of thermal conductivity with particular attentio 

* F. L. Wilkinson, Jr., “ A Suggested Design for Mechanical Engineering C 
1 ,” Mechanical Engineering, Vol. 63, p, 681 (1941). 
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practical insulations. Unsteady state equations are made palatable i 
the student by some graphical representations. Dimensional analys 
is dealt with in the most elementary manner. Some applications of he£ 
transfer in experimental engineering and the interdependence of hej 
transfer and surface friction form the subject of the last three chapter 

In order to facilitate the solution of the problems, the Contents 
supplemented by a list of tables and a list of figures needed for numeric: 
calculations. These are followed by a list of symbols. Answers to tl 
problems will be printed separately. 

Space and purpose of the book forbid ample literature reference 
Only a very small selection of them is given at the end of the chapter 
However, every effort has been made to give credit to individuals f( 
use of their materials. Special thanks is due to Mr, H. E. Hollensb« 
Editor of Industrial Power, for the release of certain copyrighted materij 
herein contained in the form of problems. 

The authors wish to express appreciation to President H. T. Heald <; 
Illinois Institute of Technology, and Dr. A. A. Potter, Dean of th 
Schools of Engineering, Purdue University, who by initiating unde: 
graduate courses in heat transfer, encouraged our task, and to Professc 
H. L. Solberg, Head of the School of Mechanical Engineering, Purch 
University, for his keen interest during the preparation of this textbool 


Chicago, Illinois 
Lafayette, Indiana 
March, 194^ 


Max Jakob 
Geokge a. Hawkins 
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Symbol Quantity Units* 


A 

Area 

ft^ 

B 

Buoyancy 

lb ft-® 

C 

Circumference; mean coil diameter 

ft 

c 

Specific heat 

B slug-1 F-i 

Cp 

Specific heat at constant pressure related to unit volume 

B fr® p-i 

Cp 

Specific heat at constant pressure • 

B slug-1 P-i 

D 

Diameter 

ft 

E 

Density of emission 

B hr-i fr® 

E 

Electrical potential 

volt 

f 

Friction factor 


Fa 

Configuration factor or area factor 


Fe 

Emissivity factor 


Fs 

Shearing force (friction) 

lb 


Gravitational acceleration 

fthr-2 


Grashof number 


H 

Fundamental unit of heat energy 

B 

h ^ 

Heat transfer coefficient 

B hr-i fr® F-i 

I 

Intensity of radiation 

B br-i fr® 

I 

Electric current 

ampere 

k 

Thermal conductivity 

Bhr-ift-ip-i 

L 

Length, Fundamental unit of length 

ft 

1 

Distance 

ft 

1 

Latent heat 

B slug-1 

M 

Fundamental unit of mass 

slug 

M 

Mass velocity 

slug bf-i ft-® 

m 

Rate of mass flow 

slug hr-i 

n 

Frequency 

hr-i 

m) 

Nusselt number 


V 

Pressure 

lb ft-® 

{Ft) 

Prandtl number 


Q 

Rate of heat flow 

Bbr-i 

R 

Thermal resistance 

B-ibrF 

R 

Electrical resistance (only used with prime sign or 



subscripts) 

ohm 

r 

Radial distance 

ft 

m 

Reynolds number 


s 

Side length of a square 

ft 


* For lengths, forces, and masses the units of the British gravitational system are 
used. Time is given in hours instead of seconds in accordance with the customary 
engineering practice in heat transfer work. Temperatures are given in degrees 
Fahrenheit (F), absolute temperature in degrees Rankine (R). 
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Symbol Quantity Units* 


T 

Absolute temperature 

R 

T 

Fundamental unit of time 

hr 

t 

Temperature 

F 

U 

Overall heat transfer coefficient 

B hr"“^ ff^ ' 

V 

Velocity 

ft hr'i 

X 

Distance; thickness 

ft 

y 

Distance 

ft 

a 

Thermal diffusivity 

ft2 hr“i 

a 

Absorptivity for radiation 



Coefficient of cubical expansion 

R-i 

7 

Specific weight 

lb fr® 

(AOw 

Log mean temperature difference 

F 

€ 

Emissivity 


€ 

Temperature coefficient of electrical resistance 

F-i 

e 

Fundamental unit of temperature 

F 

d 

Temperature difference 

F 

K 

Temperature coefficient of thermal conductivity 

F-i 

X 

Wavelength 

ft 


Dynamic viscosity 

lb ff”^ hr 

V 

Kinematic viscosity 

ft2 hr-i 

P 

Density 

slug ft""® 

P 

Reflectivity for radiation 


<r 

Constant in Stefan-Boltzmann's law 

B hr~^ fr^ R- 

r 

Time 

hr 

r 

Transmissivity for radiation 



* See footnote page xv. 



CHAPTER I 
INTRODUCTION 

I-l Importance of Heat Transfer and Insulation in Engineering 

The laws which govern heat transmission are very important to the 
engineer in the design, construction, testing, and operation of heat 
exchange apparatus. Heat transfer problems confront investigators in 
nearly every branch of engineering. 

Electrical engineers apply their knowledge of heat transfer to the 
design of cooling systems for niotors, generators, and transformers. 
Chemical engineers are concerned with the evaporation, condensation, 
heating, and cooling of fluids. The mechanical engineer deals with 
problems of heat transfer in the fields of internal combustion engines, 
steam generation, refrigeration, and heating and ventilating. An 
understanding of the laws of the flow of heat is important to the civil 
engineer in the construction of dams and structures, and to the architect 
in the design of buildings. 

In the latter instance the question of effective heat insulation will 
prevail. However, problems involving the use of insulation likewise 
occur in each of the fields of engineering. Actually every engineer will 
be confronted from time to time with the question of how to transmit 
heat in the most effective way, or of how to protect a construction most 
efldciently against heat or cold losses. 

1-2 The Three Kinds of Heat Transfer 

There are three different types of heat transfer: conduction, con- 
vection, and radiation. (Ref. I-l.) They have in common that temper- 
ature differences must exist and that heat is always transferred in the 
direction of decreasing temperature. On the other hand they differ 
entirely in the physical mechanisms and laws by which they are 
governed. 

Heat conduction is due to the property of matter which allows the 
passage of heat energy, even if a physical body is impermeable to any 
kind of rays and its parts are not in motion relative to one another. 

Heat convection is due to the faculty of moving matter to carry heat 
energy such as transporting a load from one place to another. 

1 
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Heat radiation is due to the property of matter to emit and to absorb 
different kinds of rays, and to the fact that an empty space is perfectly 
pei-meable to raj’s and that matter allows them to' pass more or less. 

If the flow of heat is independent of time, it is spoken of as steady 
or stationary state. Engineering problems of steady flow of heat by 
conduction include the flow of heat through furnace walls and pipe 
insulation where the surface temperatures are maintained almost 
constant, and in general where changes of conductive heat flow with 
time are negligibly small. 

The fundamental relation for the steady flow of heat by conduction 
originates from the French physicists Biot and Fourier and can be 
expressed by 

gk = kA^ [I-I] 

This formula is based on the assumption of a homogeneous substance 
in which a constant temperature difference A^t is held between the 
points of a plane area A and any points at a short perpendicular dis- 
tance Ax from this area. Then a steady heat flow per unit time 
(rate of heat flow) occurs in the direction of decreasing temperature. 
The factor of proportionality k is called heat conductivity or “ ther- 
mal conductivity.’^ 

Whenever heat flows by conduction in the transient or unsteady 
state the temperature of a fixed point within the material does not 
remain constant; that is, the temperature within a material undergoing 
cooling or heating varies with the time. Some of the industrial prob- 
lems which involve this particular type of heat transmission are: the 
annealing of castings, the vulcanizing of rubber, and the heating or 
cooling of the walls of buildings, furnaces, and ovens. 

In unidirectional flow of heat by conduction in the transient state, 
in addition to Eq. I-l the following relation for the rate of heat stored 
qst has to be considered : 


Here A^t means the average temperature increase of the substance in 
the small volume A • Ao; in a short time interval At, when p and Cp are 
the density and specific heat, respectively, of that substance. 
f>Cp{A • Ax) Art apparently is the heat energy stored in the volume 
A • Ax in the time At. 

Heat transfer by convection occurs on walls of rooms, on the outside 
of warm and cold pipes, and between the surfaces and fluids of all types 
of heat exchangers. 
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BASIC UNITS 


For this type of heat transmission the following equation which goes 
back to Newton, is in general use: 

qh = hA- At [1-3] 

It simply states that an invariable temperature difference At between 
the area A of a surface and a fluid in contact with it, causes a steady 
heat flow q^. The factor of proportionality defined by this equation, 
is called film coefficient of heat transfer ” or simply “ coefficient of 
heat transfer.^' 

Heat transfer by radiation is important in boiler furnaces, billet- 
reheating furnaces, and other types of heat exchangers. Solar radiation 
plays an important part in the design of heating and ventilating systems. 

The rate of heat, qrj radiated from an area, A, is given by the equation 

qr = [1-4] 

in which T is the absolute temperature, e is a factor depending on the 
kind of surface and on temperature, and is a constant of nature, inde- 
pendent of both surface and temperature. It will be shown that for a 
perfectly black surface 6 = 1. In this case Eq. 1-4 simplifies into 
Stefan-Boltzmann’s law, so called in honor of Stefan who found it 
empirically and Boltzmann who proved it theoretically. 

It is important to note that the majority of industrial problems deal- 
ing with heat exchangers do not involve a single mechanism of heat 
transfer, but a combination of two or more. In a steam condenser the 
transfer of heat is both by ^conduction and convection from the con- 
densing steam to the cooling water. In the furnaces of large steam gen- 
erators heat is transferred by radiation, convection, and conduction. 

In the following chapters the basic equations governing each type 
of heat transfer and cases dealing with a combination of two or more 
of the ways of heat transfer mil be considered. 

1-3 Basic Units 

Consistency of units in numerical calculations is even more important 
in heat transfer than in most other fields of engineering. This is due 
to the complex arrangement and combination of units in equations. 
For instance, many formulas are given as products of power functions 
of the variables involved which often have fractional exponents. There 
are some texts and articles of good standing which contain equations 
and tables of mixed units, the use of which leads to serious errors unless 
in each case careful individual consideration is given to the conversion 
of units. 

As usual in engineering calculations on heat transmission, in this text 
foot and hour are generally used as units of length and time. In 
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addition the use of the units degree Fahrenheit and British thermal unit 
is obvious. The main difficulties arise in the use of the pound force 
and the corresponding unit of mass, the slug. The latter is that mass 
which when acted upon by a force of one pound will be accelerated 
1 ft/sec“. 

Sometimes, however, the poimd mass will be used as a unit. This 
will be denoted by Ib;^? whereas lb without subscript or with the sub- 
script/ means pound force in this text. 

Confusion in the use of force and mass units must be avoided first of 
all in dealing with the density, p, and the specific weight, y. Because 
the numerical value of the former in [Ib^/cu ft] is equal to that of the 
latter in [lb//cu ft], they often are considered as identical. This how- 
ever is erroneous. In the mme system of units they differ by the 
gravitational acceleration p as a factor according to the equation 

7 = P • g [1-5] 

The unit of density in the foot-hour-pound force system is the slug 
per cubic foot. According to Eq. 1-5 the specific weight which is ex- 
pressed in pound force per cubic foot may be converted to the density 
unit of [slug/cu ft] by dividing by 32.16 ft/sec^. So, the density of 
water at 68 F is 62.305/32.16 = 1.937 slug/cu ft. 

Other physical properties in which force or mass is involved are 
the specific heat and the viscosity. The units of these will be discussed 
in Sects, IV-1 and VI-2, respectively. 

The physical centimeter-gram-second system and some other sys- 
tems are used only occasionally to show conversion or in practice 
problems. 

REFERENCE 

I-l. M. Jakob, A Survey of the Science of Heat Transmission,” Research Series 
68, Engineering Experiment Station Bulletiny Purdue University, Vol. XXIII, 
No. 4a, 1939, 
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THERMAL CONDUCTIVITY 

II-l Dimensions of Thermal Conductivily 

In Eq. I-l the subscript k was used to distinguish the heat flow by 
thermal conduction from that by convection and radiation, and x was 
used to distinguish the temperature difference in the direction x from 
the temperature change in time, denoted by subscript t in Eq. 1-2. 

. If only heat flow by conduction independent of time is considered, 
then Eq. I-l may be written without any subscripts: 

q = kA^ [II-l] 

Ax 

The factor of proportionality k in this relation represents a physical 
property of the substance through which heat is conducted and is 
called “ heat conductivity ” or “ thermal conductivity,” as already 
mentioned. 

When physical equations like Eq. II-l are used it is imperative that 
all the constituent parts be expressed in consistent physical units. 
The most satisfactory way to check the physical dimensions in an 
equation is to assume that the numerical magnitude of each term is 
unity and then to consider only the remaining units. The units may 
be operated on as though they were the symbols of an algebraic equation. 
Writing the units in square brackets and expressing the physical dimen- 
sions of k one obtains* from Eq. II-l : 

l[B/hr] = lMl[sqft]^ [II-2] 

or 

[k] = [B hr-i ft-i F-^] [II-3] 

Equation II-3 gives the physical dimensions of the property k in 
British technical units. 

It is important to note that, whereas in Eq. II-2 sq ft and ft were 
used for the area A = 1 and the distance Aa: = 1, respectively, these 

* B means “ British thermal unit ” and F means “ degree Fahrenheit,” according 
to the “ American Standard Abbreviations for Scientific and Engineering Terms,” 
1941. 
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are combined into ft”^ in Eq. II~3. This simplification is due to the 
use of consistent units, namely, sq ft for areas and ft for distances in the 
case considered. 

If another system of units is to be employed, it is only necessary to 
replace the units in Eq. II-3 b}- their equivalents in the other units. 
For example: 

1 B = 252 cal (abbreviation for physical calorief) 

1 hr = 3600 sec 
1 ft = 30.48 cm 

IF = (f) C (abbreviation for degi'ee centigrade) 
and by substitution in Eq. II-3: 

1[B hr-i ft-i F-i] = r 252 cal —A 1 

L 3600 sec 30.48 cm (f )Cj 

= 0.00413 [cal sec""^ cm~^ C““^] 

This is the conversion factor of the thermal conductivity from British 
technical units into the usual physical units. By similar procedure, 
it is possible to find the following conversion factors; 

I B hr”*^ ft"“^ F""^ = 1.488 kcal hr^^ m“^ C”^ (metric technical units) 

= 0.0173 watt cm"“^ C^’^electro-physical units) 

11-2 Factors Which Influence Conductivity 

The numerical values of k for different substances vary from almost 
zero for gases under extreme vacuum conditions to about 7000 
B hr ^ ft ^ F which has been observed for a natural copper crystal at 
the very low temperature of -422 F. The values of k for a substance 
depend on the chemical composition, the physical state and texture, 
and the temperature and pressure. 

From the viewpoint of molecular physics thermal conduction is con- 
sidered as a transportation of heat energy due to molecular "motion, 
whereas heat transmission connected with the ordinary flow of a fluid 
belongs .to thermal convection. 

The molecular transmission of heat is smallest in gases. It may be 
compared with the transportation of energy by billiard balls. In a gas 
space, as on a billiard table, the efficiency of the exchange of energy is 
excellent, but the field is almost eihpty, compared with the size of the 

* Concerning the use of sq ft for the area and inches for the thickness of a plate 
see Sect. III-2. ’ 

I I cal is the heat quantity required for heating 1 gram of water by 1 C, 1 kcal 
(kilo-calorie) is the heat quantity required for heating 1 kilogram of water by 1 C. 
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balls, and since the balls naove about in a zigzag manner the total 
energy exchange for unit space in a given direction is small. Thus the 
heat conductivity of gases is lower than for any other substances. In 
liquids and in electrically insulating solid bodies the heat energy is 
considered to be handed down from places of higher temperature to 
those of lower temperature by elastic oscillations similar to the propa- 
gation of sound. Because for both kinds of substances matter is packed 
much more closely in a given space than for gases, the heat conduction 
of liquids and electrical insulators is much better than that for gases. 
However, the oscillatory transportation of energy is hampered by the 
irregular arrangement of the molecules and atoms in liquids and in 
so-called glassy or amorphous solid substances. In contrast to this, 
the crystalline structure of a substance, comparable with a chain of 
springs, affords a much more, effective propagation of heat energy in 
the direction of an existing temperature decrease, and therefore crys- 
talline substances are much better heat conductors than liquids or 
glassy bodies. In the metallic substances, other important carriers 
of energy enter into the action, which are known as free electrons. 
Experience and theory have shown that the thermal and electrical 
conductivity of pure metals at the same temperature are approximately 
proportional. 

Details of this complex picture, particularly the theoretical consider- 
ation of the influence of temperature, are beyond the scope of this text. 
Only two items of practical significance will be mentioned. In solids 
used for insulating materials, the insulating effect is due mainly to the 
air in the pores which as a gas is a poor conductor of heat. The apparent 
density of a porous or fibrous body which is defined as the ratio of its 
total mass to its total volume,* is small because of the large number of 
air voids. This accounts for the decrease in the heat conductivity of 
insulating materials as the apparent density becomes smaller. If 
water penetrates into a porous substance and fills the pores, for instance 
due to excessive atmospheric humidity, then the body can no longer 
be considered as a combination of solid material and air, but a part of 
it is water with a corresponding higher thermal conductivity value. If 
instead of fine pores, wide air spaces are considered, it might be con- 
cluded that the insulating effect would be the same as for a substance 
having small pores and the same apparent density. This is not true, 
because of air circulation within the large free spaces which decreases 
the insulating effectiveness of the substance. In addition, for the 
same temperature range, more heat is radiated through a substance 

* In the tables of this chapter “ apparent specific weights ” are used instead of 
“ apparent densities.” See the remark in Sect. 1-3. 
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composed of wide spaces than through the same volume with narrow 
spaces. This fact can easily be derived from the laws of radiation. 
Radiation through the pores is generalh” included in the value for ky as 
is convection if the pores are wide enough so that this, too, must be 
considered. 

These theoretical remarks may be sufficient for a general under- 
standing of the large variations in the heat conducthdty values of 
various materials which are presented in the follovdng sections for 
information and use in numerical calculations. The sequence is, as 
far as possible, from smaller to larger values of k. Where the physical 
dimensions are not indicated explicitly, the dimensions for k are 
B hr”^ ff"^ Temperatures are in degrees Fahrenheit. 

n-3 Gases and Vapors 

Table II-l shows the thermal conductivity at 32 F and 212 F for a 
number of gases and vapors along with the molecular weight. The 
table relates to a range in which, according to theory and experience, 
the pressure has no -appreciable influence on the thermal conductivity. 
This is the case for moderate vacuum conditions, and also for approxi- 
mately atmospheric pressure, except for those substances which are 
liquids or not far removed from saturation at this pressure. From the 
table, it is apparent that the thermal conductivity of gases or vapors 
with higher molecular -weight is smaller than that for the lighter ones. 

0.022 

0.020 

Ll. 

r 0.018 

'V. 

QQ 

£ 0.016 

!> 

0 

1 0.014 

o 

0 

1 0.012 

a> 

sz 

H- 

0.010 

0.008 

-150 -100 -50 0 50 100 .150 200 250 300 350 

Temperature, F 

Fig. II-l. Thermal conductivity of some gases and vapors. 
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Figure II-l shows the variation of the thermal conductivity with 
temperature for some gases and vapors. In general, the thermal con- 
ductivity increases with an increase in temperature. 


TABLE II-l 

Thermal Conductivity of Gases and Vapors at Moderate Vacuum 

Pressures 


Gas or Vapor 

1 

Chemical 

Formula 

Molecular 

Weight 

Thermal Ci 
[B hr”^ 
a 

32 F 

onduotivity 
ft-i F-i] 

.t 

212 F 

Hydrogen 

H 2 

2 

0.099 

0.124 

Helium 

He 

4 

0.082 

0.097 

Methane 

CH4 

16 

0.0175 


Ammonia 

NHs 

17 

0.0124 

0.0171 

Water vapor 

H 2 O 

18 


0.013 

Nitrogen ^ 

N 2 

28 

0.0141 

0.0175 

Ethylene 

C 2 H 4 

28 

0.0097 

0.0154 

Air 


29 

0.0141 

0.0177 

Ethane 

C 2 H 6 

30 

0.0104 

0.0184 

Oxygen 

O 2 

32 

0.0141 

^ 0.0181 

Carbon dioxide 

CO 2 

44 

0.0081 

0.0121 

Ethyl alcohol 

C 2 H 5 OH 

46 

0.0081 

0.0121 

Benzene 

CsHe 

78 

0.0050 

0.0101 

Carbon tetrachloride 

CCI 4 

154 


0.0050 


Liquids 


Accurate data are only available for a few liquids 
Fig. II”2 show the thermal 
conductivity of water as a 
function of temperature. It 
is seen that a maximum oc- 
curs at about 250 F. Water 
. is the best conductor of all 
non-metallic liquids. The 
heat conductivity of organic 
liquids at ordinary tempera- 
ture is of the order of mag- 
nitude of 0. 1 B hr“^ ft“^ F~^ 
with variations of seldom 
more than ±20 per cent. 


Table II-2 and 



Fig. II“2; Thermal conductivity of water. 
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TABLE II-2 

Theemal Conductivity op Water for Various TEs^rpEEATUREs 


Temperature 

Thermal Conductivity^ 

F 

B hr-i ft“i F-i 

40 

0.327 

50 

0.333 

60 

0.336 

70 

0.344 

SO 

0.349 

90 

0.354 

100 

0.360 

no 

0.364 

120 

0.366 

130 

0.371 

140 

0.377 

150 

0.381 


taken from A. Eucken and M. Jakob, Der Chemie-Ingenieur, Vol, I, Part 1, Chanter VI 
Akademieche Verlagegesellschaft, Leipzig, 1933. 


n~6 Insulating Materials 

Several of the most used types will be described here in some detail. 
In the selection of a suitable insulation such factors as insulation 
efEciency,t resistance against temperature, shock, or corrosive chemical 
fumes, cost, and reuse value should aU be considered. In general, all 
heat-insulating materials may be grouped into the following general 
classes. 



Fig. II-3. Cork block insulation., 
t Defined in Sect. Xn-3. 
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(a) Low-Temperature Inmlations 

In this class the insulators are used to prevent freezing or sweating 
and to insulate ice-water pipes and other refrigeration equipment. 
The materials are manufactured so as to have a low conductivity and a 
durable structure under moisture conditions. The problem of keeping 
moisture away from the cold surfaces and at the same time keeping it 
outside of the insulation itself is of paramount importance in the design 
of low-temperature insulation. Cattle hair, wool felt, cork, and various 
combinations of these serve extensively as raw materials. Cork is 



Fig. II-4. Pipe covering used to prevent sweating. 


used in many forms such as sheets, blocks, or bulk. Figure II-3 shows 
a section of cork block. In Fig. II-4 is to be seen a combined low- 
temperature pipe insulator designed to prevent sweating and to protect 
against moisture. 

(6) Insulation for Building Purposes 

This type of insulation is used primarily for reducing the heat loss, 
from building walls, roofs, etc., during cold weather and the heat pene- 
tration during warm weather. Aside from cork, many special insu- 
lations as rock wool, slag wool, and glass wool are used. The latter 
three consist of fine fibers which resemble the general appearance of 
raw cotton. The fibers are made by blowing the molten material into 
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Fig. II-6. Glass wool. 
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thread form. These insulators are applied in bulk form, blankets, and 
special reinforced batts. Samples of rock wool and glass wool are 
sho\\m in Figs. II-5 and 6 respectively. 

Metallic or reflecting type insulation is also used in this field. It 
consists essentially of very thin metal foils which, for instance, can be 
fastened to building paper. The insulating effect of such metallic sur- 
faces will be explained later. (See Sect. XII-4.) 



Fig. II-7. Asbestos paper pipe insulation. 


(c) Insulation for Heating and Process Work 

Owing to the relatively low temperatures encountered in this field 
(150 to 300 F) asbestos paper structures are extensively used. They 
have the advantages of a rather small conductivity and low cost. 
Asbestos paper is not suited to much higher temperatures because of 
the decomposition of the binder. Figure II~7 illustrates a type of 
pipe insulation made up of layers of asbestos paper having thin 
air' spaces between layers. The structure is to be seen clearly in the 
figure. Figure II-8 shows a section of corrugated asbestos paper 
insulation board. The same type is also made in the form of pipe 
insulation. 
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Fig. II-8. Corrugated asbestos paper board insulation. 


(d) Insulation in the Power Generation Field 

In this field the temperatures encountered are usually between 300 
and GOOF. Because the resistance of the insulation to mechanical 

vibration and shock must be high, 
an insulation composed of 85% 
magnesia and 15% asbestos called 
‘‘85% magnesia” is often used. 
It is very effective and has high 
reuse value. Asbestos sponge felt 
and many other materials are like- 
wise on the market. Figure II-9 
shows a piece of “ 85% magnesia ” 
pipe insulation. Figure II-9a gives 
the thermal conductivity of the ma- 
terial independent of temperature. 

(e) High-Temperature Insulations 
Diatomaceous earth and asbestos 
are used primarily in this range 
(600 to 1900 F) because of their re- 
sistance to decomposition and high 
temperatures. Diatomaceous earth is excellent in insulating value, 
temperature resistance, strength, and durability. It consists of tiny 
silica skeletons of “diatoms,” microscopic plants which have been 
deposited millions of years ago. Layers of this material have been 



Fig. II-9. Pipe insulation made of 85% 
magnesia and 15% asbestos. 
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Fig. II-9a. Thermal conductivity of an insulation consisting of 85% magnesia and 
15% asbestos having an apparent specific weight of 16.9 Ib/ft®. 

Data from R. H. Heilmaa, Industrial and Engineering Chemistry, 28, 7S2 (1936). 


found which are 1400 ft thick. It is used either in natural form 
(Celite) or in the form of insulating bricks which are made by grinding 
the earth, pugging, pressing, and firing it in kilns. An insulating brick 
cut from the material as found in nature is shown in Fig. II-IO. Saw 
marks on the sides and layer marks on the top surface are visible. 

Asbestos is a silicate mineral having the characteristic properties of 
fibrous structure and fire resistance. Diatomaceous earth and asbestos 
are often mixed together to form effective insulations for this tempera- 
ture range. 



Fig. II-IO. Insulating brick cut from natural Celite. 
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An insulation consisting of one layer of diatomaceous earth and one 
layer of '' 85% magnesia ” serves for high-temperature pipe coverings. 
In the combined form the high-temperature insulation is used next to 
the hot surface and in sufficient thickness to avoid temperatures of more 

than 600 F on the outer sur- 
face. ‘ A layer of “ 85% mag- 
nesia is then placed around 
the high-temperature insula- 
tion to decrease further the 
heat loss. Figure II-l 1 shows 
a piece of combined high- 
temperature and “ 85% mag- 
nesia insulations for use on 
high-temperature pipes. 

Many other types of insu- 
lations are on the market, 
but a complete discussion is 
beyond the scope of this text. 

In Table II~3 are given 
data relative to the conduc- 
tivity of some insulating ma- 
terials together with the 
Fxg. II-ll. Combined pipe insulation. apparent specific weight and' 



TABLE II-3* 


Thermal Condtjctivitt of Some Insulating Materials 


Insulation 

Apparent 
Specific 
Weight 
Ib/cu ft 

Thermal Conductivity 
[B hr-i ft-i F-i] 
at 

20 F 1 100 F ] 200F 

Corrugated asbestos (4 plies 
per inch) 

13.0 

0.0354 

0.0430 

0.0525 

Glass wool 

1.5 

0.0217 

0.0313 

0.0435 


4.0 

0.0179 

0.0239 

0.0317 


6.0 

0.0163 

0.0218 

0.0288 

Hair felt 

8.2 

0.0237 

0.0269 

0.0310 


11.4 

0.0212 

0.0254 

0.0299 


12.8 

0.0233 

0.0262 

0.0295 

85% magnesia ” 

16.9 

0.0321 

0.0344 

0.0373 

Rock wool 

4.0 

0.0150 

0.0224 

0.0317 


8.0 

0.0171 

0.0228 

0.0299 


12.0 

0.0183 

0.0226 

0.0281 


* Data taken from R. H. Heilman, Ind. and Eng. Chemistry, 28 , 782 (1936). (See als6 Fig. II-9o.) 



II-6] REFRACTORY MATERIALS 17 

mean temperature. The thermal conductivity of powdered diatoma- 
ceous earth for various apparent specific weights can be taken from 
Table II-4. Examination of the data in these tables indicates that 
the thermal conductivity of insulating materials increases with tempera- 
ture and apparent specific weight. The increase due to temperature 
is greater for the lower apparent specific weights. 


TABLE II-4* 

Thermal Conductivity of Powdered Diatomaceous Earth 


Apparent 
Specific 
Weight 
Ib/cu ft 

100 F 

200 F 

Thermal C( 
[B hr-i j 
a 

300 F 

inductivity 
Ft-' F-'] 

.t 

400 F 

500 F 

600 F 

10 

0.024 

0.029 

0.034 

0.038 

0.044 

0,048 

14 

0.030 

0.033 

0.036 

0.039 

0.043 

0.046 

18 j 

0.038 

0.040 

0.043 

0.045 

0.048 

0.049 


* Data taken from Diatomaceous Earth by R. Calvert, Chemical Catalog Co., Inc., 1930. 


The size and number of the spaces within the insulating material as 
well as the arrangement of the fibers in fibrous materials have a marked 
influence on the thermal conductivity. 

II-6 Refractory Materials 

In the design and construction of furnaces and related apparatus 
refractory materials are used owing to their ability to withstand high 
temperatures without serious deterioration. It is sometimes desirable 
to select refractory materials of low thermal conductivity, whereas in 
other cases, good conductors of heat are required. Table II-5 shows 
the thermal conductivity of twelve different types. 

From the data on kaolin bricks in this table it is apparent that for a 
given temperature the thermal conductivity decreases as the porosity 
increases. This comes from the insulating effect of small air spaces 
which has already been mentioned. ' 

The thermal conductivity of most of the refractories quoted in 
Table II-5 increases with temperature. The opposite is true with the 
magnesite brick and the silicon carbide brick. The reason for this 
different behavior must be sought in the inner structure of the re- 
fractory materials (Ref. II-l). It has been found, theoretically as 
well as by experiments, that the thermal conductivity of crystals varies 
as the reciprocal of the absolute temperature, whereas amorphous 
(glassy) substances show an increase of the conductivity with the 
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*Data taken from F. H. Norton, J. Am. Ceram. Soc., 10 , 30 (1927). 

t The values in parentheses are round means of values, given in " Modern Refractory Practice ’’ by the Harbison-Walker Refractories Co. 
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temperature. So in a complex refractory material the crystalline com- 
ponents act in the one direction, the glassy components, the air in the 
pores, and the radiation through the pores act in the opposite direction 
as concerns the influence of temperature on the apparent heat con- 
ductivity. The two mentioned bricks apparently consisted mainly of 
crystalline constituents. The great value of the thermal conductivity 
of the magnesite brick at 392 F is a consequence of this behavior. For 
the same temperature, heat conductivities of magnesite bricks up to 
> 3, and for silicon carbide bricks up to ^ > 11 have been observed 
by reliable investigators. 



Specific Weight, Ib/cu ft 


Fig. 11-12. Thermal conductivity of Douglas fir at a mean temperature of 75 F. 

Data taken from “ Thermal Conductivity of Building Materials,” by F. B. Rowley and A. B. 
Algren, Bulletin 12, Engineering Experiment Station, University of Minnesota, 1937. 


II-7 Building Materials 

The thermal conductivity of wood depends upon such factors as 
moisture content, density, and temperature. Figure 11-12 shows the 
variation of the thermal conductivity of Douglas fir for various specific 
weights and moisture contents at a mean temperature of 75 F. From 
the figure it is seen that for a given moisture content the thermal con- 
ductivity varies linearly with the specific weight and that it increases 
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with the moisture content for a given specific weight. Figure II-13 
shows the thermal conductivity of various woods of different specific 
weights at a constant moisture content of 12 per cent. A straight line 
represents the results fairly well* 
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jz 

OQ 

’I 0.08 


£ 0.07 
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0.05 


0.04 

20 25 30 35 40 45 50 

Specific Weight, !b/cuft 

Fig. 11-13. Relation between thermal conductivity and average specific weight 
for twenty different woods at a mean temperature of 75 F and 12 per cent moisture 

content. 

Data taken from " Thermal Conductivity of Building Materials,” by F. B. Rowley and A. B. 
Algren, Bulletin 12, Engineering Experiment Station, University of Minnesota, 1937. 
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Thermal conductivities for various other materials are given in 
Table II-6. Again the strong influence of the moisture content is 
striking. 

* The most recent results of experiments of J. D. MacLean on various woods 
[see Mech. Eng. 63, 734 (1941)] came to our knowledge too late for use in preparing 
this text. 
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TABLE II-6 

Thermal Conductivity for Various Building Materials 


Material 

Apparent 
Specific Weight 
Ib/cu ft 

Thermal 
Conductivity 
Bhr“^ft-^ F-i 

Asphalt 

130 

0.40 

Brick, dry 

105 

0.30 

1% moisture by volume 

105 

0.40 

2% moisture by volume 

105 

0.60 

Concrete, dry 

120 

0.45 

10% moisture by volume 

140 

0.70 

reinforced 

140 

0.75 

Glass 


0.5 to 0.6 

Gypsum, dry 

80 

0.25 

Limestone, fine grain, dry 

105 

0.40 

fine grain, 15% moisture by volume 

105 

0.55 

coarse grain, dry 

125 

0.55 

Linoleum 

75 

0.11 

Plastic, with 2% moisture *by volume 
Rubber, vulcanized, soft, with 40% pure 

115 

0.5to0.6 

rubber 


0.17 

with 90% pure rubber 


0.10 

Sand, dry 

95 

0.20 

10% moisture by volume 

100 

0.60 

Soil, dry 

125 

0.30 

fresh clay, 28% moisture by volume 

125 

1.35 


^-8 Metals and Alloys 

Pure metals are the best heat conductors. Because of their crys- 
talline structure, the thermal conductivity increases inversely to the 
absolute temperature as with n9n-metallic crystals. Impurities, how- 
ever, disturb the free conduction; they have not only the effect of 
reducing the heat conductivity appreciably, but also they likewise 
may change the decrease of the conductivity to an increase with tempera- , 
ture in the higher temperature range. Alloys generally behave like 
amorphous substances having relatively small thermal conductivities 
which increase with temperature; but a decrease is sometimes observed 
as well. 

Table II-7 gives values of the thermal conductivity of entirely or 
rather pure metals. The mentioned influence of temperature and the 
effect of even small impurities will be noticed. 

Values for a few non-iroii alloys are given in Table 11-8. Here- the 
increase of the conductivity with the temperature is general, and it is 
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TABLE II-7 


Thermal Condtjctivttt of Different Metals 


Metal 

Impurities 

% 

-58 F 

Thei 

[I 

32 F 

:mal Conduc 
i hr-i ft-i I 
at 

392 F 

jtivity 

752 F 

1112 F 

Aluminum 

0.3 

134 

132 

131 

131 



1.0 

116 

115 

114 



Copper 

0.1 

» . . 

220 

215 

210 

206 

Iron 

0.1 


43 : 

36 

28 

22 

Lead 

0.0 

22 

21 

19 



Nickel 

0.1 ! 



42 

34 



1.0 

34 

35 

34 

34 

34 

Platinum 

0.05 

41 

40 

42 

45 

47 

Silver 

0.0 

243 

242 





0.1 

236 

233 

217 

204 

225 

Zinc 

0.2 

67 

66 

62 

54 



TABLE II-8 


Thermal Conductivity of. Some Non-iron Alloys 




Thermal Conductivity 


Material 


[B hr-i ft-i F-i] 
at 



-58 F ! 

32 F 

392 F 

752 F 

Brass: 70% Cu, 30% Zn 

55 

60 

75 

80 

“ChromelA”; 80% Ni, 20% Cr 

i 


9 

11 

“ Chromel P ” : 90% Ni, 10% Cr 
Gunbronze: 86% Cu, 9% Sn, 4% 



12 

14 

Zn 



i 38 

44 

Monel; 67% Ni, 29% Cu 



18 

20 


further seen that some alloys conduct the heat less than any of their 
constituents. 

Table II-9 contains values for some iron alloys, all of them heat- 
treated. It is seen that increasing chromium and nickel contents re- 
duce the conductivity to about one-half or one-third the value for 
wrought iron. The decrease of the thermal conductivity with increasing 
temperature turns to an increase only for the high-grade chromium and 
nickel steels. 

The range in values of the thermal conductivity for various gases, 
vapors, liquids, and solids is represented in Fig. 11-14. 
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THERMAL CONDUCTIVITY [II-8 



Absolute Temperature, R 

Fig. II-14. Thermal conductivity of gases, liquids, and solids. 
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CHAPTER III 

CONDUCTION OF HEAT IN THE STEADY STATE 


in-1 Definition of Steady State 

A steady state transmission of heat by conduction exists whenever 
the flow of thermal energy during equal intervals of time is constant, 
and the temperatures at various positions within the material remain 
fixed. In order to apply the basic law (Eq. I-l) to different specific 
cases, it may be expressed in the following differential form: 


The negative sign indicates that the temperature decreases as the 
distance from a reference point in the direction of the heat flow increases. 

This formula may also be used if the area and thermal conductivity 
both vary with x. By integrating the equation, there results a number 
of equations which can be used for solving problems which occur in 
engineering practice. In this chap- 
ter only cases where the surface 
temperatures of the bodies under 
consideration are known will be dis- 
cussed. 

III-2 ConductionThrough a Homo- 
geneous Plane Wall 

It will be assumed that the ma- 
terial composing the wall is homo- 
geneous and that the thermal con- Fig. III-l. Homogeneous plane wall, 
ductivity of the wall material is 

independent of the temperature. The direction of the heat flow together 
with the two surface temperatures is indicated in Fig. III-l. The first | ^ ; 
step necessary in the analysis of the problem is the separation of the] S 
variables in Eq. III-l. This leads to 

, -kA'dt 
ax == 
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CONDUCTION OF HEAT IN THE STEADY STATE 


Integrating this relation gives 


—kAt 

X — r o 

5 


where C represents a constant of integration which may be determined 
from the following boundary conditions (see Fig. III-l). 

When X == xi, t = h and when x = X2, t = h. Substituting these 
conditions in the integrated relation gives the following equations; 


Xi 



+ C 


+ C 


The constant C is eliminated by subtracting the lower equation from 
the upper: 

kAt2 

q q 

Solving this equation for q gives a relation useful for calculating the 
heat flowing through the wall : 

q = kA “ [III- 2 ] 

In this equation (ti — ^2) and (x2 — xi) represent the temperature 
difference and the thickness of the wall respectively, and A and k refer 
to the area and thermal conductivity , of the material. Introducing 
At = ti — t2 and Ax = X2 — xij Eq. III -2 takes the form from which 
the development started, namely 

7 A 

,~U- m-ii 

but now without the restriction that Ax and At be small. 

If 'one expresses q in B/hr,* A in sq ft. At in F,* and Ax in ft, then in 
order to have consistent units, k must be in B hr~^ ft”^ F^^ The unit 
ft ^ comes from A/ Ax, which indicates the necessity to convert the 
thickness Ax into foot units when it is given in inch units. 

In practice k is often expressed in 

B/sq ft 
r — — ~ units 
hr • F/m. 


* See footnote, p. 5. 
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2 


where the cross-sectional area is expressed in sq ft and the thickness i 
inch units. This practice is objectionable in that it leads to compl: 
cations and errors when bodies other than plane plates are to be deal 
with. 

Example III-l. The interior of an oven is maintained at a temperatui 
of 1500 F by means of suitable control apparatus. If the walls of the ove 
are 9 in. thick and constructed from a material having a thermal conductivit 
of 0.18 B hr“^ ft“^ F“h calculate the heat loss for each square foot of wa 
surface per hour. Assume that the inside and outside wall temperature 
are 1500 F and 400 F respectively. 

Solution: Equation III-2 applies here. Substituting the necessary value 
and dimensions, particularly Ax = K 2 ft, in the equation, gives 


= 0.18(1) 


1500 - 400 

q/iq 


= 9«A n /L 


^ q/A = 264 B hr'^ ff 


This answer indicates that 264 British thermal units will flow through eacl 
square foot of wall area pe 


hour. 

III-3 Conduction Throng] 
a Composite Plane Wall 

Equation III-2 will no\v 
be used to establish a re- 
lation for the heat flow 
through a composite wal 
shown in Fig. III-2. As n 
the previous case, it wil 



Firebrick 

ki2 

Insulat- 
ing 
‘ brick 
kzz 

Red 

brick 

^34 







h 

tz\ 

h 



be assumed that the ther III-2. Composite plane wall 

mal conductivities for the 

various materials composing the wall are independent of the 
peratures. The heat flow through an area A is the same for 
section and may be represented as: 


tern- 

each 


ki2A{ti — tz) 


for the firebrick 


^2zA{t2 ^ is) 
(Ax)23 


for the insulating brick 


^34-d (tz ^^4) 
(Aa:)34 


for the red brick 
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Solving for the various temperature differences: 

fciaA- 


$(Aa:)23 


q{£^)zi 

t3 ““ 4 — : 

k^4,A 


Then by adding these three temperature differences the following is 
obtained: 



q r (Aa;)i 2 
A _ ^12' 


(A3:)23 (Aa:)34 ~ 

^23 ^34 - 


or 


A(ti-U} 

^ (Ax) 12 ^ (Aa:)23 ^ (Aa;)34 

ki2 ^23 ^34 


[III-3] 


This equation may be generalized for any number of layers. Calling 
the total temperature difference 2^ (At) and using the summation sign 
53 also in the denominator, gives 


AJLiAt) 



Example III-2. The wall of a furnace is made up of 9 in. of firebrick, 
5 in. of insulating brick, and 7.5 in. of red brick. If the average values of 
the thermal conductivities for the firebrick, insulating brick, and red brick 
are 0.72, 0.08, and 0.5 B hr~^ ft“^ F“^ res^pectively, calculate the temperatures 
h and is at the contact surfaces of the firebrick and insula^g brick, and of 
the insulating brick and red brick respectively. Assumemat the inner and 
outer surface temperatures h and 4 of the wall are 1500 F and 150 F respec- 
tively, and that the resistances of the mortar joints are negligible. 

Solution: The rate of heat flow-through the wall is obtained from Eq. 
III-3. 


1(1500 - 150) 

9/12 5/12 7.5/12 

0.72 0.08 0.5 


- 180 B/hr 


Since this value of q is equal to the heat flow through any section, Eq. 
III-2 will be used to calculate the temperature 4 by applying it to the fire- 
brick section. 

q = 0.72(12/9) (1500 - h) 
which yields h = 1312 F. In like manner 

- 375 F 
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:-4 Heat Resistance 

The heat resistance or thermal resistance, is defined as the ratio 
he length of the heat flow path to the product of the thermal con- 
ivity and the area. Units of the heat resistance are hr F. 
rms of this concept, Eq. III-2 may be written in the following 
rxer: 


ti — t2 

^ “ R 

E 

lis relation is similar to Ohm’s law, I = — . The analogous terms in 

Re 

3 two equations are electrical potential {E) and thermal potential 
— the electrical resistance {R^) and heat resistance (E), and 
3 electrical current (I) and heat flow (g). Equation III-3 reduces to 
3 following if the resistance term R is used: 


<1 = 


h ^4 

i?12 + ^23 + -B 34 


[III~3a] 


cordingly, the solution of Ex. III-2 on this basis becomes 


1500 - 150 
1.04 + 5.21 + 1.25 


= 180 B/hx 


“5 Conduction Through a Homogeneous Cylinder Wall 

Jor this case a cylindrical section of insulation surrounding the pipe 
3wn in Fig. III-3 will be used in the analysis. It will again be assumed 
t the" insulating material is homogeneous 
that the thermal conductivity of the 
jrial is independent of temperature. The 
of the heat which flows radially through 
.ifferential volume 2x7* • dr * L is 

q=^ -k-2wr--^Ly ' [III-4] 

X/ d/T 

lis relation r represents the variable ra- 
of the pipe and L the length of the sec- 
For simplification, a linear foot of pipe 
be selected so that L equals unity. In 
ler to keep the equation consistent in units, g' = q/L in B hr”“^ ft~^ 
ntroduc^H ‘RmiRtimi TTT-4 will now he nrenared for inteeration bv 



Fig. III-3. Homogeneous 
cylinder wall. 
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grouping the terms as follows: 

r 

Integrating, ^ 

t= - -^\nr + C [III-5I 

Zirk 

The constant of integration C vail be determined from the following 
boundary conditions (see Fig. III-3). 

When r = ri, t = h and when r = t = ^ 2 - Substituting this in 
Eq. III-5 one obtains 

h = ■- In ri + C 

C will be eliminated by subtracting the lower equation from the upper: 


^2 


Solving this relation for gives the equation for calculating the heat 
passing radially through a cylindrical section of pipe insulation 1 ft in 
length : 




/ 


2Tk{ti — ^2) 
In To/r 


[III-6] 


Example III-3. Calculate the heat loss from 30 ft .of 2-in. nominal pipe 
covered with 1 in. of an insulating material having an average th^mal con- 
ductivity of 0.0375 B hr"’- ft~^ F”h Assume that the inner and outer sur- 
face temperatures of the insulation are 380 F and 80 F respectively. 

Solution: The outside diameter of a 2-in. nominal pipe is 2.375 in. Add- 
ing 2 in. to the pipe diameter, the insulation diameter is ^qual to 4.375 in. 
Applying Eq. III-6, the heat loss per linear foot is equal to 

2x0.0375(380 - 80) 

^ 4.375/12 
“ 2.375/12 


The total heat loss for 30 ft of pipe is equal to (116) (30) = 3480 B/hr, / 

ni-6 Conduction Through a Composite Cylinder Wall 

Equation III-6 will be used to establish a relation for calculating the 
flow of heat per linear foot through the insulation shown in Fig. III-4. 
The heat flow through each section is the same in the steady state and 



III-7] 


PEOBLEMS 


33 


By subtracting, C is eliminated, and by rearranging the following is 
obtained: 


^ ko 1 + “ (ii + fe) ^ = 

2 JX2 — Xi 


h ~~ t2 

^X2 Xt 


[in-9] 



of the thermal conductivity. Therefore, if a linear function exists be- 
tween the thermal conductivity and the temperature as given by 
Eq. III-8, the value of A: in Eq. III-2 may be replaced by a value of k 
taken at the mean temperature. It has been shown (see p. 40 of 
Ref. I-l) that the same holds for Eq. III~6 and also for the case of 
a sphere (see Problem III-16). 


Example III-5. Calculate the heat loss through 1 sq ft of a flat slab of 
85% magnesia 15% asbestos insulation 2 in. thick if the surface temperatures 
are 260 F and 180 F respectively. 

Solution: For a mean temperature of 220 F, Fig. II-9a yields 

_ - 0.0377 B hr-i ft”^ F-^ ^ 

and by use of Eq. 111-9 the heat loss becomes 

q” = 0.0377 — ° = 18.1 B hr-i ft"® 

2/12 

or 

q = Ag" = 18.1 B/hr 


PROBLEMS 

III-l. If the heat loss per hour through 1 sq ft of a furnace wall 18 in. thick is 
520 BY^determine the outside surface temperature of the wall material. Assume an 
inside surface temperature of 1900 F and an average thermal conductivity of 0.61 
B hr-i ft"^Sr ^ 

III-2. Ttie wall of a cold-storage room is mad^ up of 12 in. of cork having a mean 
thefEBcfal conductivity of 0.028 B hr“^ ft~^ If the inside and outside surface 
temperatures are 20 F and 70 F respectively, calculate the temperature at a plane 
8 in. from the cold inner surface. 

III~3. Calculate the heat loss from a composite furnace wall made up of 9 in. of 
firebrick, 6 in. of insulating brick, and 4 in. of red brick. The inside and outside 
surface temperatures of the wall are 2000 F and 110 F respectively. The mean 
thermal conductivity values for the firebrick, insulating brick, and red brick are 
0.7, 0.08, and 1.0 B hr”^ ft“^ F~^ respectively. 

III-4. Verify the following: In order to convert in B hr”^ ft*"^ F"^ to 
watt cm~^ C~^ multiply by 0.0173. 

T> / „ ~ fx 

III-5. Convert a value for Jc of 60 , — hito B hr""^ ft“^ F~^. 

hr F/m. 



CHAPTER IV 

CONDUCTIOK OF HEAT IN THE UNSTEADY STATE 

[V-l General Equations for Unsteady State 

If the rate of heat flow and the temperature at any point of a system 
change with time, the condition is termed “ unsteady state ” or “ tran- 
sient state.” Whereas until now only heat conduction in the steady 
state has been considered, this chapter deals with some of the most 
mportant cases of heat conduction, under unsteady state conditions. 
Many industrial problems of heat transmission eome^iijder this class 
Df heat flow. One familiar example is the flow of heat through a building 
(vall during the daily twenty-four-hour heating and cooling cycle. 
Other examples include the annealing of castings, cooling of ingots, 
Duming of bricks, and the heating and cooling of slabs and furnace wall 
sections. 



x=o x=L 


Pig. IV-1. Front view of Fig. IV-2. Cross section of plate with 
plate. incoming and outgoing heat flow. 

The general equation for the flow of heat by conduction in the un- 
steady state will be derived for the simplest case, the warming up of a 
plane plate which is heated from one side and cooled from the other. 

Figures IV-1 and 2 show a front view and a cross section of the plate. 
Its area may be A and its thickness L. It is assumed that the heat 
flow be uniformly distributed over the area A and that the edges of the 
plate, given by abed in the front view and by oe and hj in the sectional 
view, are perfectly insulated against heat losses or gains. Then heat 
will flow only in the direction x. 
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The rate of flow of the heat which enters at the left surface may be 
called go? that which leaves at the right both being unvariable 
with time, and go > gn. 

At a certain instant of time r, the distribution of temperature t across 
the plate may be given by the curve a in Fig. IV-3a. The rate of heat 
flow qx at any distance x from the left-hand surface can be found from 
this curve by means of Eq. III~1 which in the present case will be 
written as 

q. = -kA^ [IV-1] 


Here the partial derivative dt/dx has been chosen instead of dt/dx in 
order to indicate that the temperature t does not only vary with the 
distance x, but likewise also with the time r. The latter change will be 
indicated by dt/dr. 

The slope of curve a is dt/dx = tan 4>, It is often called the temper- 
ature gradient.” If t decreases when x increases, dt/dx is negative. 
Then, according to Eq. IV-1, g^, becomes positive. This is the case of 
heat conduction in the direction of increasing x. 

'For X = 0 and x ^ L the slope dt/ dx can be calculated by substituting 
the given constant values go and qi for qx in Eq. IV-1. Using the sub- 
scripts 0 and L for the gradient at x == 0 and x == L, this equation 3delds 


and 



?0 

kA 

1 [IV-2] 

1 

qL 

kA 

^ IIV-3] 


Since it was assumed that qo > Ql, it follows that curve a will be steeper 
on the left than on the right surface, as shown in Fig, IV-3a- 

By taking tan <t> from this figure for different distances x one finds the 
distribution of the temperature gradient dt/dx across the plate. It is 
represented by the curve b in Fig. IV-36. Since t decreases from left 
to right, the ordinates of curve h are negative throughout. According to 
Eq. IV-1 each ordinate of curve 6 must be multiplied by (— M) in 
order to obtain the heat flow qx as a function of x. Obviously, all 
values of qx become positive in the present case and because curve h is 
rising from left to right, the heat flow qx will decrease from left to right. 
(Curve c in Fig. IV-3c.) 

The heat flow, however, according to the law of conservation of 
energy can decrease only when heat is to be stored in the wall. To find 
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Fig. IV-3a. Temperature distribution (Ex. 
IV-1). 


0 — X=L 


.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 


Fig. IV-3b. Temperature gradient (Ex. 
IV-1). 



80 

60 





0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 


Fig. IV-3d. Sbpe of temperature gradient 
(Ex. IV-1). 
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this storage; a piece of the wall of differential length dx may be con- 
sidered. 

From Fig. IV-Sc it is seen that 


^2 = ^ 1 + (tan ^pl)dx = gi + 



[IV-^] 


where the subscript 1 relates to the place xi and 2 to the place X 2 = 
xi + dx. 

From Eq. IV-1 it follows by differentiation that 

Bx dx^ 


When this is substituted in Eq. IV-4 and the subscript of the derivative 
is omitted 


dH 

^2 = gi ~ M ^ dx 


[IV-5] 


or 


dH 

dqz = qi — q 2 = hA o dx 
ox 


[IV-6] 


Because qi is the rate of heat flow entering the differential piece at 
xi and q 2 that leaving at X 2 , their difference; which is called dq^ in 
Eq. IV-6, is just the heat energy stored in imit time in a slice of the 
plate of thickness dx. 

As t, dt/dXj and qx, likewise dqs can be represented graphically. It 
is only necessary to determine graphically in Fig, IV-36 the slope 

d{dt/dx) dH 


for different values of x and then use these values of tan x as ordinates 
of a curve d in Fig. IV-3d. If the plate is considered as cut in slices of 
equal thickness dx, the factor (kA • dx) in TEq. IV~6 is constant. It 
follows that the rate of heat energy stored in the slice dx is found from 
curve d simply by multiplying each ordinate by the constant factor 
QcA • dx). 

The rate of heat storage dqz can also be found in a different way. 
Multiplying the volume {A^dx) by the density p gives the mass^ 
Multiplying the mass {A • dx • p)hj the specific heat Cp gives the heat 
stored if the temperature increases by one degree; so if it increases by 
dt in the time dr, the heat stored is 

dt 

dqs = A - dx • pCp. “ 
dr 


fIV-71 
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Equating dqz, as given 


or 


where 


by Eqs. IV-6 and 7, leads to 


A-dx- pCp = kA 

dT dx^ 

[IV-S] 

dt k dH dH 

dr pCp dx^ . ^ dx^ 

[IV-9] 

k 

[IV-9a] 

a = — 

pCp 


is a property of the wall material which is called “ thermal diffusivity.” 

By contemplating curve d in the light of Eq. IV-9, it is seen that by 
multiplying the ordinate dHldx^ by the magnitude a the change of 
temperature in time dt/ dr is given directly at any distance x and at the 
fimp. T. In the present case dt/dr is positive. This means that the 
temperature all over the cross section of the plate increases with time, 
which is obvious because it was assumed that more heat enters than 
leaves the plate. 

The physical dimensions of the thermal diffusivity a can be derived 
from 

“ [k] = [B hr-^ ft"’- F-i] 


[p] = [lb„»ft-3] 


[cp] = [Blb„-il^i] 

Therefore 

[«] = [— 1 = [ft"hr-i] 


The remarkable feature of this relation is that neither a heat unit nor 
a temperature unit appears. 

Further it is seen that the product [pCp] does not include mass or 
force units, but is expressed by [B ft"*'^ F”^]. This, of course, holds 
also if p is expressed in [slug/cu ft] and Cp in [B slug"”^ F“^]. With 
these units the numerical value of p would be only 1/32.2 of the usual 
value of density, as mentioned in Sect. 1-3, and the numerical value of 
Cp would be 32.2 times as large as the usual value. For water, for 
instance, p = 1.937 slug/cu ft has been found (see Sect. 1-3), and the 
specific heat of water is 32.2 B slug~^ F“^ instead of the familiar value 
1. To bring them back to the usual values one can introduce the spe- 
cific weight 7 instead of the density p. Then, according to Eq. 1-5, 
7 = Pfif = (1.937) (32.2) = 62.3 lb//cu ft and pCp = ycp/g. This means 
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that using 7 instead of p requires the use of c^jg instead of Cp in the 
product pCp. For water the numerical value of Cp/p = and the 

physical dimension of this magnitude becomes — = 

[ft sec i 

[B Ih/"^ F“"^]. It is seen that in the pound force system the latter 
expression is not the imit of the specific heat as often is believed and 
claimed, but the unit of Cp/p. Thus, in this system Eq. IV-9a becomes 


k 


y(cp/g) 


[IV-96] 


Dealing with the product pcp- only, and not with its single factors, 
the whole difficulty can be avoided by using a special symbol Cp == pCp = 
7 (cp/^) with the dimension [B ft“"^ that is a specific heat at 
constant pressure related to unit volume instead of unit mass. Then 
Eqs. IV-9a and b assume the common form 



[IV-~9c] 


Example IV-1. The distribution of the temperature t across a large plane 
dry concrete wall 9 in. thick, which is heated from one side, is measured by 
thermocouples inserted in holes in the wall, and it has been found that at a 
certain instant r, the temperature can be represented approximately by the 
equation 

i = 90 - SOx + 16a;2 + 320:^ - [IV-IO; 

where lengths are in feet and temperatures in F. 

If the area of the wall is 5 sq ft, calculate (a) the heat entering and leaving 
in unit time, (6) the heat energy stored in the wall in unit time, and (c) the 
temperature changes per unit of time. 

Solution: The specific weight 7 of concrete is assumed to be 136 Ib/ft®, the 
specific heat Cp = 0.202 B = 0.202(32.2) B slug“^ F”^ and the therma 

conductivity k = 0.44 B hr-^ F-h Then Cp = (136/32.2) (0.202) (32.2; 
= (136) (0.202), and the diffusivity a = k/Cp = 0.016 ft^ hr“h The thick- 
ness of the plate must be expressed in feet, 

L = ^ ft = 0.75 ft 

By substituting different numerical values for x (from a; = 0toa; = L = 
0.75) it will be seen that the temperature distribution is that represented ii 
Fig. IV“3a. Also Figs. IV-36 to d are drawn to scale for the present example 

Differentiation of Eq. IV-10 leads to 

^ = -80 + 32a; + 96a:2 - 102.4a:® 


[IV-11 
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and 

AU 

— - 32 + 192a; - 307.2a;2 [IV-121 

dx^ 

From Fig. IV-3d it is seen that dH/dx^ has a maximum value. According 
to the rules of calculus the distance x at which it occurs is found by forming 
the derivative of the last equation and equating it to zero. Thus 192 ~ 
614.4 X = 0 and x - 0.312 ft. At this point dH/dx^ = 32 + 192(0.312) — 
307.2(0.3122) = 62.1 F/ft^. 

(а) According to Eqs. IV~2 and 3 the values of {dt/dx)Q and idt/dx)^ are 
needed for the calculation of qa and gz,. Substituting a; = 0 and a; = L = 
0.75 in Eq. IV-11 

= IsOF/ft 

\dx/o 

= -80 + 24 + 54 - 43.2 = -45.2 F/ft 

Xdx/L 

and 

50 = -kA = -(0.44)5(-80) = 176 B/hr 
\oxJq 

qL = -^(^1 = -(0.44)5(45.2) = 99.44 B/hr 
\dxj L 

(б) The heat stored is 

ga - go - = 176 - 99.44 - 76.56 B/hr 


It can also be found from Eqs. IV-8 and 12 by integration: 


■r 


dqz — kA 




(32 + 192 * - 307.2 dx 


Jx~o Jx=o Jx=0 

= hA{Z2 L + 96 - 102.4 1?) = (0.44)5(24 + 54 - 43.2) 

= 7^56 B/hr 


(c) From Eqs. IV-9 and 12 the following is obtained: 

— = a(32 + 192a; -307.2a:2) = 0.016(32 + 192* - 307.2*^) 
dr 

This means that the ordinates of curve d (Fig. IV-3d) must be multiplied 
by a = 0.016 in order to get dt/dr. Thus at 

X = 0 dH/dx^^ = 32 F/ft2 dt/dr = 0.51 F/hr 

X = 0.312 - 62.1 F/ft2 0.99 F/hr 

X = 0.75 = 3.2 F/ft2 0.05 F/hr 

The temperature of the wall does not change by the same amount at each 
point in the wall, but by about 0.5 F/hr at the heated surface, 1.0 F/hr at a 
depth of about 3 in. and 0.05 F/hr at the cooled surface. 
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Example IV-2. Answer the same questions as in Ex. IV-1 for a temper- 
ature distribution, given by the equation 


Solution: 


i = 90 - 8a; - SOx^ 

[IV-13] 

dt 

-r- == — 8 — *160a; 
ax 

[IV-13a] 

dH 

= -160 

[rv-136] 


(a) According to Eqs. IV-2 and 3: 


50 

5L 



-(0.44)5(~8) = 17.6 B/hr 


-(0.44)5(-128) - 281.6 B/hr 


(6) 53 = 5o ~ 5L = —264 B/hr 

The negative sign shows that no heat is stored, but more heat energy is given 
up at a; = L than received at a; = 0. 

This can also be found from Eqs. IV-8 and 136. 

53 = r dqz - -0.44(5)160 r ^ dx -0.44(5)160(0.75) 

Jx~0 Jx=0 

= -264 B/hr 

(c) From Eqs. IV-9 and 136 the following is obtained: 


•— = 0.016(-160) = -2.56 F/hr 
or 


Here the temperature of each part of the wall decreases equally. 

As in Figs. IV-Sa, 6, and d, the distribution of dt/ dx, and dH/dx^ is rep- 
resented in Figs. IV-4a and 6. 

It is easily understood that the difference in these examples is due to 
differences in the curvature of curve a in Figs. IV-3a and 4a. In the 
former the curve is concave, seen from above, and in the latter it is 
convex. Therefore, the temperature gradient at the left is greater than 
at the right in the first, and smaller in the second case. . This indicates 
that in Ex. IV-1 more heat is supplied than given up, in Ex. IV-2 the 
opposite is true, and by this dt/dr is positive in the former, negative in 
the latter. That the curve d simplifies to a horizontal line in Ex. IV-2 
comes from the simple form of Eq. IV-13. 

Finally it must be kept in mind that dt/dr as found in Ex. IV-1 holds 
only for the instant r. At a later instant the temperature distribution 
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across the plate will be changed. Assuming, however, that dt/dr is 
■practically constant over a reasonable time interval, for instance 1 hr 
in Ex. IV-1, it is possible to find a temperature curve like curve a in 
Fig. IV-3a for a later instant. After 1 hr, according to this calculation, 
t will have increased by 0.5 F at a: = 0; 1 F at rr — 0.312; 0.05 P at 
X = 0.05. Using the new cutve, which is higher, and repeating the 
construction given by Figs. IV-3a to d, the temperature distribution 
for any later instant may be found. 



Fig. IV-4a. Temperature distri- 
bution (Ex. IV--2). 



0 

-20 
-40 
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-80 
- 100 ^ 
.- 120 ^ 
-140 I 
-160 ' 




Fig. IV-46. Temperature gradient 
and its slope (Ex. IV-2). 


Because Ex. IV-2 led to dt/dr = —2.56 F/hr, independent of x, 
curve a in Fig. IV-4a is shifted down in one hour by —2.56 F as a 
whole, and the curves b and d in Fig. IV-46 obviously will remain un- 
changed. Thus, in this case the temperature of any point in the wall 
would decrease continuously by the same amount. This, however, as 
will easily be understood, is only possible if the amounts of heat supplied 
in time unit at the left and rejected at the right surface remain un- 
changed. 

Equation IV-9 has been derived for heat flow only in the direction x. 
For heat flow in any direction in a homogeneous body, and temperature 
changes in the directions of the three coordinates cc, y, and z, a similar 
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derivation leads to the more general equation: 

dt dH dH 


IV-2 Steady State as a Special Case of the General Equation 

Since Eq. IV-9 is entirely general for conductive heat flow in one 
direction if the thermal diffusivity a is invariable, the equation must 
hold also for steady state conditions. 

In this case the temperature at any point is constant with respect to 
time, i.e., dt/dr = 0. Therefore, total differentials can be used again 
and Eq. IV-9 becomes 


dH 

dx^ 


= 0 


[IV-15] 


Integrating once gives 

dt 


I C' being the constant of integration. A second integration gives 
I t = C'x + C" [IV-16j 

iC" being a second constant of integration. 

I The boundary conditions are as follows (see Fig. III-l): When 
li* - xi,t = h, and when x = X 2 , t = tz. 

1“ Substituting the boundary values in Eq. IV-16 gives ' 

I . ti = C'xi + C" 


^Giving the two equations for C' and C'' 

C' = and 

- X2 


.^1 - 


I Substituting these relations in Eq. IV-16 gives the following relation 
ior determining the temperature at any position x within the wall, which 
|s called the temperature field equation: 


(tiX2- -x(ti-t2) 


[IV-17] 


I The solution of Eq. IV^9 is, in general, not so easy as in the very 
limple case of steady flow conditions just considered. The integration is 
(beyond the level of this text on the elements of heat transfer, There- 
i|)re, in the following, only a few of the solutions will be given, but no 
lerivations. ' . / J 
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IV-3 Sudden Change of the Surface Temperature of a Thick 
WaU 

First an infinitely thick wall may be considered which is 
at temperature U throughout. The surface temperature cha_^^ 
suddenly to ts and remains constant thereafter. . -Pere the integratioj 
of Eq. IX-9 leads to 

(ti - 

In this equation t — the temperature in the wall at a distance x frol 
the surface and at the time r after the sudden temperature change ^ 

the surface, and f( — ) is a function which is called Gausses errg 

\2 V ar/ 



Fig. IV-5. Gauss's error integral. 


integral. Numerical values of this function are given in Table IV- 

and plotted in Fig. IV-5. It is seen that — = is a dimensionless ei 

2 V ar 

pression if x, r, and a are taken in consistent units, for instance x inf 
T in hr, and a in ft^/hr. 

The rate of heat fiov/ through an area A in the direction from 
surface -into the infinitely thick wall at a distance x and at the time 
may be found from the equation 

- t,) ^-=- [iv-l 


Values of e are plotted in Fig. IV-6. 
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TABLE IV-1 

Values of Gausses Error Integral f/ — ^ 

2V ar 


X 1 

J ^ \ 

X 


X 

J ^ \ 

2V^ 


2VZ 

\2^ or) 

2V^^ 

\2V^r) 

0.00 

0.0000 

0.4 

0.4284 

1.3 

0,9340 

0.01 

0.0113 

0.5 1 

0.5205 

1.4 

0.9523 

0.02 

0.0226 

0.6 

0.6039 

1.5 

0.9661 

0.04 

0.0451 

0.7>'. 

0.6778 

1.6 

0.9763 

0.06 

0.0676 

0.8 i 

0.7421 1 

1.8 

0.9891 

0.08 1 

0.0901 

0.9 

0.7969 

2.0 1 

0.9953 

0.10 

0.1125 

1.0 

0.8427 

2.2 

0.9981 

0.20 

0.2227 

1.1 

0.8802 

2.5 

0.9996 

0.30 

0.3286 

1.2 

0.9103 

3.0 

1.0000 


Taking a; = 0 in Eq. IV-19 one finds the heat flow into the surface 
of the wall at the time instant r: 


IcA. (ti fg ) 


[IV-20] 



0.4 0.8 1.2 1.6 2.0 2.4 2.8 3.2 


4a r 

Fig. IV-6. The e-f unction in Eq. IV-19- 

The heat energy which has entered the wall through the surface in 
the whole time r, beginning with the instant of sudden change of the 
surface temperature will be 

/»T=T 

Q ' dr 
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By substituting q from Eq. IV-20 and by taking all constant fact 
out of the integral, one obtains 




— P dr 
V T 


and by integrating 


Q = —kA(ti — 



[IV- 


Example IV-3. Calculate the temperature in a plane 7 in. from the s 
face of a very thick wall, and also the heat flowing into 10 sq ft of this pi 
10 hr after the surface temperature of the wall changes suddenly from 7 
to 1500 F and remains constant thereafter and also find the total heat ene 
taken up by the wall in 10 hr. Assume that the average thermal conducts 
and diffusivity values of the wall material are 0.5 B hr“^ ft“^ F“^ and ( 
ft^ hr“^ respectively - 


a; ^ 7/12 

2\/o.03(10) 


0.532 


With this value for — and according to Fig. IV-5 the value of the fu,^ 
2Var 

tion f( — ) is equal to 0.55. Substituting in Eq. IV--18 gives 

\2VaT/ 

t = 1500 + (70 - 1500)0.55 = ISOO^- 786 = 714 F 


at the required position and time. In order to determine the heat flow, coi 
putea;Y4 ar = 0.284 and here with, by referring to Fig. IV-6, find _ q - j 

Substituting the various known values in Eq. IV-19 gives 

0 75 

q =: -(0.5)10(70 - 1500) 4g : .. = 5520 B/hr 

V7r0.03(10) 


The total heat energy taken up by the wall is equal to that which ! 
entered the surface. For r = 10 hr, Eq. IV-21 leads to 


Q = -(0.5)10 (70 - 1500) 2. 




600 B- 


If the considered wall is of finite thickness L and — 7= > 0.6, * 

2 V ar 

same equations as those for the infinitely thick wall may be used. 
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i Sudden Change of the Surface Temperature of a Sphere or Cylinder 

or a sphere or a cylinder of diameter D, originally at temperature ii, 
ughout, and sudden change of the surface temperature to a penna- 



ACXT 

Fig. IV-7. Function F of Eq. IV-22. 


•t value Williamson and Adams (Ref. IV-1) presented the follow- 
equation: 

= h + (U - Q f (^) 


[IV-22^ 




COXDUCTIOX OF HEAT IN THE UNSTEADY STATE [IV -5 




In this equation tc is the temperature in the center of the sphere, or in 
the center line of the cylinder, r hours after the sudden change of the 

surface temperature. Any value of the function i^ay be ob- 

tained from Fig. IV-7 as ordinate at the abscissa 4 ar/D^. In this figure 
curves are given for spheres and cylinders of length L = D and L = x . 
For a cylinder of finite length L > D^m interpolation will give approxi- 
mate values. 


Example IV-4. Calculate the temperature at the center of a steel sphere 
16 in. in diameter 20 minutes after the surface temperature changes from 
80 F to 1200 F. An average value a = 0.31 for steel may be assumed. 
Solution: 

4ar _ 4(0.31) (20/6 0) ^ ^ 

D2 (16/12)2 

From Fig. IV-7 


0.20 


Substituting in Eq. IV-22, the final temperature at the center is found to be 
t = 1200 + (80 - 1200) 0.20 = 976 F 

IV-5 Periodical Change of Surface Temperature 

A thick plane wall where the surface temperature changes according 
to the sin-function, as shown in Fig. IV-8, will now be considered. 
Owing to the resistance of the wall material to the flow of heat, a definite 
period of time will be required before the surface temperature change 
affects the temperature distribution within the body. The time lag or 
the time required for the temperature at a given point within the body 
to be influenced by a surface temperature change may be found from 
the following relation: 


[IV-23] 

airn 

In this relation x represents the distance from the surface, n the number 
of complete changes per unit of time, and a the thermal diffusivity of 
the material. In Fig. IV-8 abscissas are values of time r, and ordinates 
are temperatures. A surface-temperature curve is drawn in full lines, 
and the temperature at a distance x from the surface is shown by the 
dotted curve. .The time lag Ar is represented as the distance between 
two corresponding positions off both curves. From Eq. IV~23 it is 
seen that Ar becomes larger the deeper the heat wave enters into the 
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wall, and that it decreases with increasing frequency n and thermal 
diffusivity a. 

The maximum deviation of the temperature from its average 
at the position x within the body, may be found by use of the following 
relation: 

= A^o.m • [IV-24] 

where Aifo.Tw represents the surface-temperature amplitude. From 
Eq. IV-24 and Fig. IV-8 it is seen that the temperature oscillates be- 
tween ta + and ta “ Atx,m where ta is its average value. The 



_^Time for One Complete Surface^ 

Temperature Variation 

► Time r 

Fig. IV--8. Temperature variation • inside a wall for sin-variation of surface 

temperature. 


amplitude of oscillation Atx,m decreases with increasing distance x from 
the surface. The greater the thermal diffusivity a of the material and 
the smaller the frequency n, the higher will be the value of 
The temperature at any point x within the body after t hours may 
be calculated from 



This actually is a combination of Eqs. IV-23 and 24. 

Surface temperature variations which approximate this condition are 
found in the field of heating, ventilating, and air conditioning. Because 
this particular type of temperature variation is important, several 
examples will be presented. 
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Example IY-5. A very thick brick waU is heated and cooled periodically 
every twenty-four hours according to a sin-variation of the temperature. If 
the surface temperature range is from 80 F to 150 F, calculate the time lag 
of the temperature wave at a point 8 in. from the surface. Assume that 
h - 0.6 B hr-ift-i F-^ Cp = 0.25 B lb;„--^F“S and 7 - llOlb/ft-^ 

k 0 6 

Solution: Substituting a; = 8/12 ft, a: = — = ftVbr and n = 


1/24 in Eq. IY~23 give 


At 


8/12 

2 


/Ti0(0.25) 

Vo, 


(110) (0.25) 


= 6.2 hr 


Example IY--6. If the range in temperature during twenty-four hours 
at the earth’s surface in a locality is —IS F to 25 F, calculate the range in 
temperature at depths (a) of 2 ft and (h) of 6 ft and the lag of the temperature 
wave at these depths if a sin-curve variation in surface temperature exists. 

Assume that ^ = 1.3 B hr“^ ft“^ F“^, Cp == 0.4 B lbm~^ F-^ and 7 = 
126.1 lb/ ft--3 


Solution: A^om = 


25 - (-15) 


= 20 F; n 


k 


1/24; a = f = 


1.3 


(126.1) (0.4) 


1 4 - 

0.0258 ftVbr; ta = - - 5 F; 




1 ^ 

1 


(a) 

—2^ 

= 20 e 

V 0.0258 

24 

= 0.22 F 


Atz,m, == 20 e 

/ ’T , 

T 


(6) 

V 0.0258 

24 

= 0.00 F 


Thus, at a depth of 2 ft the temperature range is from 4.78 F to 5.22 F and 
at 6 ft the temperature is constant and equal to 5.0 F. 


(a) 


(jb) 


At 


At 


2\0.i 

2^0 j 


24 


02587r 


= 17.2 hr 


24 


0,02587r 


- 51.6 hr 


Example IV-7. The twenty-four-hour range in temperature at the earth’s 
surface in a given locality is from — 10 F to 10 F. Determine (a) the ampli- 
tude of temperature oscillation at a depth of 1 ft, (6) the time lag of the tem- 
perature wave at a depth of 1 ft, and (c) the temperature at a depth of 1 ft, 
five hours after the surface temperature reaches the minimum' temperature. 
Assume ^ = 0.3 B hr”*^ ft“^ F"\ Cp = 0.47 B lb^“^ F~^ and 7 == 100 lb/ ft"^. 
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Solution: 


ta 


10 % 10 


OF 


Aio.7 


10 - (- 10 ) 


10 F 


(a) 

(h) 


n = 


2 ^ 

24 




10 6 
At 



1 

i 0.00638 

*24 


= ^ I- 

2\/0j 


24 


OOOSStt 


= 10e~^-^ 0.11 F 

= 17.3 hr 


(c) By substituting a; = 0 in Eq. IV-25 one obtains the temperature devia- 
tion from the average A^o = Ato^m sin 2'K'nr, and from this A^o = 0 for r = 0 
in conformance with Fig. IV~8. After that, the minimum temperature on 
the wall occurs when A^o = -Aio.m and this occurs the first time when 
3 3(24) 

27rnrm = f tt or at = — = 18 hr. Five hours later r = 18 + 5 

4n 4 

= 23. At this time and at x - 1 ft according to Eq. IV-25 


Atx = 10 e sin 



10(0.011) sin 1.49 


= 0.11 sin 85V = 0.11(0.996) = 0.11 F 
The temperature tx — ta At^ = 0 + O.ll — 0.11 F 


PROBLEMS 

/ 

IV-1. A reinforced concrete foundation slab 28 in. thick is to be poured into a 
foundation form, the bottom of which is earth at a temperature of 10 F. If the 
concrete is originally at a temperature of 65 F, calculate approximately the distance 
the freezing temperature will advance before the slab sets. Assume that the initial 
set takes place in sixty minutes, the final set in ten hours, and the heat generated in 
the concrete during setting is neglected. Assume that = 0.74 B hr”^ ft“^ F“^, 
Cp = 0.211 B lbm“^ F-^ and 7 = 141.7 lb/ frl 

IV-2. Calculate the temperature which the center of a steel sphere 0.6 ft in 
diameter will have attained at a time fi^ minutes after the sphere has been plunged 
into a heat-treating bath at 100 F. Assume that the surface of the sphere immedi- 
ately reaches the temperature of the bath. Assume the initial temperature of the 
sphere to be 1100 F. Assume that h = 21.5 B hr~^ ft"^ F~\ Cp = 0.132 B Ibnj.'”^ F"“\ 
and 7 = 490 lb/ ft"^. 
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IV~ 3 . Calculate the temperature at the center of a steel sphere 16 in. in diameter 
twenty minutes after the surface temperature changes from 80 F to 1200 F. Assume 
a = 0.31ft2hr“^ 

IV-4. Calculate the temperature at the center of a long cyhnder 18 in. in diameter 
after being placed in the air at 70 F for one hour. Assume that the initial tem- 
perature of the cylinder was 1000 F and that the surface temperature immediately 
reached the temperature of the air. Assume that k = 22 B hr“^ ft“^ F"^, Cp = 
0.14 B IhnT^ F"^ and 7 = 490 lb/ fri 

IV-5. The wall of a large furnace is made up of refractory material 21 in. thick 
having a thermal conductivity of 0.7 B hr“^ ft"^ F”^ and a diffusivity value of 0.022 
ft^ hr"’^ A thermocouple 3 in. from the inside surface indicated a temperature of 
1250 F seven hours after the firing-up period. If the gas temperature on the inner 
surface is 2000 F, calculate the theoretical temperature at the thermocouple position. 
Assume that the innermost layer of the refractory assumes the temperature of the 
hot gases immediately after the firing-up period begins. The initial temperature 
of the wall is assumed to be at 65 F. 

IV-G. The twenty-four-hour range in temperature at the earth’s surface in a 
given locality is from -10 F to 10 F. Determine (a) the range in temperature at 
a depth of 1 ft, ( 6 ) the time lag of the temperature wave at a depth of i ft, and (c) 
the temperature at a depth of 1 ft, five hours after the surface temperature reaches 
the minimum temperature, k = 1.34 B hr”^ ft*"^ F”\ Cp = 0.2 B lbm“^ F”\ and 
7 = 126 lb/ fr^ 

17-7. Calculate the temperature of the earth 2 ft below a large rectangular gas 
conduit after fourteen days from the time the hot gases start to flow through the 
duct. Assume that the gas temperature is 800 F and the initial temperature of the 
soil is 70 F. = 1.4 B hr“^ ft“^ Cp = 0.2 B Ibrn”^ F“\ and 7 = 126 lb/ frl 

IV- 8 . In a recently published article (Ref. IV-2) experimental results are re- 
ported on the time lag of the temperature wave through a common brick wall 12 % 
in. thick. The actual time required for the temperature wave to advance from a 
point % in. from the outer surface to a position % in. from the inner surface was 
nine hours. Using a diffusivity value of 0.0194 ft^ hr"\ check the experimental 
time lag of the temperature wave by use of the method suggested in this text. 

If the range in temperature (defined as twice the amplitude) is 8 F at a plane 
11 % in. from the outer face of the wall, calculate the temperature at a plane 11 % 
in. from the outer face after fourteen hours from the time the plane under considera- 
tion reaches the maximum temperature. Assume that the mean temperature of 
the plane equals 76 F. 
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CHAPTER V 


STEADY STATE HEAT CONDUCTION IN BODIES 
WITH HEAT SOURCES 

V-1 Temperature Distribution in a Plane Plate in Which Heat is 
Produced Homogeneously 

The transient storage of heat energy has been considered in Chapter 
IV. The subject of heat sources within a heat-conducting body will 
now be discussed. Cases of this type occur in exothermic chemical 
reactions such as the combustion of fuel in the fuel bed of a boiler 
furnace. Another example is an electrical coil wherein heat is produced 
and is dissipated by thermal conduction. 

Only the simplest type of this process will be dealt with in this text. 
Figure V-1 shows a perspective sketch of a vertical section taken 
through a cylindrical coil. The rate of heat produced homogeneously 
in the electric wires per unit volume of the coil may be calculated from 
Ohm’s law. This unit volume mcludes the volume occupied by the 
insulation and the wires. For purposes of simplification the height of 
the coil will be assumed to be such that only heat conducted in radial 
direction must be considered. If the influence of the curvature of the 
coil is negligible, then the coil shown in Fig. V~1 may be replaced by a 
plane plate with heat sources uniformly distributed inside the plate 
(Fig. V-2). Here heat will be conducted only in the directions +x and 
—z. The two surfaces will be assumed at a constant temperature 4 
and are represented hy x = L and x = —L. 

When no heat sources existed, the heat energy stored in unit time in a 

slice of the plate of thickness dx was found to be kA — ^ ^ a^ccording to 

(/IC 

Eq. IV-6. If the rate of heat q'” per unit volume* is produced in the 
coil an amount q'" {A • dx) will likewise be stored. However, according 
to Eq. IV-15 the total rate of heat storage under steady state conditions 
must be zero. The following equation results. 

kA^dx + q'"A • cfo: = 0 


.-,+,"' = 0 tv- 

* In this text g' is in B hr“^ ft"h q" is in B hr ^ ft ^ and g'" is in B hr“^ ft"^. 
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where k represents the heat conductivity of a fictitious substance having 
the samp composition of insulation and wire as an actual volume of the 
coil. The value of k may easily be measured or calculated. The total 
differential is used in Eq. V-1 because in steady state conditions the 



Fig. V-1. Section through a cylindrical coil. 



x=-L .:r=0 \x=^L 


Fig. V-2. Plane plate as simplified representa- 
tion of a cylindrical coil. 


temperature t does not depend on the time r, but only on x. If 0 = 
t — U is introduced in Eq. V-1, that is, counting the coil temperature 
from the surface temperature as zero point, the equation is converted 
into the form '' 

Since t, is supposed to be constant, <feldx^ = dHfdx^. Integration of 
Eq. V-2 gives 

= - \rx^ + Mx + JV 


[V-3] 
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The constants of integration M and iV may be found from two boundary 
conditions. The maximum temperature, 0^, in the coil will occur at the 
mid-plane. So at x - 0, = 6^. Substituting this in Eq. V-3 gives 

N [V-4] 

d must decrease symmetrically in the direction of x and —x from the 
maximum value. This condition may only be satisfied if in Eq. V-3, 

M = 0 [V-5] 

Substituting Eqs. V-4 and 5 in Eq. V-3 gives 

/// 

[V-6] 


Since t = U or 6 = Q 2 i;t x = L and at a; = — L 



which when substituted in Eq. V-6 gives 



[V-7] 


This means that a parabolic temperature distribution exists across the 
cross section of the coil with the apex at the median plane. 


V-2 Maximum Temperature in an Electric Coil 

According to Vidmar (Ref. V-1) there exists a simple and practical 
relation between the maximum and average temperature of a coil. 

The electrical resistance of copper and other pure metal increases 
linearly with the temperature. For the sum of all the lengths of wire 
inside a slice of unit thickness of the plate 

R' = Rj{l+a) [V-8] 

where 

R' = the resistance of the sum of the wire lengths at a temperature t 
R/ = the same at a temperature ts, and 

€ = the temperature coeffiSfeiR of electric resistance of the metal 

wire., (e == 0.0024 F“^ ^-copper.) 

H For a slice of thickness dxj then 

E' = E/(l + €0) ‘ do; 

The rate of heat energy produced in the coil by a current I in the coil 
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according to the laws of Ohm and Joule is 

• dx = I^R'dx = PRJil + ed)dx 


Integrating over the thickness of the coil 


g = g'"A2L = • dx = + 

^x=—L ^x=-L 


ed)dx [V-91 


In electrical engineering the average temperature da of a coil when a 
current I is flowing is found by measuring its electrical resistance Ra 
and also measuring its resistance Rs at a very small current such that 
the coil does not become heated sensibly above ts. This gives another 
relation for the heat produced electrically in the coil. 


q = PRa = PR,(1 + eB^) = I^Rj2L{l + jV-lO] 

In this equation R^ = Rj2L, since Rj was defined as the resistance 
for the thickness 1, and the coil is 2L thick. 

Equating q from Eqs. V-9 and 10 it follows that 

^ r 9dx 


Substituting 6 from Eq. V-7 gives 




or 


or 


6m 

^ ~ (ia 



[V-11] 

[V-12] 


Thus by measuring the surface temperature ts of the coil by means of 
an attached thermocouple and by determining the average*- temperature ^ 
ta or da ^ ta — ts from the increase of electrical resistance when a current 
is passing through the coil, the maximum temperature inside the coil 
at that current may easily be found. This temperature by which th%. 
durability of a coil is limited can scarcely be measured by any other way ^ 
in practice. Placing thermocouples inside the coil is rarely possible 
owing to the difficulty of electric insulation of the thermocouples against 
the main current. 

Example V-1. A vertical cylindrical electric coil 1% in. thick having ^ 
an inside diameter of 10 in. is heat insulated on the bottom and top so that 
no appreciable amount of heat will flow in a vertical direction. 
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When located in an oil bath the surface temperature of the coil was found 
to be 70.5 F at a current and voltage of 0.009 amp and 0.1 volt respectively* 
At 2.5 amp the resistance of the coil was 12.1 ohms. However, the surface 
temperature remained unchanged. Calculate the maximum temperature 
inside the coil when the current was 2.5 amp. Assume that the wires were 
made of copper. 

Solution: The ratio of the radii 10/13.5 of the coil is so great that the cur- 
vature may be neglected. A current of 0.009 amp is negligible as compared 


with 2.5 amp because by the latter according to Ohm’s law 


(--Y 

Vo. 009 / 


77,000 


times as much heat is produced than by the former one. Thus the starting 
temperature of the entire coil was ~ 70.5 F. According to Eq. V-10 

Ra = -Ks(l + edo) 


Rs — 


0.1 volt 


0.009 amp 
€ = 0.0024 F-i 


11.1 ohm 


Thus 

and 


Ra = 12.1 ohm 


12.1 = 11 . 1(1 + 

0.09 

0.0024 ” 


0.0024^a) 
37.4 F 


According to Eq. V~ll 

= 1.5^^ = 56.1 F 


Therefore, the maximum temperature of the coil was 

+ = 70.5 + 56.1 = 12^F 

REFERENCE 

V-l. M. ViDMAB, Suggestion of an Addition to the Test Codes on Temperature 
■Rise,”* Elektrotechn. u. MaschinenhaUj 26, 49, 65 (1918). 

* Title translated by the authors. 



CHAPTER VI 


INTRODUCTION TO THE DIMENSIONAL ANALYSIS 
OF CONVECTION 

VI-1 The Nature of Heat Convection 

Heat conduction is considered as a molecular exchange of energy, 
and therefore it is governed by molecular laws, some of which have been 
mentioned in Chapter 11. The motions of molecules involved in this 
process cannot be seen by our eyes nor measured by ordinary hydraulic 
or aerodynamic instruments. In contrast to this, heat convection is a 
transportation of heat energy by fluids in ordinary motion which is either 
directly visible or can be visualized by the usual instruments for the 
measurement of the flow^ of fluids. Heat transfer by convection is the 
exchange of heat energy between moving parts of the fluid or between 
these and surfaces of different temperature. 

Since in this process heat is conveyed mechanically from one particle 
to another, it is obvious that the transfer of energy depends upon the 
motion of the fluid and is governed by the laws of fluid dynamics, in 
addition to the laws of heat conduction and heat storage which must be 
considered at the same time. 

From this it will be understood that heat convection is a very complex 
process and that the simplicity of Eq. 1-3 is delusive. As a matter of 
fact, the film coefficient h, defiped by that equation, is a function of 
many variables, such as shape and dimensions of the surface, kind, 
direction, and velocity of the flow, temperature, density, viscosity, 
specific heat, and thermal conduptivity of the fluid. 

The differentia] equations winch describe convective heat transfer 
belong to the most difficult class in theoretical physics, and only for a 
very few simple cases has it been possible to solve them under simplifying 
assumptions. 

The empirical treatment proved to be likewise unsatisfactory. 
Because of the many variables mentioned above, it was almost impossible 
to find any more general relations, i.e., equations which could be applied 
to other arrangements and conditions than just those%sed by a particu- 
lar investigator. 

Thus, the knowledge of heat transfer by convection was extremely 
meager until about thirty years ago when the application of the so-called 

60 
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*»£ 2 •**“ •« ‘b. 

. Before explaining this principle and showing how it is used in problems 

•*“ oS 

and the 

oTfS 

dffierent law. Whenever the fluid particles Lw in Shsta^lt^fal! 
lei to the axis of the pipe without radial components, the flow is Sled 
VISCOUS, ammar, or streamlmed. The condition where radial com- 
onents or eddies exist together with the fluid motion parallel to the pipe 

"’ll f Reynolds investigated tL 

imblem of the » smilar flow of fluids. In this study he considered L 
rcumstances which would make the floWdiffemntfluids throu^ 
ircular tubes of different diameters with differenrvdocfe “ similar ” 
n a famous paper published in 1883 (Ref. VI-1) he showed*S^he 
low of fluids IS similar when a certain » dimensionless ” group of van! 

Jmber Ms'""'’ ™ ^®3^°lds 


(Re) = 


Vgl^P VaDy 


pg ^ 


m-i] 


In these ratios D represents the diameter, Va the average fluid velocity 

nd p, 7, and p the den^ty, specific weight, and dynamic viscosity of the 
md respectively. ^ ^ ^ one 

The dynamic viscosity p is defined by Newton’s equation 


dy 


[VI-2] 


here R, is the shearingforce or fluidfriction between two parallel layeis 
a fluid, which have equal areas A, are separated by the distance 
t and are moving, one parallel to the other, with the velocities of p 
Ld y + & respectively. According to this definition dv/dy is the veloc- 
r slope or gradient perpendicular to the direction of the flow of the tw'O 
Lid layers. If A = 1 and dv/dy = 1,F, = p. 

» The two fct letters of Reynolds’ name are In general used as a symbol for 
iynolds number. It is recommended to use this symbol with parentheses. 



62 IXTRODCCTIOX TO THE DLMEXSIOXAL AXALYSIS 

The physical dimensions of m are found by Eq. VI-2 in the same w 
as were those for k by Eqs. II-2 and 3. 

I[ibl - IW .IM 

or 

M = [lb sec] [VP 

In referring to Sect. I~3 it must be kept in mind that [lb] in the Briti 
system of units is a force. Instead of [lb] one can substitute the mi 
unit [slug] using the equation of definition 

l[lb] = l[slug] l[ft sec""^] 

and one obtains 

[ju] = [slug ' . [VP 

The unit [sec] can be replaced by [hr/3600] in these equations. 

The equation analogous to Eq. VI-4 in the physical system of units 

[n] = [gram cm""^ sec*^^] [VP 

This unit is called a poise in honor of the French physician Poiseii 
who was one of the first to state the law of pressure drop in streamli 
flow in a tube. The unit poise is often used in physics, particularly i 
hundredth part which is called a centipoise. Water at 68.4 F has a vi 
cosity value of exactly 1 centipoise. Values for water at other temper 
tures in lb ft"*^ hr units are given in Table VI-1 and values for air 
Table VI-2. 

TABLE VI~1 

Dynamic Viscosity of Water foe Various Temperatures 


Temperature 

Dynamic Viscosity’ 

F 

10-9lbft-2hr 

40 

8.99 

50 

7.59 

60 

6.49 

70 

5.66 

80 

5.00 

90 

4.43 

100 

3.96 

110 

3.58 

120 

3.25 

130 

2,97 

140 

2.73 

150 

2.51 


* Data taken from “ International Critical Tabl^/’ Vol. 5, 1929. 
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TABLE VI-2 


fNAMic Viscosity of Aie for Vaeioi 

IPERATURES AT ATMOSPHERIC PrESSUI 


Temperature 

Dynamic Viscosity 

F 

10“^ lb fr^ hr 

40 

0.101 

50 

0.102 

60 

0.104 

70 

0.105 . 

80 

0.106 

90 

0.108 

100 

0.109 . 

no 

0.110 

120 

0.112 

130 

0.113 

140 

0.114 

150 

0.115 

160 

0.117 

170 

0,118 

180 

0.119 

190 

0.121 

200 

0.122 


* Data taken from A. Eucken and M. Jakob, Der Ch&mie-Ingenieut, Vol. I, Part 1, Chapt. IL 
Akademische Verlagsgesellschaft, Leipzig, 1933, 


The ratio ju/p which occurs in Eq. VI-1 is called kinematic viscosity 
and is denoted by v. Because the unit of the density p is [slug ft“”^] = 
[lb ft“*^ sec^ ft”^] = [lb ft~^ sec^], one obtains from Eq. VI-3, as well as 
from Eq. VI--4, 

[,] = = [ftVsec] [VI-6] 

-Pj 


Since in this equation only kinematic magnitudes (length and time), 
but no dynamic magnitudes (force or mass) occur, the kinematic 
viscosity is a very convenient name for v. In contrast to this, /x is called 
“ d 3 mamic viscosity in this text (Ref. VI-2) instead of the name 
“ absolute viscosity ” which seems to be less specific and rather in con- 
trast to relative viscosity ” than to “ kinematic viscosity.” 

Substituting v from Eq. VI-6 in Eq. VI-1, the latter simplifies to 

{Re) = — [VI-71 


I and it is easily shown that the Reynolds number {Re) is really a dimen- 
sionless magnitude. 

The advantage of using Eq. VI-7 is that neither force nor mass units 
are involved. Many authors, however, express {Re) in terms of the 
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so-called mass velocity M = VaP, with which Eq. VI-1 converts into 

m-7a] 

Using this definition has the advantage that M is independent of pres- 
sure, temperature, and phase of the fluid. Since M has the dimensions 
[slug ff^^ sec""^], ju must be expressed in [slug ft“^ sec""^] units according 
to Eq. VI-4. 

Because the unit “ slug occurs in both the numerator and denomina- 
tor of Eq. VI-7a, both may be divided by 32.2. Then, since 1/32.2 of a 
slug equals 1 Ibm, the following equations of units may be assumed: 

[M] = [Ib^ ft-2 sec*-'] 

and 

[ 11 ] = [lb,n ft“' sec~'] [VI-761 

The latter unit is only 1/32.2 of the value of the unit according to 
Eq. VI-4 or 3. 

If Eq. VI-76 is found in books and articles without the subscript m 
or “ mass in connection with the unit lb, it must be kept in mind that 
the pound used therein means the mass of a pound, and the force unit 
in this system is the so-called poundal, equal to 1/32.2 of the gravi- 
tational force which acts on a pound. 

However, because ji does not appear alone in Eq. VI-1, but only in the 
term m/Pj ^ similar conversion is possible as was noted for the product 
Cpp in Sect, IV-1. By multiplying both sides of Eq. VI-76 by g it 
follows 

[gn] = [lb/ ft“' sec“'] [VI-7c] 

Thus, the latter expression is not the unit of p, as sometimes is believed, 
but rather the unit of gfi. 

On the other hand it is allowable in the numerical calculation of (Re) 
according to Eq. VI-1, to replace p[lbw^/cu ft] by the numerical equal 
value of 7 [lb//cu ft] and p[\hm ft””' sec”'] by the numerical equal value of 
{gix) in [lb/ ft”' sec”']. This procedure simplifies the calculation if /xis 
given in the units used in Eq. VI-76 as often is the case. In this book 
however, in general, the unit [lb ft”^ hr] or the unit [slug ft”' hr”'] 
will be used, both belonging to the British technical system with the 
hour instead of the second as unit of time. 

Example VI-1. Calculate the Reynolds number for water at 80 F flowing 
through a tube 2 in. in diameter, if the average water velocitj^ and specific 
weight is 10 ft/sec and 62.4 fb/ft® respectively. * 

— - ^ 


(Re) = 


MD 
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Solution: According to Table VI-1 the value for the viscosity is 5.00 
(10"®) lb ft"2 hr.^,..-By^nverting to mass units 

K /X = 5.00(10-9)36002 = 0.0648 slug fr^hr-i 


p — - — = 1.94 slug/ft^ (see Sect. 1-3) 
32.2 


Substituting in Eq. VI-6 gives 

0.0648 


1.94 


0.0334 ft2/hr 


Then from Eq. VI-7 


. M<3W)(A) . 

0.0334 ’ 


It may be noticed that seconds, wherever they occur, are converted into 
hours. 

Example VI-2. If the mass velocity for a fluid at 80 F flowing through 
a 1-in. diameter tube is 30,(100 slug hr"^ ft“2, calculate the Reynolds number. 
Solution: Substituting in Eq. VI-7a gives 


_ (30,000) (A) 
0.0648 


38,600 


Since the velocity distribution is the same if the Reynolds number is 
the same, it is obvious that the latter is decisive for all characteristics 
of the flow of'^a fluid which depend on the velocity distribution, as for 
the pressure drop in a tube. 

If, for instance, it is desired to find the pressure drop at twice the Rey- 
nolds number of a specified state, this can be accomplished by doubling 
either the velocity or the pipe diameter. Or if^xperimental data 
dealing with the pressure drop for a fluid having %;^ematic viscosity vi 
flowing through a tube of diameter Di and velocity' are known, then 
it^jgossiblej^detgmm drop for a fluid having a kinematic 

viscosity V 2 == 2vi under the same conditions by conducting experiments 
with the first fluid using a velocity V 2 — vi/2 or using a pipe having a 
diameter D 2 - Di/2 because, according to Eq. VI-1, the result leads to 
the same Reynolds number. 

Because the velocity distribution in the cross section depends on the 
magnitude of {Re), Reynolds further guessed that the transition from 
the laminar type of flow to the entirely different conditions of turbulent 
flow might take place at a definite Reynolds number. Indeed his experi- 
ments showed that turbulence starts at Reynolds numbers of about 2000 
to 2100. These experiments have been repeated frequently; the best 
known led to {Re) ~ 2300. This is called the critical Re^molds number. 
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VI-3 General Basis of Dimensional Analysis 

The usefulness of the Reynolds number for correlating experimental 
results for the flow of fluids through pipes cannot be overemphasized. 

A general method by which groups like the Reynolds number are 
attained in different branches of science is dimensional analysis. The 
basic principle upon which this is founded is dimensional homogeneity, 
that is, dissimilar quantities cannot be added together to form a valid 
physical relation. Mathematically, this means that in any physical 
equation both sides of which can be written as power functions or as 
algebraic sums of power functions, the sum of the exponents of the basic 
units must be the same on the left and right sides. For instance, 
according to the laws of freely falling bodies, 

V = V^l [VI-8] 

where v is the velocity when the body has fallen through the distance I 
and g is the acceleration due to gravity. As a numerical example 

30 A = ,/2(32.2)^14ft 

sec \ sec^ 

Considering the basic units alone 

[ft^ sec“^] = sec~^ ft^^^] 


Obviously, the above rule for the sum of the exponents yields 

1 = 2 +^ exponents of ft 

and 

— 1 = — 1 for the exponents of sec 

Now assume that Eq. VI-8 is not known, but that it is only postulated 
that V depends on g and 1. The latter fact may be expressed by 

z; = Cgn^ [YI-9] 

where a and h are unknown exponents, and C is an unknown constant. 

Representing the fundamental imits of length and time by L and T 
respectively, an equation for the dimensions of Eq. ¥1-9 is 

and the rule of the sums of the exponents yields 
for L 1 = a + h 
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Solving these two equations for a and h 


r. — X 

a ~ 2 


and by substitution in Eq. V-9 

cV^ 

Therefore, the law, given by Eq. VI-8, has been found by the method of 
dimensional analysis and only the constant C remains unknown and 
must be determined by experiments. According to Eq. VI-8, 0 = ^/2, 

VI-4 Application of Dimensional Analysis to Free Convection 

Exactly the same method can be applied for the much more compli- 
cated cases of heat convection. This was done first by W. Nusselt who 
in 1909 developed basic equations for forced convection, and in 1915 
equations for free convection. The method will be shown by the appli- 
cation of the principle to the example of free convection of a fluid in 
contact with a warm vertical wall. Before any valid relations may be 
established, it is essential to know what variables are involved. 

Correct selection of the variables comes about through experience. 
For this particular case, the variables involved are listed together with 
the dimensions in Table VI-3. 

The dimensions of the heat transfer coefficient h can be determined 
from Eq. 1-3 exactly in the same way as those of the heat conductivity k 


TABLE VI~3 

List of Variables and Their Dimensions for Natural Convection 





Exponents of units of 


Variable 

Symbol 

Heat 

Time 

Length 

Mass. 

Temp. 



H 

T 

L 

M 

0 

Heat transfer coefficient 

h 

1 

-1 

-2 

0 

-1 

Height of vertical wall 

Temperature difference between the 

1 

0 

0 

. 1 

0 

0 

surface and the fluid 

At 

0 

0 

0 

0 

1 

Density of the fluid 

p 

0 

0. 

-3 

1 

0 

Dynamic viscosity of the fluid 
Coefficient of thermal expansion of 

IX 

0 

j 

-1 

”1 

1 

0 

the fluid ' 

' /S 

0 

0 

0 

0 


Acceleration due to gravity 

g 

0 

-2 

1 

0 

0 

Specific heat of the fluid 

Cp 

■1 

0 

0 

-1 

-1 

Thermal conductivity of the fluid 

h . 

1 

-1 

-1 

0 

-I 
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have been derived from. Eq. II-l in Eqs. II-2 and 3. Thus in British 
technical units 

[h] = [B hr-i F-'] [VI-10] 

All other dimensions of Table VI-3 have been derived before or are 
well known. 

A general relation which includes all the variables of Table VI-3 is 
h - [VI-11] 

The corresponding equation of dimensions is 

CL®e^ (Mir ^ y y {h 


The equations fpr the sums of the exponents become 
H 1 =n + p 
y — 2m ~ p 

L — 2 = a — 3/— -i+m — p 

M 0 f + i - n 

By means of these equations, five of the exponents may be expressed 
in terms of the others, say j, m, n. By simple algebraic operations the 
following is obtained: 

a = — 1 + 3m 

6 = i 

/ = 2m 
i = — 2m + n 
p = 1 - n 

By substituting these five exponents in Eq. VI-11, rearranging and 
collecting the terms containing like exponents 

This equation can be further simplified if the relation between g and /3 is 
considered. Free convection is due to a thermal buoyancy effect. If 
in a fluid at temperature t\ and density pi a unit of volume is heated to 
temperature ^ 2 ^ and, as a result, has a smaller density P 2 , then cording 
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to Archimedes’ law, the buoyancy B = pig ~ p^g acts upon the men- 
tioned unit of volume. However, with respect to the definition* of the 
coefficient of thermal expansion pr-= P 2 [l + dfe -.h)] and by substi- 
tution B = p 2 ^g{h - h)- From this, since ^ and g enter the equatidfils^ 
owing only to the buoyancy, they will always appear in the form of a 
product (/3fir). According to Eq. VI-12 this is possible only if j = m. 
Therefore, it simplifies to 

For ideal gases it can be proved exactly that ^ = l/Ti, where Ti is 
the absolute temperature of the main bulk of the gas-f 

It may be mentioned that as early as 1881, L, Lorenz by integration 
of the differential equations of fluid flow and heat convection came to an 
equation for free convection of gases on a vertical wall which is a special 
case of Eq. VI-13 with m = n = 34. The constant C was determined 
by Lorenz for different gases; it was found to be an individual constant 
for each of them. 

It may further be mentioned that Eq, VI-13 holds also for cylinders 
(wires or tubes) located horizontally in a fluid at rest. For these I is 
the diameter of the cylinder. 

VI-5 Main Advantage of Dimensional Analysis 

Examination shows that each of the three groups of magnitudes in 
this equation is dimensionless, like Reynolds number, which as explained, 
has the advantage that one can change any involved magnitudes 
together with the other magnitudes in such a way that the numerical 
value of the dimensionless group remains the same. By this the result 
of any variation of each of these magnitudes can be obtained by s im ple 
calculations if the influence of the change of only one magnitude in a 
group is found from experimental investigation. The three constants 
C, m, and n must be determined by experiments. 

What work and time may be spared by setting up an equation of 
dimensionless groups like Eq. VI-13 can be seen by the following esti- 
mation: Assume that in order to find h for all thinkable variations, one 
would proceed by changing first taking, for instance, five different 
cylinder diameters, then conducting experiments with each of these and 
with fluids of five different thermal conductivities k, further varying p 
five times, and so on each of the seven independent variables Z, k, p, 

AZ, IX, Cp. This would mean 5^ = According to Eq.. VI-13 

* This definition is slightly different from N^elt^s original definition. 

t This does not hold for Nusselt's original definition of /3. 
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on the other hand it would be sufficient to make experiments with five 
different values of the first group at the right, each with five values of the 
second group, that is, 5^ = 25 experiments, in order to get about the 
same result as with more than 78,000 tests with the direct method. 

Actually, of course, more tests than the minimum of 25 would be 
conducted; oh the other hand, not all the variations of the mentioned 
seven independent variables would be found necessary or even possible. 
However, the saving of work and time would still be enormous. As a 
matter of fact, thousands of random experiments conducted before with- 
out the use of the principle of similarity have not led to any reliable 
knowledge of the behavior of heat transfer, whereas nowadays, relatively 
small numbers of tests yield very elucidating results. 

However, it must not be forgotten that simplifying assumptions are 
inevitable in this method. For instance, since the temperature changes 
across the fluid, the same holds for the viscosity and thermal conductiv- 
ity which depend on the temperature; but it is almost impossible to con- 
sider this influence exactly in the dimensional analysis. For this and 
other reasons the results of dimensional analysis may not be very exact 
and often only show a qualitative agreement with experience. 

Similar dimensionless ratio equations as Eq. VI-13 have been worked 
out for many cases of heat transfer and are extremely useful for solving 
problems. 


VI-6 The Most Important Dimensionless Groups 

Dimensionless groups occur over and over again . The most important 
are the following:* 


Reynolds number 
Grashof number 


m 

(Gr)^ 


vDp vDy 

= or or 

M pg ' 

• At 


vD 

V 


or 


MD 


Prandtl number 
Nusselt number 


(Pr) = 


(iVw) = 





Here a diameter D is used as a characteristic length. In the case of a 
vertical wall, for instance, the height must be taken instead of this in the 
Grashof number (see Eq. VI~13). 

It may be mentioned that the Prandtl number depends only on the 
physical properties of the fluid, and therefore itself is a physical property 

* As (Re), also thfe other dimensionless groups iti' .general, are. denoted by the two 
first letters of the name of an investigator who is to be honored by the denotation. 
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of the flowing substance, whereas the other groups depend also on the 
geometrical arrangement, the flow, the temperature, etc. 

With the symbols, given above, Eq. VI-13 may be written in the 
abbreviated form 

(Nu) = C{GrT{PrY [VI-14] 

The use of dimensional analysis for forced convection will be dealt 
with in Chapter VIII. 


PROBLEMS 


VI“1. Verify the following statement: Divide the surface JSlm coefBcient in 
B hr"*^ ft“^ units by 0.2048 to convert it to kcal hr"^ C“h 
VI-2. Verify by dimensional analysis the equation given below for the heat 
transferred from a rigid body submerged in a flowing fluid. It is to be assumed that 
the body is maintained at a definite temperature which is higher than the fluid 
temperature. The variables other than the rate of heat transfer ( 5 ) are a linear 
dimension of the body (Q, the velocity of the fluid, (zi), the temperature difference 
between the body and the fluid (Ai), the specific heat of the fluid based on a unit 
volume (Cp), and the thermal conductivity of the fluid [k). 

. k J 

VI-3. Verify the equations given below which have been used for correlating 
the data for heat transfer for fluids flowing in turbulent motion through the inside of 
pipes. 

(Nu) = C{Re)\Pr)^ 



VI- 4 . What is the significance of G in the Mowing equation of definition of the 
Reynolds number: 


{Re) 
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CHAPTER VII 

HEAT TRANSFER BY FREE CONVECTION 

VH-l Some Peculiarities of Free Convection 

As already explained, the motion of the fluid particles in natural or 
free convection heat transfer is due to the differences in densities of the 
parts of the fluid. A domestic gravity hot-air heating system is an 
exampld of heat transfer by natural convection. Heated air rises 
through the furnace pipes, passes through the registers, is cooled at the 
walls, windows, etc., of the room, and then returns through the cold-air 
ducts to the furnace. Another example of natural convection is the 
heat transferred from the surface of an insulated pipe to the still air of 
a room. In this case, some of the heat is transferred by radiation. 
On the other hand, the transfer of heat from hot flue gases to the tubes 
of an economizer or water heater serves as an example of forced con- 
vection, since the flue gases are made to pass through the unit by fans. 

In Chapter VI it has been shown that natural convection heat trans- 
fer theory is complicated because of the many variables involved in the 
solution of a particular problem. It has also been mentioned that a 
complete mathematical solution is known for only a very few cases. 

The experimental data on natural convection for some cases have 
been correlated by dimensional analysis and these correlations are 
useful in the solution of problems. Approximate simplified equations 
which apply to a limited range of experimental data have been developed 
for several cases. Some of these equations will be given in this chapter. 
They should be used with considerable reservation since they are valid 
only for a limited range. 

There is another item which requires consideration. At in the basic 
Eq. 1-3 is called the temperature difference between the area A of the 
surface and the fluid in contact with it. The temperature of the surface 
is well defined and can be measured by attached or inserted thermo- 
couples. However, there is no unique temperature of the fluid. Ac- 
cording to experience, the temperature of the fluid’is equal to the surface 
temperature at the places where the fluid is in contact with the surface; 
then, if the surface is warmer than the environment, the temperature 
decreases sharply in a thin fluid layer close to the surface and finally it 
becomes almost constant at some distance from the wall. If the surface 
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is cooler than the environment, the behavior is similar, with a sharp 
temperature increase of the fluid beginning at the surface. In natural 
convection, At generally is defined as the difference between the temper- 
ature of the surface and that of the fluid at some distance, say several 
feet, from the surface. 

VII-2 Empirical Equations for Heat Transfer on Horizontal and Ver- 
tical Surfaces 

The following empirical equations (Ref. VII~1) can be used in the 
approximate solution of problems of heat transfer between horizontal 
plates and still air under free convection conditions. 

h = 0.38 [VII-1] 

for plates facing upward and 

h = [VII-2] 

for plates facing downward. 

It may be mentioned that At here stands for At degree Fahrenheit/ 
1 degree Fahrenheit because only dimensionless numbers can bear an 
exponent. This then means that for Ai = 1, = 0.38 B hr""^ ff"^ 

according to Eq. VII-1. . 

Similar instances of the use of exponents in this text must be under- 
stood in an analogous manner. 

Example VII-1. Calculate the heat loss by natural convection from a 
horizontal heated plate at 300 F facing upwar^d to the still air of a room at 
80 F. The plate is 7 in. by 9 in. 

Solution: 

h = 0.38(300 - 80)0-25 ^ 1 47 g p-i 

The total loss is 

1.47 (300 - 80) = 142 B/hr 

144 

For calculating approximately the heat transfer by free convection 
between vertical plates and still air the following formulas may be used 
(Ref. VII-1): 

U = tVII-3] 

for plates of height L > 1 ft, and 

h = 0.28 {-] {VII-4] 


for plates of height L < 1 ft. 
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Example VII-2. Calculate the film coefficient for free convection from 
a vertical heated plate 18 in. high at 200 F to still air at 60 F, 

Solution: 

h = 0.3(200 - 60)0*25 = 1.03 B hr-^fr^F-i 

For free convection between horizontal pipes and still air the following 
approximate formula is recommended (Ref. VII-1): 

/^A 0.25 

/i = 0.22(^-j [VII-5] 

where D is the external diameter of the pipe in feet units. 

Example VII-3. Calculate the heat loss by natural convection from a 
4^in. nominal horizontal steam pipe at 200 F to the still air of a room at 70 F. 
Solution: 

h = 0.22 poo 7 , 95 B iir-i p-i 

\ 4.5/12 

The heat loss per square foot of surface is 

qf' At ^ 0.95(200 - 70) - 123.6 B hr^i ff^ 

VII-3 Use of a Dimensionless Correlation 
The results of many investigators on the determination of the natural 
convection coefficients for horizontal cylinders in still gas have been 
correlated by dimensional analysis. This has shown that for gases of 


iGr) (Pr) 

10 “"^ 10”2 100 10 "^ 10 ® 10 ® 



Log iCr) iPr) 

Fig. VII— 1. Correlated data for determining coefficients of free convection between 
horizontal cylinders and diatomic gases. 
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equal numbers of atoms one can assume m = n in Eqs. VI-13 and 14. 
Thus the latter simplifies to 

(Nu) ^ C[(Gr) (Pr)r [Vll-fi] 

The data for diatomic gases are represented by the line drawn in Fig. 
VII-1 where both coordinates are dimensionless. This chart is very 
significant in that it represents the correlation of test data for fine 
horizontal wires in gases as well as for large pipes up to 10 in. in diam- 
eter, and temperature differences from cylindrical surfaces to the gas up 
to about 2000 F for wires and 700 F for the large pipes. For this par- 
ticular plot, all the fluid properties occurring in Eq. VI-13 are taken at 
an average temperature between the fluid and the surface and the value 
of the coefficient of cubical expansion jS is assumed to be equal to the 
reciprocal of the average absolute temperature of the gas, as is true for 
ideal gases (see Ex. VII-4, footnote). 

It is seen from the graph that when (Gr)(Pr) increases from 10"“^ 
to 10'^^, that is, by 12 decimals, Nu increases only from about 0.5 to 15, 
that is about 30-fold. 


Example VII-4. Calculate the heat loss by natural convection per square 
foot of outer surface from a horizontal 4-in. nominal steam pipe to the still 
air of the room if the air and pipe surface temperatures are = 70 F and 
= 200 F respectively. 

Solution: The outer diameter of the pipe D — 0.375 ft. The temperature 
at which the physical properties of the fluid (with the exception of are 
determined is equal to 


Un 


200 -f 70 


135 F 


At this temperature the involved physical properties of the air are: 

= 0.00208 slug ft-5 

32.2 

H = 0.1135(10-®)360(V= 0.00147 slug ffi hr-i (from Table VI-2) 
Cp = 0.238(32.2) = 7.64 B slug-i F"! 
k - 0.0161 B hr“^ ft“^ ,F”^ (from Fig. II-l) 

Further;* 




1 

Pi 


1 

460 -h 70 


0.00189 F"i 


* In Sect. VI-4 it is mentioned that for ideal gases /9 = 1/Pi where Ti is the 
absolute temperature of the main bulk of the gas. 
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Wherever pounds occur, they are converted to the mass unit slug, by dividing 
by 32.2. This is necessary because in the derivation of the Grashof number 
mass units are used (Eq. VI-11, etc.). The usual value of the specific heat 
is related to 1 lb mass; in order to relate it to 1 slug it must be multiplied by 
32.2. (See Sect. IV-1, the paragraphs before Ex. IV-1.) 

In addition to the properties of the air, the gravity constant g occurs in 
{Gr), Because the unit of time used above is the hour, it is necessary to take 
g = 32.2(3600)^ ft hr"^. Otherwise the consistency of the units in the di- 
mensionless group would be disturbed. 

Now, substituting all these magnitudes, one obtains 


m - 




= 10.75(10^) 


(Pr) = = 0.697, and 

h 

{Gr) iPf) - 7.49(10®) 

log {Gr) (Pr) = log 7.5(10®) = 6.875 


Using this number and referring to Fig. VII-1 one finds log {Nu) == 1,42 and 
from this {Nu) = 26.3. But {Nu) -\D/h Thus 


26.3(0.0161) 

0.375 


1.13 B hr-i ft“2 F“i 


and the heat transfer per square foot is equal to 

- 1.13(200 - 70) = 147 B hr"! fr^ 

Referring to the empirical Eqs. VI-1 to 5, it may be mentioned that these' 
too are special cases of Eq. VI-13. In all of them the result, already found 
by L. Lorenz and mentioned in Sect. VI-4, m = n = K, is used. 

VII~4 The Concept of a Gas Film 

The disadvantage of dimensional analysis is that it does not help to 
explain why such dimensionless groups as the Grashof number occur. 
A more conspicuous theory is based on the assumption of the existence 
of a gas film. Langmuir had already visualized the idea of a gas film in 
his doctor’s dissertation in 1906 (Ref. VII-2) and in 1912 he presented 
this concept in more detail (Ref. VII-3) as follows: Close to a heated 
wire practically no free convection occurs, and, as a result, a thin layer 
or film exists through which heat passes by conduction and radiation 
alone. It is only necessary to consider convection outside of the film. 

Langmuir came tg.|his concept through observations which indicated 
that close to the gloMng threads of electric lights convection was very 
smalL^ This he explained on the basis of the well-known increase of 
viscosity of gases with the temperature. 
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Rice (Ref. VII-4) combined this theory with dimensional consider- 
ations and arrived at formulas which are rather similar to those found 
by dimensional analysis alone. 

PROBLEMS 

VII-W Calculate the film coefficient of heat transfer by free convection and the 
total heat transferred from a heated vertical plate 1.5 ft by 2 ft at 200 F to still air at 
50 F. 

VI^2. Determine the total heat transferred by natural convection from a hori- 
zontal plate facing upward at a temperature of 300 F to air at 100 F. 

VII-3. Calculate the heat transferred per square foot by natural convection from 
a 6-in. pipe at 400 F in contact with air at 70 F by use of Fig. VII-1. 

VII-4. Calculate the heat transfer by natural convection from a horizontal wire 
0.05 in. in diameter and 16 in. long at 800 F to nitrogen gas at 100 F. 
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CHAPTER VIII 


HEAT TRANSFER BY FORCED CONVECTION 

Vm-l Some Peculiarities of Forced Convection 

. Forced convection is probably the most important method of heat 
transfer employed in engineering. It is used in almost every type of 
heat exchanger, at least for one of the fluids and often for both. Forced 
convection heat transmission also occurs in devices which are not 
classified as heat exchangers, such as furnaces with artificial draft 
and thermal engine coolers. It is only at high temperatures that heat 
transfer by radiation may be more effective. 

As in free convection (see Sect. VII-1) the significance of the tempera- 
ture difference At in Eq. 1-3 must again be considered. In forced con- 
vection the temperature difference is the difference between the well- 
defined surface temperature and the temperature of the fluid. The. 
latter temperature varies throughout the fluid. In a gas being cooled 
as it flows through a pipe, the fluid temperature is at a maximum value 
at the surface and at a minimum value at the axis of the pipe. No 
unique agreement exists as to whether the temperature at the axis, the 
linear average temperature along a diameter, or another mean tempera- 
ture should be used. However, the mixing-cup temperature tm is con- 
sidered probably the most convenient to use. This could be formed 
by perfectly mixing the fluid passing the cross-sectional area of the tube 
in a “ cup ” and measuring its average temperature. For turbulent 
flow the temperature of the fluid increases or decreases sharply in a very 
thin layer next to the surface, and is practically constant throughout the 
bulk of the fluid; for this the fluid temperature may be taken at any 
point in the cross-sectional area except close to the surface. The method 
for determining the average temperature with respect to the tempera-- 
ture change in the direction of the flow of a heated or cooled fluid will 
be discussed in Sect. IX-6. 

For free convection it has been shown by means of dimensional analy- 
sis that the Grashof number (Gr) and Prandtl number (Pr) are two 
groups which may be used as independent variables in setting up 
empirical equations and constructing graphs, and the heat transfer 
coefficient may be represented in terms of the dimensionless Nusselt 
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fluids in turbulent flow through pipes 


I amber (Nu), For forced convection {Nu) may be expressed as a 
Bction of Reynolds number (Re) and again the Prandtl group (Pr), 
In general, all the dimensional equations which represent the indi- 
■idual coefi&cient h for forced convection may be expressed by the form 

(Nu) - C{Re)^{Pr)^ [VIII-1] 

! where C, m, and n are constants. 

; Vin-2 The Heating of Fluids in Turbulent Flow through Pipes 

A good general correlation (Ref. VII-1 and VIII~1) of the data on 
heating of any fluids in turbulent motion (Re > 2300) through pipes is 

.. (f) 


m 


[vin-2] 


For cooling of fluids it can be likewise used approximately, lo 
The physical quantities are all determined at the average temperature 
of the fluid (Ref. VIII-1). 

Example VIII-l. Calculate the average film coefficient of heat transfer 
j at the water side of a single-pass steam condenser. The tubes are 0.902 in. 
inside diameter, and the cooling water enters at 61.4 F and leaves at 69.9 F. 
The average water velocity is 7 ft/sec. 

i Solution: h will be found by means of Eq. VIII-2. D = 0.902/12 — 0.0752 
!ft; V - 7(3600) = 25,200 ft/hr. The pertinent physical properties of water 
' at a mean temperature of 65.7 F, as taken from Tables II-2, VI~1, and from 
the Keenan and Keyes’ Steam Tables (Ref. VIII-2), are: 


P = ^ = 1-94 slug ft-’ 

, e.oiiOT’) lb ftr' ht = 6.o(io-’) (3600)’ slug, ft-' hrr' ^ 
= 0.0778 slug ft-' hr-' ' "" 

Cp = 1 B lb„-' F-' = 1(32.2) B slug-' F"' 
k = 0.341 B hr-' ft-' F-' 


^ 25,20m0_752)lj_4 ^ 

IX 0.0778' ’ 

(Pr) .afE 0-0778(32.2) ^3^ 

k 0.341 

(Nu) -■ — = 0.023(47, lOO)®-® (7.34)“-^ 279.5, and 


279.5(0.341) 


1270 B hr-' fr’ F"' 



80 


HEAT TRANSFER BY FORCED CONVECTION 


In the case of gases the Prandtl number chanps very little. It is 
about 0.67 for monatomic gases, 0.70 for diatomic gases, 0.89 for tij. 
atomic gases. 

The corresponding values of (Pr)°-^ are 0.852, 0.867, and 0.955. 
Taking approximately (Prf'^ = const. = 0.9 which corresponds t« 
_ 0.794 the error in (Nu) is in the range of only ±5 per cent for 
these three kinds of gases and Eq. VIII-2 simplifies to 


W 

k 


/'vDp\°-^ 

0.0207 


[VIII-3] 


Example VIII-2. If air enters a pipe of 4-in. inside diameter at 65 F and 
^f^eaves it at 300 F, what is the coefficient of heat transfer? Assume an avet- 
!> age air velocity of 15 ft/sec and atmospheric pressure. The^iscosity and 
' +.>»orma,1 conductivity values for a mean temperature of 183 Fjas taken from 
Table VI-2 and Fig. II-l are 0.119(10-8) ib ft-2 hr and 0.0171 B hr-i ff^ 
respectively. The specific weight is 0.0617 lb ft 

Solution: With D = Vi 2 it,v = 15(3600) ftAr, P = 0.0617/32.2 slug/{t», 
and ja = 0.119(10-8)3600^ slug.ft'’- hr-'^ one obtains 


A. 


m 


= !^ = (15)3600(^) 


10 ® 


0.0617 


32.2 0.119(3600^) 


= 22,400 


By substitution in Eq. VIII-3 


Nu = 


0.0171 


and h = 3,23 B hr"! ft“2 


VIIl-3 The Ffeatitig of Liquids in Streamline Flow Through Pipes 

The equations which are used for streamline flow in pipes wifii 
viscous liquids like oils are more complicated than those for turbulent 
flow. 

Sieder and Tate (Ref. VIII-1 and 3) in experiments with viscera 
petroleum oils in tubes with inner diameters ranging from 0.39 to 1.57 k. 
and heated lengths ranging from 3 to 11.6 ft, arrived at the following 
formula: 



[Vin^l 


(Nu) = 
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where h is related to the difference between the mean surface tempera- 
ture ts and the arithmetic mean of the “ cup temperature 
of the liquid at entrance and exit; 

L is the heated length of the tube; 

is the d3mamic viscosity of the liquid at temperature and 
jUs is the same at mean surface temperature ts. 

In the experiments ix^/iis ranged from 0.004 to 9.8. 

All factors in Eq. VIII^ are dimensionless. Comparison with 
Eq. VIII-1 however shows that there are two more dimensionless 
terms D/L and involved. Geometrical similarity of tubes 

obviously exists only if D/L is a constant. This is the treason why 
D/L appears in the equation as a parameter. In turbulent flow, the 
exponent of D/L is about Ks and (D/L)^^^^ approaches a practically 
constant value if L is relatively great compared udth D. Thus, in 
turbulent flow in relatively long tubes Eq. Vin-2 can be used without 
taking care of the small variability with D/L in that range. This is 
not true with streamline flow where, according to experience, the ex- 
ponent is only The dimensionless ratio fXs considers the influence 
of temperature on viscosity which already has been mentioned in Sect. 
VI-5. By means of this term Eq. VIII-4 can be used for heating and 
for cooling of viscous oils in tubes, whereas without this term different 
equations for these two cases would be necessary. 

VEH-i The Heating of Fluids Flowing Normal to Single Wires and 
Tubes 

For this case, as for gas. flow inside of a tube, one can use the equation 

(Nu) = C(Rer [VIII-5] 

Different investigators found rather different values for C and m. As 
an average C = 0.3 and m = 0.57 may be assumed. 

Example VIII-3. Calculate the heat transfer coefficient for the case of 
air flowing normal to a 1-in. pipe at a Reynolds number of 8000 if the average 
surface temperature of the pipe is 183 F. 

Solution: D = ft, and k = 0.0171 B hr"'^ ft”^ F^^ (see Ex. VIII-2). 
By substitution of {Re) = 8000 in Eq. VIII-5 

(Nu) = 0.3(8000)0*57 = 50.0 

and 

* ^ = 50.0 _ (0. 0 1^ ^ 3 g 
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Since the constants of Eq. VIII-5 are not known exactly, it does not 
introduce an appreciable error if fc is taken at the surface temperature 
instead of at an average between the gas and surface temperatures. 

For the heat transfer from tubes to liquids, Eq. VIII-1 may be used 
with C = 0.6, m = 0.5, and n = 0.31 in the range from = 50 to 
Re = 10,000 (Ref. VIII-4). 

VIII-5 The Heating of Fluids Flowing Normal to Banks of Staggered 
Tubes 


Here Colburn (Ref. VIII-5) recommends the following relation for 
correlating the variables which influence the heat transfer coefiScient; 

,- 0.4 

[VIII-6] 


XVmPCpJ \ k / \ M / 


In this formula D is the outer diameter of the tubes, and v^, is the 
mean velocity in the smallest free cross-sectional area between the 
tubes passed by the fluid. The viscosity and thermal conductivity 
values are determined at the average film temperature, 

tf = tm + ** “ [VIII-7] 


where and ts are the mean fluid and surface temperatures respectively. 
Equation Vin~6 can be converted to 

or [VIII-8] 

(Nu) = 0M(Remf-^(Pry^^ 

as can be easily checked (Ref. VIII-l). 

Vin~6 The Surface Film Theory of Heat Transfer 

The equations appearing in this chapter are based on dimensional 
analysis. However, as for free convection, there exists a surface film 
theory based on more prominent facts. 

In 1904 L. Prandtl first directed attention to the different behavior 
of the bulk of a fluid in turbulent motion and of a thin layer on the 
solid walls which limits the flow. In the former mixing movements 
prevail, whereas in the latter gas friction is characteristic for the motion. 
In 1910 Prandtl showed in an analogous manner that heat is transferred 
in the bulk of the flow by mixing movements only, and in a surface film 
by thermal conduction alone. Using Reynolds' idea that gas friction 
and heat transfer are governed by the same laws (see Chapter XY), 
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Prandtl developed a theory concerning this analogy which is known as 
the surface film theory of forced heat convection. This theory is still 
in a state of development. 

As a result of Langmuir’s concept of a gas film, which incidentally is 
very similar to Prandtl’s assumption, Rice developed formulas for forced 
convection (Ref. VIII-6) which, like those for free convection, are not 
very different from the relations found by dimensional analysis. 

PROBLEMS 

VIII-1- A cylinder 1 in. in diameter is placed in a hot air duct at right angles to 
the direction of the air fiow. If the air and cylinder temperatures are 100 E and 
300 F respectively, calculate the value of the film coefficient for a mass flow velocity 
of 8000 Ihm hr"^ ft~^ (see Sect. VI-2). 

VIII-2. Water at a rate of 5000 lb /hr flows through a 2-in. internal-diameter 
tube. The temperatures of the water entering and leaving are 80 F and 120 F 
respectively. Calculate the heat transfer coefficient at the water side. 

VIII-3. Air flows through a tube having an internal diameter of 4 in, at a volume 
rate of 10,000 ftVhr. If the entrance and discharge temperatures are 60 F and 
200 F respectively, calculate the heat transfer coefficient. Assume an average air 
pressure of 40 Ib/sq in. 

VIII-4. Determine the coefficient of heat transfer for oil flowing through a 
heated tube 8 ft long and 1 in. in diameter at an average velocity of 1.5 ft/sec. The 
entering and discharge oil temperatures are 120 F and 140 F respectively. Assume 
that the tube surface temperature is constant and equal to 160 F. The specific 
gravity,* average thermal conductivity, and average specific heat values for the oil 
are 0.88, 0.08 B hr-'^ ft“' F“h and 0.5 B lbm“^ F-^ respectively. The variation of 
the dynamic viscosity of the oil for various temperatures may be obtained from the 
following data: 


Temperature 

Dynamic Viscosity 

F 

10-2 j]3 ft-2 seQ 

120 

1.3 

130 

1.1 

140 

0.9. 

150 

0.7 

160 

0.6 


VIII-5. Calculate the heat transfer coefficient for water flowing over a bank of 
tubes 1 in. in diameter at a mean velocity through the smallest free area of 8000 
ft/hr. Assume that the average film temperature is 1 18 F. 
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CHAPTER IX 


HEAT TRANSFER BY THE COMBINED EFFECT OF 
CONDUCTION AND CONVECTION 

K-l Cases of Combined Conduction and Convection 

Thus far only such cases have been treated where either conduction 
and convection occurs alone or the effect of the one prevails so much 
that the other can be neglected. This is not generally true. 

There are two main cases in which the combined effect of conduction 
and convection must be considered. The one deals with the projections 
of a surface into the environment, as for instance, the handle of a hot 
vessel or fins on the cylinder of a combustion engine. Here heat is 
conducted from the root to the free end or edge of the projection, and at 
the same time heat is given up to the environment by convection and 
by radiation. The radiation effect will be neglected or simply con- 
sidered as included in the convection effect in this chapter. 

The second case occurs at the wall of a building and in heat exchangers. 
In a building, heat energy may be transferred from the inside air to the, 
wall, through the wall by conduction, and finally given up to the outside 
air. Thus, there is a series of one convective and one conductive and 
another convective transfer of the same rate of heat, similar to the heat 
flow by conduction through a series of different layers of insulation of a 
composite wall which has been dealt with in Sect. III-3 and 6. Con- 
sider the transfer of heat from the hot water to the cold water in a 
tubular water cooler. The temperature of the hot water is reduced as it 
passes through the apparatus owing to the gradual dissipation of thermal 
energy of the hot water particles to the cooler surface of the walls which 
separate the hot and cold water. The thermal energy then flows 
through the metal wall by conduction and finally into the cold water by 
convection. 

IX-2 Heat Transfer from a Rod Heated on One End 

The first of these cases will be treated by using as an example the heat 
transfer from a rod of arbitrary constant cross section protruding from 
a vessel at constant temperature. 

In Fig. IX-1 a rod whose cross section has an area A and a circum- 
ference C protrudes by the length x = L from a perfectly heat-insulated 
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liquid bath at constant temperature <o into a room in which the air 
temperature is ta- The rod is made of a material having a thermal con- 
ductivity of k, and the film coefficient of heat transfer from the rod 
surface to the air may be h, both k and h being assumed as constant in 
the range of temperatures considered. It is further assumed that 
temperature differences in any cross-sectional area of the rod are negli- 
.gibly small. 



Fig. IX-1. Rod protruding 
from a liquid bath. 



Fig. IX-2. Temperature 
distribution along the rod, 
shown in Fig. IX-1. 


First consider the heat balance for a differential element contained 
within two planes perpendicular to the rod axis at the distances x and 
X + dx from the body at constant temperature to. There is a heat flow 
by conduction qa, originating from the liquid bath, and another one qh 
in the direction of the free end of the rod. Furthermore, an amount qc 
is given up by convection from the rod to the air. Under steady state 
conditions, there will be a heab balance of incoming and outgoing heat 
energy 

qa = qb + qc [IX-1] 

Denoting the temperature at the distance x from the vessel by t, then 
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according to Eq. III-l 


1 

1! 

[IX-2] 

Furthermore, according to Eq. IV-5 


M ^ dx 

[IX-3] 

Finally, according to Eq. 1-3 


qc = h{C • dx) {t - ta) 

[IX-4] 


becaxise the surface of the differential slice through which heat is given 
up to the environment is C ■ dx, and the difference between its tempera- 
ture and that of the air is t — ta. 

Substituting qt, and q^. from Eqs. IX-3 and 4 in Eq. IX-1 one obtains 


dH 

Ua = qa — kA ^dx + hC • dx(t — to) 


or 


dH 

dx^ 




[IX-5] 


Let d - (t — ta) be the temperature of the rod, referred to the air‘ 
temper ature a s zero, and simplify by introducing the symbol 
m = 's/hC/hA, 

Then, because ta = const., = d^djdx^. Thus Eq. IX-5 

simplifies to 

J2a 

^ = [IX-6] 


The general solution of this differential equation is 
6 = Me-^^ + 


[IX-7] 


where M and N are two arbitrary constants of integration. It is easy 
to check Eq. IX-T. Differentiating once gives 


dx 


= -Mme-^ + Nme^^ 


[IX-8] 


Differentiating a second time gives: 
d^d 


dx^ 


+ Nmh^ = + Ne^^) 
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and by substitution from Eq. IX-7 


d^9 

dx^ 


= m^B 


[IX-6] 


The constants M and N will be found from the so-called boundary 
conditions, two of them being necessary and sufiSeient. 

If, for instance, the free end of the rod in Fig. IX-1 were perfectly 
insulated by a pad P, then qa would be zero at x = L, that is, no heat 
could leave the rod in the axial direction at x = L. Thus according to 
Eq. IX-2 

or 



Substituting this and x = L in Eq. IX-8 gives a first boundary condition. 
0 = -f NmeT^^ 


or 


^2mL 


M 

N 


[IX-9] 


A second boundary condition is found from the fact that owing to the 
perfect insulation of the liquid bath, at x = 0 the rod temperature is 
t — to, or 8 = Bq. Substituting this in Eq. IX-7 gives 

00 = Me° + Ne° = M +N [IX-10] 

From Eqs. IX-9 and 10 it follows 

^ = and A = tK-ll] 

and by substitution in Eq. IX~7 

® = ^0 (l + + 1 + 

Example IX-1. A cylindrical rod of %-m. diameter and 9 in. long pro- 
trudes from a well-insulated steam vessel at 210 F into air of 80 F. The free 
end of the rod is insulated so that the heat losses from this end may be negli- 
gible. Determine the temperature at the free end if the rod is made of (a) 
silver, (6) iron, and (c) glass. Show the temperature distribution along the 
rod in these three cases. Take the heat conductivities of the three substances 
equal to 256, 36, and 0.64 B hr”^ ft“^ F”’- respectively, and in all three cases 
assume == 1.44 B hr”^ ft"'^ F'-h 
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dl = Bq ■" 


Solution: Substituting a. = L in Eq. IX-12 yields the temperature at the 
end of the rod ^ \ 

dl = ^ 1 + e2>"i/ 

This can be simplified by reducing to a common denominator: 

g-mt(l + e”^(l + 

9l = do' (1 + e^”^) 

2(gmt 4- e-”^) pX-13] 

^mL 4- 2 + e-2”^ e”^ + 

flo = to - <o = 130 F 
i = 9/12 = 0.75 ft 
'yn = " n / hC /hA. 
h = 1.44 B hr”^ iirVE~’^ 

cA 




r/3^Y 

t \12/ 


IIX-131 


k = 64 ft-i 

A ^J± 

12 

It is necessary to convert aU inch units into foot units smce foot is 
usS i a unit length in h and fc. It may be shown that the dimension of 
m is ft-^ so that the exponents ma and mi m the above equations beconie 
ZlioSess As mentioned in Sect. VII-2 the exponents of a physi- 
Sr^d mathematically exact equation must be dimensionless. 

With the above values 

I l.44(6§ ^ 


and by substitution in Eq. IX-13 for case 


(o)' fc = 256 


(5) fc = 36 


= — = 0.6ft-i 
16 

^ = ^ = 1.6ft-‘ 

D 


wJj “ 0.45 


mL = 1.2 


(c) h - 0.64 


mL = 9 
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Furthermore : 

(a) = 1.568 

• ( 5 ) gwiL - gi.2 = 3.285 

(c) = e® = 8104 


and by substitution in Eq. IX-12 

(a) Ol = 130 — — = llSFand^L = + ia 


1.568 + 


1.568 


118 + 80 = 198F 


(6) Ol - 130 = = 72.5 F and II = 72.5 + 80 = 152.5 F 


3.285 + 


3.285 


(c) di = 130 0.032 F and ti = 0.03 + 80 = 80.03 F 


8104 + 


8104 


From these computations it is seen that the temperature of the free end of 
the silver rod is only 12 F below that of the heated end, that of the iron rod is 
almost at the mean between the steam and air temperatures, and that of the 
glass rod is practically equal to the air temperature. 

The temperature distribution in the rods may be calculated by substi- 
tuting different values for x in Eq. IX“12. 

The calculation will be performed for xi = 0.25 ft and X 2 = 0.5 ft. 

(a) = 1.111 

gm(0.5) ^ gO.3 ^ 1^244 


g2mL = gmCl.S) gO.9 ^ 2,460 


/ 


Bi = 130 


1.111 


■ + 


V 2.460 


1.111 
1 + 2.460 


= 130 


/ 0.9 
b-407 


h = 205.1 F 


Mli \ = i 25 .ir 
3.460/ 


g = 130 )= 130(0.571+0.360) = 121.0 F 

\l.407 3.460/ 

k = 201.0 F 
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(J) g«(0.26) = ^0.4 == 1.492 

^m(0.5) = gO.8 = 2.226 
gimL = ^2.4 = 11.02 


. iao | '■‘® + 

' , 1 + 11.02 

i i T" 

11.02 


0.671 1.492 

1.091 12.02 


^6 2_^6 

\L091 12.02, 


ti = 176.2 F 


= 130(0.412+0.185)- 77.5 F 
t 2 = 157.5 F 


(c) = 20.09 

^m(0.5) ~ g6 _ 4Q3 5 

g2mL ^ gl8 ^ 65,700,000 


h = 86.5 F 


,, . 130 = 130(0.00248 + 0) - 0.32 F 

k = 80.3 F 

The temperature distributions, according to these figures, are shown 
in Fig. IX-2. All the curves become tangent to the horizontal axis, 
which is true only for a perfectly insulated pad P as showm in Fig. IX-l. 

If the end of the rod was uninsulated, then with the value of h at 
X = L, according to Eq. IX~-2 


(qa)x=L = -kA = -kA 

\dx/x^L \d 


as before. However, for this case it is not equal to zero, and according 
to Eq. 1-3 

{qa)x=L ~ hA{tx=:L ia) ” hAO^ssL 
Equating both expressions one obtains 

./de\ 


0 = 


[IX-14] 
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But from Eq. IX-8 it follows that 

\dxJx=.L 


and from Eq. IX-7 

d:,^L = Me~^^ + 


By substitution in Eq. IX-14 the first boundary condition is obtained 
- Ne^^) = + Ne^^) [ 1 X^ 15 ] 


and the second boundary condition is given by Eq. IX-10 as before. 
Equations IX~10 and 15 may be solved for the two unknown quantities 
M and N, and the resulting expressions are substituted in Eq. IX-7. 

Thus, the procedure is exactly the same as described before; only the 
calculation is more cumbersome. 

The total rate of heat given up by a rod is found by considering that 
all this heat must cross the area A at 3 ; = 0. Therefore it is necessaiy 
to compute 


QxssO — 



[IX-16! 


or, by substitution from Eq. IX-8, using the denotation go instead of 

QxssO 

go = kAm(M - N) [IX-17] 

Example IX-2. Calculate the rate of heat given up by the iron rod in 
Ex. IX-L- 

Solution: Substituting M and N from Eq. IX-11 gives 


M 


i—) 


This can be simplified by employing a common denominator. Then 

g2mL_ Q-imL • ^ _j.„ 


M --N 


1 ^ Q-2mL ^ ^ 


— vu T. — ^ 

_j_ g—wL 

and by substitution in Eq. IX-17: 


^ g>-mL)2 

= Sq tanh(wL) [IX-IS] 


go = 


gwL ^ Q—mL 


= Mm^o tanh(mL) 


[IX-191 
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Substituting the data given for the iron rod: 

go = 36 1.6 (130) tanh 1.2 22.96 tanh 1.2 = 22.96(0.834) 

= 19.1 « 19 B/lir 


Finally, it may be mentioned that the above equations are not 
restricted to cylindric rods, but hold likewise for any projections from a 
heating surface, for instance, fins on a plane plate where A represents 
the area of the cross section through the fin, and C the circumference of 
this area. 


IX-3 Heat Transmission Between Two Fluids Through a Plane Wall 


The combined effect of conduction and convection in heat exchangers 
may be considered in a manner similar to heat conduction through a 
composite wall as shown in Sect. III-S. 



Fig. IX-3. 


Fig. IX-4. 


Counterflow through channels with plane 
separating walls. 


Temperature distribution in the system of 
Fig. IX-3. 


Figure IX-3 represents a cross section through two channels 1 and 4 
with a common plane separating wall, the surfaces of which may be 2 
and 3. Fluids at different temperatures are assumed to flow in the two 
channels, for example, in opposite direction (counter-flow) as indicated 
by the two arrows. For the present, the length of the heat exchanger is 
considered such that the temperature of each fluid may be assumed to 
be constant. Let the temperatures of the fluids be ti and ^4 (see Fig. 
IX-4), and the surface temperatures of the separating wall h and tz, 
its thermal condi|.ctivity And the film coeflScients of heat transfer hi 2 
and /i 34 respectively. Similarly as in heat conduction through a com- 
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posite wall, the rate of heat flowing perpendicularly to the surfaces l 
and 2 through an area A is 

q = hi2A(ti — h) — k2zA ^ . = h^iAitz — t^) 

From this: 


k-h 


g 

/2-12-4. 


h h 


g(Ax)23 


fe — ^4 


/134A 


By adding these three temperature differences it follows that 

-1 


q r 1 
=-l — 


Aik 


, (ha:)23 , 1 1 

“I ; r 


12 


or 


Q = 


A(ti h) 

1 ^ (Aa:)23 ^ 1 

hi2 ’ k2S hz4: 


^23 

= TJA (i^i — ^4) 


[IX-20] 


which corresponds to Eq. III-3. For simplification a new symbol 

1 

hi% 


U - 


(Aa:)23 ^ J_ 

^23 ^34 


[IX-211 


has been introduced. The magnitude U is callec the overall coefficient 
of heat transfer. Its physical dimensions are the t ame as those of the 
film coefficient of heat transfer. 

Since in most practical cases the temperatures ai; d ^4 of the two 
fluids are known or may be measured or calculated easily, whereas the 
measurement of the surface temperatures and k is difficult or even 
impossible, the overall coefficient has considerable practical importance. 


Example IX-3. A heat exchanger w: ^1 consists of a copper plate 5^ 
in. thick. If the two surface coefficients am 480 and 1250 B hr”^ ft~^ F"^ 
respectively, calculate the overall heat transfer coefficient. Assume that h 
for the given problem is 220 B hr“^ ft”^ 

Solution: 


U = 



0.375/12 


+ 


T 

1250 


1 

0.00302 


331 B hr"i ft”2 F-i 


220 
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As in heat conduction through a composite wall (Sect. III-4) heat resist- 
ances may be used. Equation IX-20 then assumes the simple form 


2 = 


^1 ” h 

Ei 2 + E23 + -S34 


p:X-20a] 


which corresponds to Eq. III-3a. 

IX-4 Heat Transmission Between Two Fluids Through a Cylindric 
Wall 

In the case of a double tube heat exchanger, as represented in Fig. 
IX-~5, the procedure is similar to that with thermal conduction through 
several layers of insulation on a cylinder (see Sect. III~6). 



Fig. IX-5. Double tube heat exchanger. 


The nomenclature is analogous to that used for Figs. IX-S and 4. 

The flow of heat per unit length of the concentric tubes of Fig. EX-S 
may be represented by 


or 


o' — /112 • 2 rr 2 {h ~ h) 
= 2rr^{h - k) 


^23 * 27r(^2 ”• ^^3) 

In n/H 


By addition 


k — — 


k ”” k “ 


Jii22rr2 

<i 

^2327r 

In rzlr 2 


k — k — 


hi2Trrz 


^ ^hi22'jrT2 


Inr3/r2 ^ 1 \ 

^2327r hs42'7rrs/ 
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or A _ h h 

® “ 1 1 InrsAa ^ 1 

hu^TTT^ hz^ 2 Trz 

[IX-22! 

Referring to a length L of the double pipe let 


q^q'L 

A .2 ~ 2 tT 2 • L 

Az = 2 rrz * L 

[IX-231 

[IX-24] 

[IX-25] 

From Eqs. IX~22 to 25 the following is obtained: 


q = C72-^2 (^i *“ ^4) = UzAzih 4) 

[IX-26] 

Here U2 and Uz are overall coefficients of heat transfer, which are 
related to inner and outer areas A2 and A3, respectively, by the equations 

1 1 A 2 larz/r 2 ^ A2 

hi 2 kzz^TfL hz^Az 

and 

[IX-27] 

As 1 Aslnrs/r-g ^ 1 

hizAz kzz^n'L hz 4 t 

[IX-281 


These may easily be checked. 

Furthermore, it may be considered that, according to Eq. IX-25, 
= 2'n'r3. Therefore the second term in the denominator of Eq. 


IX-28 equals 


rs In ra/r^ 


Referring to the first term, it is seen that 


— / = — • So this term becomes 7-^ > and 

A 2 27rr2 • L r2 /ii2r2 


Uz 


rs ^ rslnra/ra ^ 1 

hizVz k2z hz4. 


[IX-28a] 


Example IX-4. Calculate the overall heat transfer coefficient based on 
the outer area for a copper condenser tube of %-m. outside diameter having 
a wall thickness of 0.1 in. Assume that the inner film coefficient is 280 
B hr~^ ft"^ F~\ the outer film coefficient is 2000 B hr”^ ft"”^ F“b s-i^d the 
thermal conductivity of copper is 200 B hr“^ ft™^ F*”^. 

Solution: From Eq, IX-28(x: 

TJ^ = ^ ^ = 1R4 n T'-l 

0.00488 + 0.0000485 + 0.0005 0.00543 
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lX-5 Types of Heat Exchangers 

Thus far the length L of a heat exchanger was considered such that the 
temperatures of the two fluids could be assumed constant throughout 
this distance. Generally, the temperatures change in the direction of 
the flow and the question arises as to what average value of h — ti will 
take into account the temperature variations. 

In order to consider the mean temperature difference, the tempera- 
ture distribution in the various t 3 p)es of heat exchangers will be con- 
sidered. The subscripts a and b refer to the two ends of the exchanger 
and the prime sign to the heating fluid. 

All heat exchangers may be divided into the following five general 


Class 1. Heat exchangers wherein a fluid at a constant temperatoe 
gives up heat to a colder fluid the temperature of which gradually 



Area 


Area 


Fig. IX--6. Temperature distribution in 
heat exchanger of class 1 (condenser). 


Fig. IX-7. Temperature distribution 
in heat exchanger of class 2 (boiler). 


increases as it flows through the device. The heating fluid can be at 
rest or moving in any direction. An example of this type would be a 
steam condenser. Figure IX-6 indicates the temperature distribution 
in the apparatus. 

Class 2. Devices wherein .a fluid at constant temperature receives 
heat from a warmer fluid the temperature of which decreases as it flows 
through the exchanger. Here the heated fluid can be at rest or moving 
in any direction. A steam boiler serves as an example of this class. 
Figure IX-7 shows the temperature distribution in this type of heat 
exchange apparatus. 

Class S. Parallel flow heat exchangers wherein the fluids flow in the 
same direction and both of them change their temperature. Many 
devices such as water heaters and oil heaters and coolers fall in this 
group. The temperature distribution for this type of apparatus is 
illustrated in Fig. IX-8. 

Class Ji-. Counter-flow heat exchangers wherein the fluids flow in 
directions opposite to one another. This possibly is the most favorable 
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kind of fluid heaters and coolers. A temperature distribution diagrain 
for this type of apparatus is shown in Fig. IX-9. 

Class 5. Cross-flow heat exchangers in which one fluid flows at an 
angle to the second one, as is the case in tube banks. This group 
not be dealt with in the present chapter. 



Area Area 


Fig. IX-8. Temperature distribution Fig. IX-9. Temperature distribution 
in heat exchanger of class 3 (parallel in heat exchanger of class 4 (counter- 
flow). flow). 


IX-6 The Log Mean Temperature Difference 

A relation for the mean temperature difference will be developed whicli 
holds for heat exchangers of any of the first four classes. Refer to 



Area, y 


Fig. IX-10 which consti- 
tutes a schematic tempera- 
ture distribution diagram 
for a parallel-flow heat ex- 
changer. 

The rate of heat trans- 
ferred through area dy 
located distance x from the 
entrance of the heat ex- 
changer is 

dq=U-dy(t' -t) [IX-291 

The rates of heat dqh lost 
by the hot fluid and dqc 
gained by the cold fluid in 


Fig. IX-10. Temperature distribution diagram ps-ssing over area dy are 
for a parallel-flow heat exchanger. equal to the product of the 


rate of mass flow rrih (or me), 
specific heat Ch (or Cc), and temperature change dt^ (or dt), respectively. 
Thus 


dqh = —mhChdt' 


dqc = +mcCcdt 


[IX-30] 

[IX-311 
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Combining these two equations and dropping the subscripts of q 
since the heat gained by the . cold fluid is equal to the heat lost by the 
hot fluid, gives 

d£ -dt= -dq (— + -i- ^ 

WaCa rriccj 

Substituting d(«' - t) for {di' - dt) and a new symbol N defined b 

^ = — +~ 

’UaCa McCc 

the Mowing is obtained: 

d(t' -t) = -dq-N 
Substituting for dq from Eq. IX-29 
dit' - t) 

(.f - t) 

In (t' - 1) 

One boundary condition is 

t' -t = ta' - ta 

for 2 / = 0, that is, at the flow entrance. Substituting this in Eq. IX-33 
yields 

C = In {ta' - to) 

and again substituting in Eq. IX-33 gives 
, t' - t 

\n——=-UNy [IX-34] 

ta 

Another boundary condition is 

t' - % 

for 2/ == A, that is, for the whole heating area of the heat exchanger. 
Substituting this in Eq. IX-34 leads to 


By integrating 


- ~UN‘dy 
= -UNy + C 


1 , ih jj , 

•—In “7 = - UA 

N ta - ta 


[IX-35] 


The rate of heat qa, transferred on the total area A, may also be 
expressed by 

qA = UA{At)^ [IX-36] 

where (At)m represents a mean temperature difference. 
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From Eqs. IX-35 and 36 it follows that 

yt A 1 , ^6 tb 

= UA = — — In —j— 

(AOm N tj - ta 


or 


(AOm = 


qAN 


tb — tb 

In —r 
ta ta 


Substituting for N from Eq. IX-32 gives 


(Af)„ = - 


qA qA 

nihPh rridCc 

, tb — tb 

ta - ta 


But 


qA 


py=A 

dqh = I dq. 

yssQ tjy—o 


or by substitution from Eqs. IX-30 and 31 


qA 


pt—lh 

—mjfik dt' = I +mcCc dt 

and by integration 

qA == -mOh(th - ta) = rricCcitb - ta) 

This substituted in Eq. IX-37 gives 

'^h^'hi^a ib ) 'fTlcCcita 

(^t)m = - 


nihCh 


nicCc 


, th “ th 
^7T 


or 


(A<)m = 


jtg' - tg) - {tb - tb) 

, tJ - ta 


[IX-3' 


[IX-3 


[lX-3 


This is known as the “ log mean temperature difference.” It is oft 
written in the following form: 


(At)a = 


(Af)„ax - (AOn 


In 


(Ai)mai 

(AOmin 


[IX- 
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SOME APPLICATIONS 


In Eq. IX-40(AO (^Omin represent the maximum, and r 

mum temperature differencesintheheatexchanger. This simple veh 
for the log mean temperature difference may be used for the types of 
exchangers discussed. It must be altered for the more complex i 
of heat exchangers, such as the shell and tube devices with any nui 
of passes on the shell side and tube side, and cross-flow exchangers 
different pass arrangements and with mixed and unmixed flow, 
these types correction factors have been worked out (Ref. IX~2) 
which the log mean temperature difference for counter-current 
must be multiplied to give the true differences. 


Example IX-5. A liquid-to-liquid counter-flow heat exchanger is 
to heat a cold fluid from 120 F to 310 F. If the hot fluid enters at 500 F 
leaves at 400 F, calculate the log mean temperature difference for the 
exchanger. 

Solution: 

/ Ki' Af\f\ 1 nA OOA TH J /A A f'AA rtf r\ f r\f\ 

^ \— vxuiu 


Then accordinff to Eo. IX--40 


mrr. = 


280 


lr» : 


IX-7 Some Applications 

The analysis of the behavior of the various types of heat exchangers 
which include heat transfer by convection depends upon the knowledge 
of the film coefficients, the overall heat transfer coefficient, and the tem- 
perature differences. The values for U and h for the solution of s 
particular problem must either be taken from experimental results, 
industrial operation records, or from dimensional correlations. Som€ 
of the latter have been dealt with in Chapter VIII. 

Example IX--6. Calculate the outside tube area for a single-pass steaa 
condenser to handle 731,300 Ib/hr dry saturated steam. ' The inlet steam 
pressure is 1.09 in. Hg and the hot-well temperature is 81.7 F. The cooling 
water enters the inside of the tubes at 61.4 F and leaves at 69.9 F. The 
tubes are of 1-in. outside diameter by 0.902-in. inside diameter, and the tube 
material has a thermal conductivity of 63 B hr"^ ft“^ F“h The average 
water velocity in each tube is 7 ft sec“^. Assume that the steam side film 
coefficient is 1000 B hr'"^ ft“2 F“^. 

■ - Solution: The average water side coefficient has been determined in Ex 
VTTT-i- waa fnund to be A = 1270 B hr”^ ft'^ F^^, 
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Since the values for the individual coefficients are known, the overall heat 
transfer coefficient may be calculated from Eq. IX-28a. 


Vs = 


i 

1 0.5 In 1/0.902 T 

1270(0.902) 12(63) ^ 1000 


516 B ft'-2 


In the further solution of the problem it will be assumed that the satura- 
tion temperature (81.7 F) corresponding to the pressure of 1.09 in. Hg pre- 
vails in the condenser. Then from Eq. IX-40 


(AO. 


(81.7 - 61.4) - (81.7 - 69.9) 


In 


20.3 

11.8 


15.7 F 


The rate of heat energy required to condense dry saturated steam at 1.09 
in. Hg and 81.7 F is equal to the product of the change in enthalpy and the 
weight rate of steam condensed. Thus 

q = 731,300 (1097.4 -- 49.7) - 765,000,000 B/hr 


The area required according to Eq. IX-26 is 


765,000,000 

516(15.7) 


94,500 ft^ 


Example IX-7. The following data were reported (Ref. IX-1) on a test 
of a shell and tube water-to-water heat exchanger with a total outside tube 
area of 48.1 ft^. The arrangement is shown and the measured temperatures 
are given in Fig. IX~11. The heat exchanger consisted of 98 brass tubes 
having an outside diameter of % in. and a 0.049-in. wall. The shell, 6.06-in. 



Baffles 

Fig. IX-11. Arrangement of a shell and tube heat exchanger with baffles (see 

Ex. IX-7). 


inside diameter, had eleven half-moon baffles located approximately 4.3 in. 
apart. The dimensions and spacing of the baffles were such as to allow a 
minimum outside water passage area between the tubes of approximately 
6.67 sq in. The rate of water flowing through the tubes and shell were 18,210 
and 21,780 lb per hour respectively. 

Calculate the overall heat transfer coefficient (a) from the test data, and 
(5) using the equations which have been established by correlation of data by 
dimensional analysis. 
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Solution: (a) In the derivation of the log mean temperature difference, 
differential elements were considered beginning with Eqs. IX~29 to 31. If 
instead of this the tube sections between any two adjacent baffles are assumed 
to be the elements, the derivation leads to the same result as can be easily 
proved. Therefore, in the present case, counter-flow is assumed to exist 
although this is not exactly true. With this assumption, according to Eq. 
IX-40, 


{AtU = 


(108.1 - 59.5) - (139.7 - 97.2) 
(108.1 - 59.5) 

(139.7 - 97.2) 


46 F 


According to Eq. IX-26, the rate of heat transferred is 
g - 173(48.1)45. 

But according to Eq. IX-38 

g = 18,210(1) (97.2 - 59.5) = 686,000 B/hr 


From the two last equations one obtains C/ = 310 B hr“^ ft”^ F~^. 

(5) The individual coefficient of heat transfer hi for the outside of the tubes 
(heating fluid) will be calculated by means of Eq. VIII-8, considering now the 
flow in the shell as a flow normal to tube banks. 

It will be assumed that the average tube temperature is = 113 F. This 
will be checked again later. 

According to Eq. VIII-7 the average film temperature is 


tf 


123.9 + 


113 - 123,9 


118.4 F 


At this temperature the pertinent thermal properties of the water in consistent 
units are 

p = — slug ft-3 = 1.924 slug ft-2 
32.2 

IX = 3.32(10-8) lb fi;-2br = 3.32(10-8) (3600)8 slug ft-i hi-'- 
= 0.043 slug ft-8 lir-8 

c, = 1 B lb„-i F-i = 1(32.2) B slug-* F-i 

h = 0.367 B lir-i ft-8 F-‘ 

Further oue needs 

Ds = 0.375/12 = 0.03125 ft, and 

= 7620 ft/hr 

6.67/144 
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^ ^ 7620(0.03125 )_1.924 ^ 

^ li 0.043 

k 0.367 

Substituting these dimensionless groups in Eq. VIII-8 


{Nu) = ^ = 0.33(10630) <>'®(3.76) 1/3 ^ 133^ 

k 


and from this 




The individual cQeffi.cient hu for the inside of any of the 98 tubes (heated 
fluid) will be determined from Eq. VIII-2. 

At the average temperature of 78.3 F the thermal properties of water ar( 

p - slug ft“3 =; 1.94 slug ft~3 
o2.2 

M = 5.1(I0-») lb ft-® hr = 5.1(10-9)(3600)® slug ft-®hr-‘ 

= 0.066 slug ft-® hr-® 

Cj, = 1 B Ibm-® F-® = 32.2 B slug-® F-® 
fc = 0.348Bhr-‘ft-®F-® 

Further 

Di = 0.277/12 = 0.0231 ft 
18,210/62.3 

^^._mo(0m«i«. 4830 

n 0.066 

iS . owaye) 

^ ^ k 0.348 

By substitution in Eq. VIII-2 


{Nu) = 0.023(4830)o-8(6,09)<»-4 = 41.95 

k 


and from this 


, 41.95(0.348) 

'*'*“■”6:^31 632Bhr-®ft-®F-‘: 


Substitution in Eq. IX-28 gives the overall heat transfer -coefScient, 
sume 58 B hr-i ft“i F~i for the thermal conductivity of brass. Then f: 
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Eq. IX~28o 

0.03125 1 0.03125 1 0.03125 1 

0.0231 632 '^ 2 58 0.0231 1570 


^ 0.00214 + 0.0000815 + 0.00064 


350Bhr-ift-“F-i 


This value deviates from the actual value by 11.4 per cent of the latter. 
The deviation is not too great owing to the simplifying assumptions which 
had to be made. Considering that the three terms in the denominator of the 
last equation are proportional to the three heat resistances involved (see Eq. 
IX-20a) it is observed that the tube wall presents only about 3 per cent of 
the total resistance. Therefore, of the average temperature drop between 
the two fluids of about 45 F, only about 1.3 F is due to the heat resistance 
of the wall. Since ts is the average temperature of the outer surface of the 
tubes, that for the inner surface would be ^3 ~ 1-3. Then equating the heat 
flow between the heating water and the wall, and that between the wall and 
the heated water 

hziAsih - 123.9) = hnA2[7S,Z ~ (k - 1.3)] 
or 

1570(0.03125)7rL(^3 - 123.9) = 632(0.0231)7rL[78.3 - (^3 - 1.3)] 

Evaluating gives tz == 113.5 F which is in excellent agreement with the value 
113, assumed at the beginning. If there had been a considerable discrepancy, 
the calculation would have to be repeated with another estimate of tz. 

IX-8 Heat Transfer from a Wall at Uniform Temperature Suddenly 
Brought in Contact with a Medium at DMerent Temperature 

In Sect, IV-3 a wall with an initial temperature of ti whose surface 
temperature suddenly changed to ts and remained constant thereafter 
was discussed. A wall at U the sides of which are suddenly exposed to a 
medium at will now be considered. In this case the surfaces do not 
immediately take on the temperature of the surroundings owing to the 
insulating effect of the fluid film next to the surfaces. As it was in the 
case of the thick slab, the complete solution of the problem is beyond the 
scope of this text. However, the final results are available as- charts, for 
instance, those published by G^roeber (Ref, IX~3). These are based 
on the following relations : ::: ; 

/ 4o't 

= ^0 + (ti — [IX--41] 

t, = k + iu- to)^c ^ [IX-42] 

. • Q = pcpAL (ti ^ toW H [IX-43] 
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The symbols used are: 

L = thickness of the wall. 

A = area of the wall, 
p = density 
Cp — specific heat 

h — thermal conductivity of the material of the wall. 

h 

a = — = thermal diffusivity 
pc 

ti = initial uniform temperature of the wall. 

to = constant outer temperature, that is, the temperature of the sur- 
rounding medium with which the two sides of the wall suddenly come in 
contact. 

r = time elapsed after the wall comes in contact with the medium, 
ts = surface temperature of the wall at the time r. 
to = center temperature, that is, the temperature at the median plane 
of the wall at time r. 

h — heat transfer coefficient. 

Q = heat energy released by the wall in the tinie interval r. 

functions represented in Fig. IX-12 with the dimensionless groups 
hL/2 k and 4:aT/L^ as abscissa and parameter respectively. 

It may be mentioned that with infinitely increasing r the function ^ 
approaches the value 1, so that according to Eq. IX-43 the heat released 



hL 

2k 

Fig. IX-12. The functions and in Eqs. IX-41 to 43. 
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in infinite time becomes 

Qca pCpAlj(ti ^o) 

This is easily seen because AL is the volume, pAL the mass, and pc^AL 
is the heat capacity of the wall and {ti the original temperature 
difference. 

Example IX-8. A wall 2 ft thick at an average initial temperature of 100 F 
is exposed on both sides to hot gases at 1000 F. If yfc = 8 B hr~^ ft”^ F“^ 
p = 162 Ibm ft Cp == 0.3 B Ib^ ^ F and h = 1.6 B hr”^ ft“^ F"“^, calculate 
(a) the temperature at the surface, (b) the temperature at the center plane 
of the wall, and (c) the heat absorbed by 10 sq ft of the wall 30 hours and 
24 minutes after the gases come in contact with the wall.* 

Solution: 


and 


From Fig. IX-12 the following values were obtained: 

- 0.37, = 0.41, ^ = 0.62 

From Eq. IX--41 

t, = 1000 + (100 - 1000)0.37 = 667 F 

From Eq. IX-42 

1000 + (100 - 1000)0.41 = 631 F 

From Eq. IX-43 

Q = (162)0.3 (10) 2(100 - 1000)0.62 = -541,000 B 

The negative sign shows that the heat is not released but absorbed 
and stored by the wall. 

* As mentioned in Sect. 1-3, the density, related to pound mass, has the same 
numerical value as the specific weight related to pound force. Therefore y = 162 
lb ft~^ could have been given as well in the above example. But since the ordinary 
numerical values of Cp are related to pound mass, one must use p. The student 
should be encouraged to read, once more, the remarks in Sect. IV-1 concerning 

Cp = pcp. 


k 




8 


pcp 162(0.3) 

4ar 4(0.1645)30.4 


= 0.1645 ft^hr- 


12 


2 ^ 


= 5 


2k 


1 - 6 ( 2 ) 

2 ( 8 ) 


= 0.2 
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PROBLEMS 

IX-l. For the conditions given in Ex, IX-1, calculate the rate of heat given up 
by the silver rod and the glass rod. 

IX~2, For the conditions given in Ex. IX-1 except that the free end of the rod is 
uninsulated, show the temperature distribution in the axial direction by a graph for 
each of the three rods, and calculate the amount of heat given up. What fraction of 
the latter is lost through the free end? 

IX-3. Calculate the overall heat transfer coefficient based on the outer area for 
a condenser tube (ib = 60 B hr"^ F“^) having an outside and inside diameter of 
1 in. and 0.902 in. respectively. Assume that the steam side and water side heat trans- 
fer coefficients are 1000 B hr'^ ft“^ F^^ and 800 B hr“^ ft~- F“^ respectively. 

IX-4. Calculate the overall heat transfer coefficient based on the inside area for a 
steel tube {k = 58 B hr“^ ft"^ F"^) having an outside diameter of 2.375 in. and an 
inside diameter of 2.0 in. Assume that steam is condensing on the outside of the tube 
and that the individual coefficient for the steam side is 900 B hr"'^ ft“^ F"\ The 
coefficient for the inside is 1100 B hr-^ ft“~^ F'h 

IX-5. 200,000 Ib/hr of water are heated in a parallel-flow heat exchanger from 
230 F to 350 F. Hot gases used for heating the water enter the exchanger at 700 F 
and leave at 400 F. Determine the overall heat transfer coefficient if the total 
surface area is 20,000 ftl 

IX--6. A parallel-flow heat exchanger is to be designed to heat 10,000 Ib/hr of 
water from 70 F to 110 F with steam condensing at 250 F on the outside of the tubes. 
Tubes of l-in. outside diameter and 0.9-in. inside diameter^ 8 ft. long are available 
for use. If the average steam side coefficient is 1000 B hr"”^ ff“^ F"^ and the thermal 
conductivity of the tube metal is 64 B hr'“^ ft”^ F“\ calculate the required number of 
tubes per pass and the number of passes. Assume that the velocity of the water enter- 
ing each tube is equal to 1.5 ft/sec. 

IX-7. Determine the heat absorbed by a waE 1 ft thick initially at 80 F after the 
surfaces have been exposed to air at a temperature of 180 F for 3.5 hours, ib = 6 
B hr-i fr^ F“\ = 12 B hr“^ fr^ F~^ Cp = 0.28 B lbm“^ F“\ and p = 150 Ibm ft"®. 

IX-8. For the conditions given in Problem IX-7 calculate the temperatures at 
the center and surface of the waU. 
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CHAPTER X 

HEAT TRANSFER IN CONDENSING AND BOILING 

X-1 Change of Physical Properties in Changes of Phase 

When convection was dealt with in the former chapters, the assump- 
tion was made, although not explicitly stated, that the chemical com- 
position and the phase of the fluid under consideration remained 
unchanged. In this chapter cases wherein a change of phase takes 
place will be briefly considered. The treatment will be restricted to 
condensation and boiling, processes of far-reaching importance in 
mechanical and chemical engineering. In these processes a sudden 
change from the gaseous to the liquid phase or the reverse change 
occurs, and by this the density, viscosity, specific heat, and thermal 
conductivity of the fluid change and heat energy becomes free or is 
absorbed. It is obvious that these changes will exert a tremendous 
influence upon the heat transfer which is connected with the process. 

X-2 Condensation 

Two types of condensation are known; one called dropwise, and the 
other film condensation. The former consists of the formation of 
drops on the cold surface which grow until removed from the surface 
under the influence of gravity and vapor friction. In film condensation 
a thin layer of condensate covers the cooling surface, increases in thick- 
ness, and flows down and off under the influence of the same forces as 
mentioned above. ' ' 

In either type of condensation the heat of condensation which be- 
comes free must be carried from the vapor to the cooling surface. 
Applying Eq. 1-3, Ai is the temperature difference between the vapor 
and the solid surface, and is a film coeflScient of heat transfer which, 
however, only deserves this name for film condensation. Coefficients 
for dropwise condensation are considerably higher than are those for 
film condensation. For steam the film coefficient for dropwise con- 
densation has been observed to be fifteen to twenty times larger than 
for film condensation (Ref. X-1). This, however, relates only to the 
individual heat transfer on the steam side. The overall coefficient will 
increase much less because according to Eq. IX-21 it is controlled by 
the smallest individual coefficient, which for this case would be on the 
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side of the cooling liquid. Dropwise condensation is unstable and, 
until this disadvantage has been overcome, it is unsuited in practice 
for condenser design work. It is a matter of interest, however, to know^ 
something about these two types of condensation (Ref. X-2). 

Film condensation occurs when clean steam condenses on a clean 
polished or rough surface. 

If the surface is contaminated by materials such as mercaptans or 
fatty acids, dropwise condensation of steam may occur. 

Often both types of condensation occur at the same time on the same 
surface; this process is known as mixed condensation. 

W. Nusselt developed in 1916 (Ref. X-3) a simple theory of film 
condensation, the main features of which are given as follows: 

A vapor condensing in a film on a cooling surface flows on the surface 
under the influence of gravity and vapor friction, but is retarded by the 
viscous forces of the liquid. The heat of condensation passes through 
the film from the condensing vapor to the wall. Because the thermal 
conductivity of liquids is small, the film, thin as it may be, presents a 
rather great resistance to the flow of heat, and, therefore, an appreciable 
temperature drop exists across the film. The temperature of the sur- 
face of the film in contact with the wall is equal to that of the wall, and 
the surface of the film in contact with the vapor is assumed to be at the 
saturation temperature. By combining the laws of laminar flow of a 
fluid and of heat conduction through it, Nusselt arrived at the follow- 
ing relation: 


h 


= C 


Jjjx • At/ 


Cx-i] 


wherein, using foot, hour, and slug units 

h == the average coefficient of heat transfer, taken for the 
whole condensing surface, in B hr“^ ft”^ ; 

g = the gravitational constant 418 (10® )ft hr"^; 

I — the latent heat of evaporation in B/slug; 

L = the height of a vertical wall or the outer diameter of a 
horizontal tube on which the condensation occurs; 
p, p, and k = respectively, the density, dynamic viscosity, and ther- 
mal conductivity of the condensate; 

C ~ 0.943 for vertical walls; 

C = 0.725 for horizontal tubes. 


This equation is valid for arbitrary vapors which are practically at 
rest, that is, no artificial means are used to move the vapor. For forced 
convection of condensing vapors see Ref. X-3. 
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Example X-1. Determine the heat transfer coefficient for stea i at 142 F 
condensing on the outside of a horizontal l-in. diameter cylinder r;.t 138 F. 

Solution: According to the steam tables (Ref. VIII-2) and Tal ^s.II'-2 and 
YI-1, the physical properties taken at the average temperature of 140 F are 

y = 61.3 Ib/ft^ corresponding to p ‘ = 1.90 slug/ft^, 

32.2 

I - 1014 B/lb or I = 1014(32.2) - 32,600 B/slug, 
k = 0.377 B hr-i ft-^ F'^ 
fi - 2.73(10-9) lb ft-^hr or jLt = 2.73(10-9) (3600)2 
= 0.0354 slug ft”^ hr-^. 

By substitution in Eq. X-1 

_ n >70,= r (418)10^(1.90)2 32,600(0.377)^ ^/^ 

L (1/12) 0.0354(142 - 138) J 

0 . 725(3865) = 2800 B hr-^ fr^ F-^ 


X-3 Boiling 

A glance into a kettle of boiling water shows that boiling is a process 
of heat transfer by pronounced convection. Until recently, very little 
was known as to the actual nature of this process, notwithstanding the 
fact that engineers have been designing and constructing boilers for 
more than a hundred years. There are two kinds of vaporization 
analogous to the two types of condensation, the one is called film boiling, 
the other nuclear boiling. The former corresponds to fihn condensa- 
tion, but usually does not occur in practice. It is observed when heating 
is so intense that a vapor layer is formed on the heating surface. Such 
a layer, owing to the low heat conductivity of gases, forms a resistance 
to the heat flow, and by this the temperature of the heating surface 
increases rapidly. At times, this type of boiling becomes dangerous. 
The phenomenon of water drops dancing on a very hot plate without 
evaporating is due to this film effect. 

The usual kind of vaporization is nuclear boiling, that is, boiling which 
starts from nuclei, tiny cells of vaporization, such as air bubbles which 
were dissolved in the liquid or absorbed on the heating surface before 
the heating took place. Such gas cells generally originate on the small 
natural roughnesses of heating surfaces, and for this reason boiling always 
begins on these surfaces. 

Investigations of M. Jakob, performed in the last decade with differ- 
ent co-workers (Ref. X-3 and 4), have revealed the mechanism of 
nuclear boiling as follows. 
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Fron (le solid heating surface the heat energy is transmitted by 
convectioi to the liquid in contact with that surface and then from the 
liquid to \e surface of the tiny bubbles which form on the solid sur. 
face. At t4e bubble surfaces more liquid is vaporized by the incoming 
heat flow, and the bubble grows until by its own buoyancy it breaks 
off from the surface and rises in the liquid. In this movement, more 
evaporation takes place on the bubble surface; it increases in rising 
and finally blows up at the level which separates the liquid from the 
vapor space above it. This mechanism is only possible if temperature 
differences exist by which the heat is driven in the described way. 
Indeed Jakob experiments have proved that the temperature of a 
boiling liquid is always slightly higher than the saturation temperature 
which exists inside the vapor bubbles and above the liquid level. Fur- 
thermore, the heating surface is warmer than the liquid. 

As a consequence of the fact that heat energy in a boiling liquid 
generally makes a detour from the heating surface over the liquid to 
the vapor bubbles, the experimental results were represented by formu- 
las which are identical with those dealt with in Chapters VI and VII 
for free convection between a solid surface and a fluid. This is valid 
up to a certain limit which is at (Gr) (Pr) ~ 10^ to 10®. Below this 
limit the movement of the vapor bubbles is of secondary influence as 
compared to the convection between the heating surface and the liquid. 
Above this limit the columns of rising vapor bubbles act like liquid 
stirrers and forced convection takes place with much more intensive heat 
transfer and vaporization. Research in this range is still in progress. 
Some of the more recent results will be found in Ref. X-3, 5, 6, and 7. 

Denoting the rate of heat flow per unit heating surface by g" and the 
unit value of this magnitude by qo" = 1, the following formulas are 
recommended for water boiling at standard atmospheric pressure ps 
under free convection conditions. In all of these equations the units 
for h are B hr~^ ft“^ 

For horizontal heating surfaces in wide vessels up to = 5000 

Bhr-^r^: 



[X-21 

from q” = 6000 to 76,000 B hr-^ 

/ // y/4 

A.. 0.64 (L) 

tK-3] 

For the inside of vertical tubes: 

h = 1.25 
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For vertical heating surfaces in wide vessels up to g'' = 1000 
B hr""^ ft""^: 

ft, = 60(^) IX-5] 

\Qo / 

from g" = 1000 to 20,000 B 

fa" 

A. = 0.7 [X-G] 

The influence of pressure p can be considered approximately by using 
the formula 

h-=liA— [X-7] 


wherein the subscript s relates to standard atmospheric conditions.* 

It may seem strange that in the above equations is divided by a 
magnitude which is always equal to 1. This was done to retain 
the dimensional soundness of the equations, q^^ /qo' is independent of 
the system of units used. Therefore, the constant factor of the equa- 
tions, for instance, 50 in Eq. X-5, is the value of the coefficient h for 
q" = qo" = 1 B hr”^ ff”^. This. formal statement, however, does not 
mean that the equation is valid down to = 1. If, for instance, it is 
desired to convert the equation into metric units, it is only necessary 
to convert the constant factor into kcal hr”^ m“^ C”^ and to take g'' 
in kcal hr""^ m""^ units, so that Eq. X-5 becomes 


, 50 

^ 0.205 Vo"/ 



tX-8] 


in the metric system. Using — 50(g'0^^^ instead of Eq. X-5 and 
then simply converting and q' into metric units would lead to 
erroneous results as can easily be checked. 


* Eqs. X-2, 5, and 6 have been set up by M. Jakob based on unpublished correla- 
tions of his own and some other experiments (see Ref. X-3 and 4), particularly those 
of Insinger and Bliss (Ref. X-6). Equation X-3 has been given by Fritz (Ref. X-8) 
and is based on experiments of Jakob and Fritz, and Jakob and Linke. Equation 
X-4 is a rough estimate by Jakob (Ref. X-5). According to the latter paper, the 
influence of pressure from p = 1.8 to p = 14.7 Ib/sq in. could be considered approxi- 

/ 

mately by multiplying h for atmospheric pressure ps with ( — ) and between 

, \Ps/ 


p = 14.7 and p == 226 Ib/sq in. by multiplying with | 




J 


Equation X-7 is due 


to Bonilla and Perry (Ref. X-7). The exponent M their equation is just the mean 
of the exponents and 
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Because 

q" = 1i{M) [X_ 9 ] 

by definition where M is the tenaperature difference between the heating 
surface and the boiling liquid, it is easy to express A as a function of M 
instead of as has been done in Eqs. X-2 to 8, 

PROBLEMS 

X-1. Compute the heat transfer coefficient for steam at 160 F condensing on the 
outside of a horizontal 3-in. diameter cylinder at 154 F. 

X-2. Determine the heat transfer coefficient for steam at 142 F condensing on a 
vertical plate 6 in. high at 138 F. 

X-3. Calculate the approximate time required to condense 500 lb of steam at 
150 F on a horizontal cylinder 1 in. in diameter and 3 ft long at 144 F. 

X-4. Determine the value of the heat transfer coefficient h in kcal hr""^ CT'i 
units for steam condensing at 95 C on a horizontal pipe 4 cm in diameter. Assume 
that the temperature difference is 10 C. 

■ X-5. If the rate of heat input to the horizontal surface of an apparatus used for 
evaporating water under standard atmospheric conditions, is 2800 B hr""^ ft'-, 
calculate the heat transfer coefficient. 

X-6. Solve Problem X-5 if the heat input had been 60,000 B hr“^ ft”'^. 

X-7. The heat input based on the outer surface area of a vertical thin wall pipe 
2 in. in diameter is 3000 B/hr for 5 ft of pipe. If water flowing slowly through the 
tube is being evaporated under standard atmospheric conditions, estimate the heat 
transfer coefficient for the boiling liquid. 

X-8. Calculate the heat transfer coefficient between a vertical plate and water 
boiling at standard atmospheric pressure if the heat input is 600 B hr~^ ft”^. 

X-9. The heat input to a vertical wall of a device for boiling water is 15,000 
B hr“^ ft"^. If the pressure in the apparatus is five atmospheres absolute, calculate 
the approximate heat transfer coefficient between the wall and the boiling liquid. 

X-10. If the heat transfer coefficient between a vertical plate and boiling water 
is 500 B hr”"^ ft“^ F“^, estimate the rate of heat input to the plate, assuming standard 
atmospheric conditions. 

X-11. Calculate the heat transfer coefficient in metric units between a vertical 
tube and boiling water if the rate of heat input is 8000 kcal hr“^ m“^, assuming 
standard atmospheric conditions. 
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CHAPTER XI 

HEAT TRANSFER BY RADIATION 

XI-1 Heat Radiation — A Type of Wave Motion 

Whenever heat flows from one part of a body to another by con- 
duction the intervening material of the body becomes heated, that is, 
the temperature of the body will vary from the hot side to the cold 
side according to a definite relation, depending upon the position in the 
body. Radiant heat transfer differs from heat flow by conduction in 
that the medium through which the flow takes place does not become 
heated. The earth receives radiant energy from the sun, the radiation 
passing through the intervening cold empty space. The warmth one 
feels from an open camp fire on a chilly autumn evening is due chiefly 
to the radiation from the fire to the body. This energy is transferred 
through the cool night air without heating it up appreciably. 

The transfer mechanism of radiant heat energy is complicated, but a 
general understanding may be had by considering the heat transfer 
from the sun to the earth. The thermal energy of the sun is first con- 
verted into an electro-magnetic wave motion called radiation which 
travels from the sun to the earth in about eight minutes and is there 
absorbed, raising the temperature of the absorbing medium. Thus, 
the mechanism of energy transfer by radiation is composed of three 
distinct components: first, the conversion of the thermal energy of the 
hot source into an electro-magnetic wave motion; second, the passage 
of the wave motion through the intervening space; and, third, the 
reconversion of the wave motion into thermal energy by absorption at 
the cold body. 

Radiant energy is the same type of wave motion as radio waves, 
X rays, and light waves except for the wavelength. It is apparent that 
radiant heat energy is governed by exactly the same laws as light: it 
travels in straight lines, obeys the laws of reflection, suffers refraction, 
may be polarized, and is weakened with the inverse square of the radial 
distance from the source of radiation. The velocities of light and 
radiant energy waves are the same, that is, about 186,000 miles per 
second. For purposes of comparison some well-known electro-magnetic 
waves are given in Table XI-1. 
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THE CONCEPT OF A PERFECT BLACK BODY 

TABLE XI-1 
Electro-Magi^etic Waves 
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Name 

Cosmic rays 
Gamma rays 
X rays 

Ultraviolet rays 
Visible or light rays 
Infrared or heat rays 
Radio 

*1 micron — 10-* meter. 


Wavelength Range 
in Microns* 
up to 1(10"*®) 

1(10"®) to 140(10“®) 
6(10“®) to 100,000(10"®) 
0.014 to 0.4 
0.4 to O.S 
0.8 to 400 

10(10®) to 30,000(10®) 


XI-2 The Concept of a Perfect Black Body 

When radiant energy falls on a body, part may be absorbed, part 
reflected, and the remainder transmitted through the body. In mathe- 
matical form, 

^ + p + r = 1 EXW] 

where a = absorptivity or the fraction of the total energy absorbed, 
p = reflectivity or the fraction of the total energy reflected, 

T == transmissivity or the fraction of the total energy transmitted 
through the body. 

For the majority of opaque solid materials encountered in engineer- 
ing, except for extremely thin layers, practically none of the radiant 
energy is transmitted through the body. If the discussion is limited to 
opaque bodies Eq. XI-1 becomes 

a + P = 1 lKl-2] 

An arrangement which will absorb ail the radiant energy at all wave- 
lengths and reflect none is called a perfect black body. Actually no 
material with a = 1 and p = 0 exists. Even the blackest surfaces 
occurring in nature still have a reflectivity of about 1 per cent (p = 0.01). 
The physicist G. Kirchhoff, however, conceived the followmg possi- 
bility of making a practically perfect black body. If a hollow body is 
provided with only one very small opening, and is held at uniform 
temperature, then any beam of radiation entering by the hole is partly 
absorbed, and partly reflected inside. The reflected radiation will not 
find the outlet, but will fall again on the inside waU. There it will be 
only partly reflected and so on. By such a sequence of reflections the 
entering radiation will be weakened so much that almost no part of it 
will leave the hole. Thus the area of the hole is like a perfectly ab- 
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sorbing surface, and an arrangement of this kind will act just as a perfect 
black body. It may be considered a measure by which the absorptivity 
of any substances may be determined. 

XI-3 Planck’s Law of Monochromatic Radiation of a Black Body 

All substances emit radiation, the quality and quantity depending 
upon the absolute temperature* and the properties of the material com- 
posing the radiating body. 



Wave Length in microns 
(1 micron = 10”^^ cm) ' 

Fig. XI-1. Monochromatic intensity of radiation for a black body at various absolute 
temperatures (Planck’s law). 


* As far as radiation depends on the temperature of the radiating body, it is called 
temperature radiation. There are other kinds of radiation which, however, are not 
of the general importance of temperature radiation and, therefore, will not be dealt 
with in this book. 
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The total radiant energy emitted by a surface per unit area per unit 
time will be called the total density of emission and denoted by the 
symbol E. It can be considered as an integral of differential values, 
dE\, according to the equation 

J ^X=oo pk=OS 

' dEx== Jx-dx [XI-3] 

x=o */x=o 


taken over the whole range of wavelengths from X = 0 to X = oo. By 
this equation another magnitude, /x, is defined which is called the 
monochromatic intensity of radiation, the subscript relating to the 
wavelength X of the radiation. For a black body, the respective sym- 
bols are Eh and Ih\. The latter magnitude depends on the temperature 
of the radiating surface and the wavelength of radiation according to 
Planck^s law of monochromatic radiation. This is represented graphi- 
cally in Fig. XI-1, with the wavelengths as abscissas and the absolute 
temperature as parameter of the single curves. Eb and hi are con- 
nected by the relation 




|X= 

x=o 


Ib\d\ 


[XI-3a] 


The areas between the curves and the axis of abscissas from X = 0 to 
X = 00 give the respective total flux densities Eb in watt cm“^]. 

It may be mentioned that a few gases such as carbon dioxide and 
steam in thick layers radiate like a black body, but not at aU wave- 
lengths, as shown in Fig. XI-l, but only in some narrow wavelength 
intervals, called bands of gas radiation. The treatment of gas radiation 
is beyond the scope of this text. 


XI~4 Kirchhoff’s Law of Radiation 

In order to develop a relation between the emissive and absorptive 
power of a surface, first consider a large heated enclosure wherein the 
radiation is between X + dX and X. After the equilibrium temperature 
has been established, it is permissible to assume that the black body 
condition exists within the enclosure. A second body, having the 
absorptivity and emitting radiation with the intensity h, is placed 
within the cavity of the large body. When ,the equilibrium condition 
has been reestablished, the amount' of energy absorbed by the small 
body will be equal to the energy emitted in the same time. Let dq 
represent the rate of energy at wavelength X Striking the enclosed body. 
The rate of absorbed energy will be equal to ax * dq. By definition the 
rate of emitted energy will be equal to Ail\ • dX where Axis the area of 
the enclosed body. Since the two rates of energy are equal 

oc)^‘ dq = Ail\ dX 


[XI-4I 
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For a black body this equation would become 

1 • dQ ~ Ai/fex d\ [XI— 4(x] 

From Eqs. XI-4 and 4a 

h = [XI-5] 

Introducing a new term, the emissivity ex of a surface, by the definition 



it follows that 

fcx = oi\ [XI-7] 

Any one of the last three relations represents Kirchhoff’s law. Ex- 
pressed in words, Eq. XI-5 means that the rate of radiant energy, which 
is emitted by a surface at any temperature and in a small range of wave- 
lengths, is found from the known rate of energy which, under the same 
conditions, would be emitted from a black surface, by multiplying with 
the absorptivity. 

Because the absorptivity, according to Eq. XI-2, must lie between 
0 and 1, it is apparent that a black body is the best emitter of energy. 
Equation XI-7 further shows that a good emitter is also a good absorber 
and that emissivity and absorptivity are identical properties of a surface. 

By integrating Eqs. XI-4 and 4a over the whole range of wavelengths 
X = 0 to X = 00 , and considering Eqs. XI-3 and 3a, it is easily seen 
that Kirchhoff^s law holds for the total radiation as well as for mono- 
chromatic radiation. In terms of total radiation it is usually expressed 
in one of the following forms: 

E-=aEb [XI-5a] 

€ = a \Xl-7a] 

Since perfectly black substances are not attainable in engineering 
practice, the equation established for black body radiation must be 
modified if other bodies are to be considered. For convenience the 
emissive power of an ordinary opaque body will be referred to that of a 
black body by the emissivity e. 

XI-5 Stefan-Boltzmann’s Law of Total Radiation 

In 1879 J. Stefan concluded from experimental data that the total 
energy emitted by a black body is proportional to the fourth power of 
the absolute temperature of the body. Five years later L. Boltzmann 
derived the same law from a theoretical thermodynamic standpoint. 
It likewise can be derived from Planck^s law by integration or by plani- 
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metering the areas below the curves in Fig. XI-1 from X = 0 up to 
X = 00 . In honor of its discoverers the law has been called the Stefan- 
Boltzmann law. It may be written either in the form 

Ei = (tT^ [XI-8] 

or by ^ 

= aAT^ pCI-9] 

where A = the radiating area 

T = the absolute temperature and 

0 - = Stefan-Boltzmann’s natural constant. 


In British technical units, according to the best knovm measurements, 
(T = 0.174(10”®) B hr”^ ft”^ 


For convenience Eq. XI-9 may be written as 

„ -0.174 A (i) 


[XI-IO] 


Using this presents the advantage that the term raised to the fourth 
power gives relatively low figures in numerical calculations. Equations 
XI-8, 9, and 10 hold for black surfaces only. 

If the emissivity of a body can be considered as independent of the 
wavelength and temperature, which fortunately is approximately true 
for most non-metallic surfaces in the ordinary temperature range, the 
body is called a gray radiator. For this kind of radiator, instead of 
Eq. XI-9, the following generalization of Stefan-Boltzmann’s law can 


be used 


Qg = igoAT^ 


[XI-11] 


where subscript g means gray radiator and eg is the emissivity of the gray 
surface. Still more general is the form 

q = e<rAT^ Pa-12] 


This formula has been mentioned in Chapter I as Eq. 1-4. It can 
be used for any surfaces; e, however, will not be a constant but a func- 
tion of T in the general case. 


XI-6 The Emissivity or Absorptivity of Different Bodies 

The emissivity or absorptivity of a body has an importance in radi- 
ation «iTT 7 i 1 a.r to that of heat conductivity in conduction. They are 
factors of proportionality in the basic Eqs. 1-4 and 1 respectively. 

Generally, it is difficjjltJf-QSLSSE^i^ estimate or indicate the 
emissivity of iTSuface with an accuracy of a few per cent because the 
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emissivity depends to some extent on the behavior of the surface, par- 
ticularly as far as metallic surfaces are concerned. 

The emissivity of non-metallic bodies does not vary much at ordinary 
temperatures. Some approximate values are given in Table XI-2. 

TABLE XI-2 

Emissivity e of Non-Metallic Bodies at Ordinary Temperatures 
Materials 


Iron oxidft, carbon, oil 

0.80 

Rubber (gray, soft), wood (planed), paper 

0.85 to 0.90 

Roofing paper, enamel, lacquer, porcelain (glazed) 1 

Fused quartz (rough), brick (red, rough) > 

0.91 to 0.94; 

Marble (gray, polished), glass (smooth) J ^ 

Asbestos slate (rough), lampblack-waterglass, ice, wat^ 

0.95 to 0.99 


It is seen that smooth and rough surfaces (like polished marble and 
rough brick) have almost the same emissivity, and that ice and water 
are close to black body conditions, as is lampblack. This, however, 
holds only for low-temperature radiation. According to Fig. XI-l the 
maximum of intensity of radiation shifts to greater wavelengths with 
decreasing temperature. Thus radiation at low temperature is mainly 
long-wave radiation. If short-wave radiation, like that of the sun, 
strikes ice or a white surface the absorptivity is much smaller than when 

TABLE XI-3 

Emissivity € of Polished Metallic Surfaces 


Metal 

100 F 

Temperature 

500 F 

1000 F 

Aluminum 

0.04 

0.05 

0.08 

Copper 

0.04 

0.05 

0.08 

Gold 

0.02 

0.02 

0.03 

Silver 

0.01 

0.02 

0.03 

Steel 

0.07 1 

0.10 

0.14 


it hits a black surface. The good reflection of sunlight by ice or white 
fabrics is well known. It is less well known that the reflection of the 
same bodies is very small for long-wave radiation. So, if a quantity 
of soil is mixed with lampblack and another is mixed with chalk, it may 
be found that the black soil will be 10 to 12 F warmer than the white 
soil when exposed to the radiant energy of the sun, After sunset, how- 
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ever, white soil will cool as fast as black soil if the issuing temperature 
is the same, indicating that there is no appreciable difference in the 
radiating properties of black or white soil for long-wave radiation. 

Metallic smooth surfaces emit very little radiation at ordinary temper- 
ature, and the emissivity increases moderately at higher temperature. 
A few values are given in Table XI-3. 

However, when the surfaces are oxidized, i.e., covered with a chemi- 
cally non-metallic layer, the emissmty increases enormously. 

For instance, heavily oxidized aluminum may have e = 0.85 at 100 F 
and € — 0.50 at 1000 F ; black oxidized copper or cast iron with cast 
skin or oxidized steel may have, e 0.80 at ordinary temperature. 

XI-7 Heat Exchange by Radiation Between Large Parallel Black 
Planes 

Thus far, emission or absorption of a surface has been discussed 
without considering its exchange of radiation energy with other sur- 
faces. In practical cases, how'Cver, there is ahvays such an exchange, 
and this often complicates the situation very much. 

The reason one feels warm when standing in front of an open fireplace 
is because the fireplace radiates more energy to the body than the body 
radiates to the fireplace. Engineeiing problems of radiation must be 
dealt with in a similar way. 

Consider, for instance, two large parallel planes. It may be assumed 
that the surfaces are perfectly black and so large that the influence of 
the edges will be negligible, that is, as though the areas w^ere infinitely 
large. Then obviously all the energy radiated by one plane will be 
received by the other. If the absolute temperature of the warmer plane 
is Ti, then the rate of energy emitted per unit area is 

Em = crTi" 

Likewise, the rate of energy emitted by the colder plane at absolute 
temperature is 

Eh2 ” 

Since all the energy received by either plane is completely absorbed, the 
net interchange between the two surfaces is 

qu’ = Eii - ^52 = - T2^) CXI-13] 

As formerly, q means the rate of heat flow, and the double prime sign 
indicates that unit areas are considered, whereas the double subscript 
(12) relates to the net heat exchange betw^een the surfaces 1 and 2 and 
the direction of the heat flow from 1 to 2. 
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Example XI-1. Calculate the net radiant interchange between unit areas 
of two parallel perfectly black planes, infinite in extent, at temperatures of 
800 F and 1200 F respectively. 

Solution: 

Ti - 460 + 1200 = 1660 R T2 = 460 + 800 = 1260 R 
Substituting in Eq. XI-13 gives 

. ... r/l660V /126OVI 

‘ L\ 100 / \ 100 ) J 

= 8,850 B hr-i ft~ 

XI -8 Heat Exchange by Radiation Between Large Parallel Planes of 
Different Emissivity 

The rate of energy emitted by unit area of plane 1 at absolute 
temperature Ti and emissivity ei is Ei = Of this energy the 

part e 2 Ei = will be absorbed by plane 2. The difference 

El — €2^1 represents the amount reflected by 2 since the transmissivity 
T is assumed to be zero. Part of the energy reflected from plane 2 will 
be absorbed at plane 1 and the remainder reflected back to plane 2 and 
so on. By similar reasoning the density of radiation E 2 = may 
be traced through the various absorptions and reflections (Ref. XI-1). 

The rates of all amounts of energy emitted or reflected minus those 
absorbed per unit area by plane 1 may be represented by the series: 

= ei(TTi^[l — ei(l — €2) — €1(1 — €2)(1 €i)(l — €2) 

— 6l(l — €2)(1 — €1)^(1 — €2)^ — *] 

— €20-T2^[61 + €i(1 — €i)(l — €2) + €1(1 — 6 i)^(l — €2)^ + •] 

or introducing the symbol 

(1 -€i)(l --€2) [XI-14] 

qi 2 ^ = eiaTi^ll — €i(l €2) — €i(l — .62)2; ~ €i(l ■— €2)2;^ — . . .] 

- €2<tT2%1 + €1^ + €i2{2 + . . .] [XI-15] 

But because 25 < 1, 

l+2+2^+2^ + ‘'* = — 

Substituting this in Eq. XI-15 

qi2^ = cicrTi^ri ^1 — €2(rY2^--^^ — 

L I - z J 1 - z 
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?i2"=i . Pa-15a] 

- + 1 

€l €2 


This is the net heat exchange between the planes because it is the excess 
of emission and reflection of plane 1 above its absorption. 

Example XI-2. Calculate the net radiant interchange per square foot for 
two very large planes at temperatures of 1000 F and 600 F respectively. As- 
sume that the emissivity of the hot plane is 0.9 and for the cold plane 0.7. 
Solution: 

Ti = 460 + 1000 = 1460 R 
Tg = 460 -h 600 = 1060 R 

According to Eq. Xl-lSa 

212' ^ (14.6« - 10.6‘) = 3710 B hr-i fr® 

09 ^ 


XI-9- Heat Exchange by Radiation Between an Enclosed Body and the 
Enclosure. 

The radiant interchange between a body 1 having a complete convex 
surface and a concave surface 2 which surrounds it may be determined by 
the following relation, derived first by C. Christiansen in 1883: 


gi2 - - T2^) 



pa-16] 


In this relation the subscript 1 refers to the enclosed body and 2 to 
ihe enclosing surface. Ai and are the respective surface areas. 


Example XI-3. Pentane at —221 F is to be stored in the inner of two 
ioncentric polished brass spheres 9 in. and 12 in. in diameter respectively. If 
he emissivity of the polished brass is 0,03, calculate the total radiant heat in- 
ierchange between the spheres. Neglect the temperature drop through the 
netal and assume that the outer sphere metal temperature is 70 F. 

Solution: 

IL c;\2 

A I == = 47r = .77 ft^ 

A 2 = 47rr2^ = 47r = 3.14 ft® 

Ti = 460 - 221 = 239 R 

Ta - 460 + 70 - 530 R 
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According to Eq. XI-16 . 

qn = (0.174)1.77(2.39* - 5.30*) — r 

- o-:^3+iT4te-V 

= -4.53 B/hr 

The negative sign of the result indicates that heat is flowing from surface 

XI-10 Heat Exchange by Radiation Between a Small Enclosec 
and the Enclosure 

Assuming that the area A 2 of the enclosure is large comparec 
the area Ai of the enclosed body, Eq. XI~16 simplifies to 

Q12 = €i(rAi(Ti^ — ^2^) \X %. 

As in Eq. XI-16 the enclosed body must be convex throughout, an 
enclosure concave throughout. 

Example XI-4. A 2-in. oxidized iron pipe at 300 F passes through a 
in which the surroundings are at a temperature of 80 F. If the emissiv. 
the pipe metal is 0.8, calculate the net interchange of radiant energy pei 
length of pipe. 

Solution: 

Ti = 460 + 300 = 760 R; ^2 = 460 + 80 = 540 R 
The surface area of one foot of 2^^ standard pipe is 
At = X2.375/12 - 0.624 ft^ 

qu' - 0.8(0.174)0.624(7.6^ - 5.4^) = 215.5 B hr^i fr^ 

The prime sign indicates that the result is given per unit length. 

XI-11 General Equations for Heat Exchange by Radiation 

All the relations used in the previous cases are similar to Eq. XI-1 
except for a term which depends upon the emissivities of the material! 
A form of equation for all is 

qi2 = FE<rAi{T,^ - [XI-., 

In this relation Fe is an emissivity factor which shows the devii. 
tion from black body conditions in so far as the emissivities are cor 
cerned. 

For many cases it is necessary to add further a configuration factor 0 . 
area factor Fa which accounts for the fact that not every point of thi 
one area can be connected with every point of the other area by t 
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straight line which does not intercept either of the surfaces. This 
situation is commonly stated by saying that the two radiating sur- 
faces do not see each other completely. Thus the general equation 
becomes 

qi 2 ^ FEFA<rAiiTi^ - T 2 ^) [XM9] 

Consider the two radiating areas dAi and dA 2 in Fig. XI-2 at tempera- 
tures Ti and T 2 . The net interchange of energy if black body con- 
ditions prevail is equal to the energy 
absorbed by dA 2 which represents a 
fraction of all the energy radiated by 
dAi, minus the energy absorbed by dAi 
which represents a fraction of all the 
energy radiated by dA 2 . It may be 
shown that this net interchange of radi- 
ant energy is 


/ cos (j>i • cos 


■dAi- di-a X ^ 

Trr^ ) 

[XI-20] 



This follows from the fact that heat 
radiation like light follows the law of 
the inverse square of the distance. The XI-2. Denotations 
evaluation of q by integration is difficult the configuration of two areas, 
for complex cases. It is beyond the 

scope of this text to show the steps in the solution of the equation. 
For a particular case, however, the final solution will be presented. 

XI~12 Heat Exchange by Radiation Between Equal Parallel and 
Opposite Squares 

For this case Hottel (Ref. XI~2) recommends the use of Eq. XI-19 
with Fb == €162 when the areas are small compared to the distance 
apart, and with 

€l €2 


when the areas are close together.^ By integration of Eq. XI-20 Hottel 
has found values of Fa us given in Table XI-4 where Fa is represented 
as a function of the ratio of the square side S to the distance D. 
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TABLE XI-4 

Aeea Factor for Parallel and Opposite Squares 


1 

D 

Fa 

0.5 

0.07 

1.0 

0.21 

1.5 

0.32 

2.0 

0.42 

2.5 

0.49 

3.0 

0.55 

3.5 

0.59 

4.0 

0.64 

4.5 

0.66 

5.0 

0.69 

5.5 

0.72 

6.0 

0.74 


Example XI-5. Calculate the radiation between the floor (15 ft by 15 ft) 
of a furnace and the roof if the two areas are located 10 ft apart. The floor 
and roof temperatures are 2000 F and 600 F respectively. Assume that the 
emissivity values for the floor and roof are 0.9 and 0.8 respectively. 

'Solution: 

= Ag = 15(15) = 225 ft2 

Ti = 460 + 2000 = 2460 R and Ts = 460 + 600 = 1060 R 


€i€2 = 0.9(0.8) = 0.72 



€l 62 


TakeF^ = 0.73. S/D = 15/10 = 1.5. From Table XI-4, = 0.32. With 

these values, Eq. XI-19 becomes 

qi2 == (0:i7)0:73(0.32)225(24.64 - lO.O^) 

= 3,240,000 B/hr 


PROBLEMS 

XI-1. If 360 British thermal units of radiant energy strike a body having an 
absorptivity of }i and a reflectivity of 34? calculate the heat transmitted through the 
body. 

■ XI-2. Calculate the radiant heat loss passing through a 2-in. diameter opening 

in a furnace door, if the furnace and outside temperatures are 1800 F and 70 F re- 
spectively. Assume black body conditions. 

, ' XI~3. Calculate the radiant heat loss from a 3-in. nominal wrought iron oxidized 
pipe (outside diameter = 3.5 in.)^at 200 F which passes through a room the tempera- 
ture of which is maintained at 76 F. 
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; XI-4. If the heat transfer by radiation from the unpainted side of a house wall to 
the plaster side is 10 B hr~^ ft“^, calculate the emissivity of the unpainted waU. 
Assume that the average surface temperatures of the unpainted side and the plastered 
side are 110 F and 90 F respectively. Emissivity of the plastered wall equals 0.93. 
\. XI-5. If the average rate of radiant energy received and absorbed at the earth’s 
• surface on a clear day is 300 B hr"* ft"-, calculate the approximate temperature of 
the sun. Assume that 100 B hr"* ft"^ is absorbed by the earth’s atmosphere. As- 
sume that the radius of the sun and the distance between the sun and the earth are 
approximately 433,000 miles and 93,100,000 miles respectively. 

/XI-6. Calculate the direct radiant heat transfer between two parallel refractory 
surfaces 4 ft by 4 ft spaced 12 ft apart. Assume that the temperatures of the surfaces 
are 600 F and 1200 F respectively. The emissivity values for the low and high tem- 
peratures are 0.6 and 0.9 respectively. 

/ XI-7. Calculate the decrease in radiant heat loss per square foot of surface area 
through the wall of a house if two sheets of building paper are replaced by two sheets 
of aluminum foil. The surface temperatures of the air layer between the aluminum 
foils remain the same as they were with the paper, and are equal to 80 F and OF 
respectively, and the emissivity values for building paper and aluminum foil are 0.9 
and 0.08 respectively. 

XI-8. The value of the area under a 27 C curve in a figure similar to Fig. XI-1 was 
found equal to 0.047 watts cm"^. Calculate the Stefan-Boltzmann constant in terms 
, of B, hr, ft, F abs (= Rankine) units. 

XI-9. The walls of a furnace are made up of refractory lining (e = 0.8) separated 
‘ by a 4-in. air space from an outer casing of firebrick (e = 0.6), Calculate the radiant 
interchange across the air space if the lining and refractory temperatures are 2200 F 
and 1200 F respectively. 

,• XI-10. An oxidized pipe through which hot water flows at 200 F passes through 
a room at 70 F. Calculate the per cent reduction in radiant heat transmission based 
on the oxidized pipe loss, if the pipe is painted with aluminum paint. .Assume that 
the emissivity values for the oxidized and painted surfaces are 0.79 and 0.30 respec- 
tively. 
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CHAPTER XII 


HEAT TRANSFER BY THE COMBINED EFFECT OF 
CONDUCTION, CONVECTION, AND RADIATION 


XII-1 Heat Conduction in Series with Convection and Superimposed 
Radiation 

The combination of heat transfer by conduction and convection has 
been dealt with in Chapter IX. The combination of these two modes 
of heat transfer together with radiation will now be considered. The 

analysis of such complex prob- 
sudace 2 lems as the heat transfer within 

a furnace of a steam genera- 
ting unit, where conduction, 
convection, and surface and 
gas radiation are all involved, 
would indeed be very di ffi- 
cult. However, there are a 
few simple methods by which 
the influence of radiation to- 
gether with conduction and 
convection may be studied. 

Consider an insulated pipe 
of length L and cross-sectional 
area as shown in Fig. XII~1. 
It will be assumed that the 
pipe is located in a room the average* air and wall temperatures of 
which are equal to Heat is conducted radially from surface 1 to 2 
according to the equation 

_ k2,L{t, -t,) 

In rs/ri 



Fig. XII~1. Insulated pipe. 


The same amount of heat passing through the outer surface of the 
insulation is given up to the environment by natural convection and 
radiation. Fortunately radiation can be superimposed on the con- 
ductive and convective heat flow through a gas owing to its free passage 
through most gases. This means that the two kinds of heat transfer 
may be simply added as though they were like parallel electric currents. 

By referring to Fig. XII-1, the rate of heat flow from the surface 2 
to the air and the walls of the room is 


q = hMh - ts) + FsFA<rA 2 IXn-2] 
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where he is the coefficient of heat transfer due to convection. The 
symbols in the last term on the right-hand side of the equation are the 
same as those used in Eq. XI-19. From Eqs. XII-1 and 2, equating 
and substituting the numerical value of a, the following is obtained. 

_ fc2:rL(ti - h) 

® ” In ra/ri 

- - « + 0.174 - [ 0 ] lXII-31 

XII--2 Combined Coefficients of Convection and Radiation 

It is convenient at times to replace the usual radiation equation 
(Eq. XI~19) by an equivalent relation of the following form: 

qr = M2(fe - fe) [XII-4] 

In this equation hr is defined as a coefficient of radiation. Com- 
parison of Eqs. XI--19 and XII-4 shows that 

hr = FeFa<t{Tx^ + (Ti + T^) [XII-5] 

Substituting qr from Eq. XII-4 in Eq. XII-3 : 

I? = = hcAiitz — h) + hrAzih — t^) 

Simplifying 

9 = ^ = {he + hr)A^{h - h) [XII-6] 

In this relation the quantity Qic + K) represents a combination of the 
convection and radiation coefficients. 

Combined coefficients may be taken from Table XII-1. 

TABLE XII-1 

Heat Teansfer from Horizontal Tubes to Still Air at Orbinary 
Room Temperature 

Combined coefficients (he + W ia B hr”^ F~^ 

Outer Temperature Difference At[F] 

Diameter I I ! I I 1 

D, Inches 50 | 100 1 150 i 200 1 250 1 300 1 350 1 400 1 450 1 500 1 550 1 600 | 650 j 700 

1.3 2.26 2.502.73 3.003.29 3.60 3.95|4.34|4.73|5.16|5.60|6.05i6.5li6.98 

3.5 2.05 2.25 2.47 2.73 3.00 3.313.6914.0314.4314.8515.2615.716.1916.66 

5.6 1.95 2.15 2.36 2.612.903.20 3.5413.9014.30 

10.75 1.872.072.29 2.54|2.82 3.12 3.47|3.84| 
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By means of this table McAdams (Ref..XII-l) has represented experi- 
mental data of R. H. Heilman and L. B. McMillan for horizontal steel 
pipes having oxidized surfaces. Since the emissivity of oxidized steel is 
high and hr represents only a part of the combined coefficient, the values 
taken from the table may likewise be used for non-metallic surfaces 

Xn~3 Heat Losses from Bare or Insulated Horizontal Tubes 

In order to apply the equations used in Chapter III it was essential 
to know the surface temperature of the insulation. Often the tempera- 
ture of a tube is known; for example, it may be assumed equal to that 
of the fluid inside the tube, but the temperature at the surface of the 
insulation is unknown. For instance, if the most suitable thickness of 
an insulation is to be calculated, the inside temperature of the tube and 
the air temperature of the room are the only known temperatures. 
Equation XII-6 and Table XII-1 may be used in the solution of prob- 
lems of this type. 

Example XII-1. Calculate the heat loss per linear foot from a 4-in. (out 
side diameter = 4.5 in.) nominal horizontal steel pipe covered with 1 in. 0 
insulation {k = 0.035 B hr"^ ft”"^ F"^), if the pipe and still air temperature 
of the roam jtre 400 F and 70 F respectively. 

Solution: The actual outer, diameter of the tube is 4.5 in. and the oute 
diameter of the insulation is 6.5 in. For a pipe section 1 ft in length, E( 
XII-6 converts to 

Solving for the combined coefficient 

h, + hr 

r 2 {k — ^ 3 ) In r^/ri 

Substituting the given values: 

0352 ( 400 ^ 

(li-IO) 

Assuming the unknown surface temperature, h, equal to 110 F and s 
stituting gives, 

Ao + - 2.55 B hr-4r2 F"i 

For At = 110 — 70 = 40 F and D = 6.5 in., according to Table XT 
he hr == 1.9. The disargeement indicates that the assumed temperai 
was too low. 

Assume next that = 120 F. Then, 

Ac + A, = 1.97Bhr-ifr2F-i 
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For At = 120 — 70 — oO, Table XII—1 yields he -i- hr = 1.9, Tiius the as- 
sumption that U = 120 F is satisfactory, and from Eq. XII-7 
3 

1.9(27r) ~ (120 - 70) = 162 B hr-^ 

12 

In the case that large temperature differences exist in the insulation, 
it is sometimes necessaiy^ to select new thermal conducti\dty values as 
new surface temperatures are assumed. This, however, is only a matter 
of routine if charts are available which indicate the thermal conductivity 
of the material for various temperatures. 

Example XII-2. Calculate the heat loss from an uninsulated 2-in. nominal 
diameter (outside diameter = 2.375 in.) horizontal pipe at 400 F to the still 
air of a room at 70 F. 

Solution: For Af = 400 — 70 = 330 F, Table XII-1 gives he -i- hr = 3.67. 
Therefore, the heat loss ner ^quare foot of pipe is 

=1^(400 - 70) = 1210 B hr-^ ft-2 

In the design of insulation the teim insulation efficiency is often 
used. This is defined as the ratio of the heat saved by the insulation to 
the heat dissipated by the bare pipe. 

Example XII-S. The heat loss from a bare pipe (1000 B hr“^ ff^) was re- 
duced to 200 B hr“^ ft”^ by adding insulation. Calculate the efficiency of the 
insulation. 

Solution: The bare pipe heat loss is 1000 B hr~^ ft“b The rate of heat saved 
by the insulation is 1000 — 200 = 800 B hr~^ ft“^. Therefore, the so-called 
insulation efficiency is 800/1000 == 0.80 or 80 per cent. 

XII-^ Heat Transfer Through Air Spaces 

The combined action of thermal conduction, convection, and radi- 
ation is particularly interesting and important in air spaces. For 
instance, it must be considered thoroughly in the design of small frame 
buildings. The ideal case would exist if the convection and radiation 
were negligible compared with thermal conduction through the air space. 
So, the comparison of the effect of convection and radiation with that 
- of conduction is of considerable practical importance. 

Since free convection is governed by the laws of friction, it is obvious 
that the convection will be smaller in narrow spaces than in wide ones. 
It has been found that the convection in vertical air layers may be 
approximately determined by multiplying the value of the heat transfer 
by pure conduction with an augmentation factor of about 

1 J for |-in. thickness of the layer, 
for 1-in. thickness of the layer, 

2| for 2-in. thickness of the layer, and 

7| for 5-in. thickness of the layer, 
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provided that the temperature difference across the layer is 10 F. The 
factors become larger 'with, increasing temperature. From this it is 
seen how important it is for the heat instdation used in the walls and 
roofs of a building to have small air spaces. There is a lower limit, 
however, in the design, since too many narrow air layers require much 
more building material which has a higher thermal conductivity. 

In a frame wall of given total thickness, fewer and wider air spaces 
mean larger temperature difference for a single layer increasing the 
unfavorable effect. At the same time, heat transfer by radiation like- 
wise increases. At an average temperature of 80 F for instance and an 
air space 1 in. wide, the augmentation factor would be 2 instead of 1}^ 
for convection alone; and for a 2-in. wide space 3K instead of 2 %. 

The influence of radiation can be reduced by using surfaces of small 
emissivity. Aluminum foil has been extensively used for this purpose. 
Even by lining only one of the radiating surfaces with aluminum foil, 
the influence of radiation will be greatly minimized. The effectiveness 
of the reflecting surfaces, howm’^er, is largely dependent upon the perma- 
nency of the surface to retain the reflective characteristic. It is claimed 
that the enoissivity of aluminum is changed only slightly by the action 
of air. 

Another means of decreasing the convection heat transfer is by 
reducing the pressure. This is done in the so-called Dewar vessels 
(thermos bottles), which are evacuated double-walled glass vessels. 
The radiation heat transfer is reduced to a very small amount by silver- 
ing the surfaces of the evacuated space. 

PROBLEMS 

XII-1. Determine the heat loss from a 6-m. nominal diameter (outside diameter 
= 6.625 in.) horizontal steam pipe 160 ft long covered with of insulation 

{k = 0.04 B hr“^ ft“^ F“^), if the pipe and stiU air temperatures are 500 F and 80 F 
respectively. 

XII— 2. Compute the heat loss per square foot of outer surface from a 4-in. nominal 
diameter (outside diameter = 4.5 in.) horizontal steel pipe at 400 F covered with 3^ 
in. of insulation (k = 0.03 B ‘hr*"^ ft'~^ F”^) to the still air of a room at 70 F. 

XII— 3. Calculate the heat , loss -per linear foot from a 4-in. nominal diameter 
(outside diameter = 4.5 in.) horizontal pipe covered with in- of insulation A 

(kA = 0.056 B hr“^ ft“^ F”b which in turn is covered with 1 in. of insulation B 
(Jcb ” 0.039 hr~^ ft~^ F~^). The pipe surface and still air temperatures are 800 F 
and 80 F respectively. 

XII-4. Calculate the heat loss from 200 ft of 3-in. nominal diameter (outside 
diameter = 3.5 in.) horizontal steel pipe at 300 F to the still air at 80 F. 

XII-5. Determine the efficiency of insulation for a 6-in. nominal diameter (out- 
side diameter ~ 6,625 in.) horizontal steel pipe covered with in. of insulation 
(k — 0.04 B hr“^ ft”^ F””^). The pipe and still air temperatures are 600 F and 80 F 
respectively. 
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XII~6. Compute the steam condensed per hour per linear foot of bare horizontal 
6-in. nominal diameter (outside diameter = 6.625 in.) pipe, if the steam and surround- 
ing air temperatures are 275 F and 75 F respectively. Assume that the metal tem- 
perature is the same as the steam temperature. 

XII~7. Calculate the eifficiency of insulation for an 8-in. nominal diameter (out- 
side diameter = 8.625 in.) steam pipe covered with IH in. of insulation A next to the 
pipe (kA = 0.052 B hr~^ ft"^ F“^) and 2}4 in. of outer insulation B {ks — 0.041 
B hr""^ ft~^ F"^) if the pipe and still air temperatures are 850 F and 70 F respectively. 

XII"8. Calculate a combined natural convection and radiation coefficient for a 
12-in. nominal diameter (outside diameter = 12.75 in.) horizontal pipe at 300 F in 
contact with still air at 70 F by the laws of convection and radiation. Compare the 
calculated value with the one obtained by use of Table XII-1. 

XII-9. Calculate the natural convection and the radiant heat losses in per cent 
of the total loss from a 4-in. nominal diameter (outside diameter = 4.5 in.) horizontal 
pipe at 200 F to the still air of a room at 70 F. 

XII-10. Rework Problem XII~9 using a pipe temperature of 500 F. Compare the 
results for the two problems. 
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CHAPTER XIII 


EXPERIMENTAL DETERMINATION OF CONDUCTIVITIES 
AND EMISSMTIES 

Xni-l specific Properties of Matter 

In dealing with the laws of heat transfer and insulation and their 
applications, various physical properties were used such as density, 
viscosity, specific heat, thermal conductivity, and emissivity. Whereas 
the first three occur in different branches of science, the last two are 
specific for heat transfer and insulation alone. There are numerous 
experimental methods by which these properties can be determined, 
either as to their absolute values, or relative to known values of any 
standard materials. To give an idea of the procedure involved in such 
tests, some of the simplest and more common methods will be briefly 
described. 

Xni-2 Measurement of the Thermal Conductivity of Metals 

For a determination of the thermal conductivity of a metal relative 
to another, two cylindrical bars of equal cross-sectional area A are 
soldered together end to end. One end of the composite specimen is 
attached to a heater and the other end to a cooling element. The cylin- 
drical surfaces are insulated against heat losses. Then heat energy at 
the same rate q flows through either of the rods. The temperature 
drop (A<)i over a length Li of rod 1, and (A<)2 for a length L2 of rod 2 
are measured by means of thermocouples. According to Eq. II-l 

7 A 7 A rvTTT n 

q = kiA — — = hA — — [XIII-l] 

Li 1/2 

where fcj and h2 are the thermal conductivities of the two materials, or 

(AO2 Li 

No measurement of 5; is needed with this method. If the thermal 
conductivity ki of the metal of one rod is known, then the conductivity 
^2 for the other rod may be determined from the temperature measure- 
ments and Eq. XIII-2. Practical arrangements are described for 
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instance in two research papers of the Bureau of Standards (Ref. XIII-l 
and 2 ) where this method has been used for determining the thermal 
conductivity of various metals up to 1100 F. 

XIII-3 Measurement of the Conductivity of Insulating and Building 
Materials 

The most common method for the determination of the conductivity 
of insulating and building materials which are available in large plates 
is the so-called twin-plate method, it consists of placing a flat elec- 
trical heating plate H between two equal plates Pi and P2 composed of 
the material to be tested. These three plates are then placed between 
two cooling plates Ci and C2 (see Fig. XIII- 1 ). The heat energy pro- 
duced in H splits into two equal parts. Half of the heat flows vertically 
upward through plate Pi and half downw^ard through plate P2. The 
cooling plates Ci and C2 usually are hollow bodies in which cooling 
water or other liquids pass to and fro through the sections of the cooler. ' 
Temperatures are measured by thermocouples arranged at the inter- 
faces of H and Pi, H and P2, Pi and Ci, P2 and C2. If the plates Pi 
and P2 are equal in thickness L, area A and thermal conductivity, and 
if the rate of heat energy q of the heater H is known by means of elec- 
trical measurements, then according to Eq. II-l 



where M is the temperature drop across either of the plates. If there 
are small differences, an average of At may be taken. 

From the edges of plates Pi and P2 heat energy would be lost if they 
were in contact with the surrounding air. In order to avoid such losses, 
guard rings are used. In Fig. XIII -1 a heater H' similar in construc- 
tion to heater P, surrounds the latter. Plates Pi and P2 are sur- 
rounded by plates P/ and P2' which are constructed from the same 
material. If by appropriate regulation of the electric heating of 
the same temperature drop across Pi and P^^ exists as across Pi and 
P2, then no heat energy from H to P', or from Pi to Pi' or from 
P2 to P2' since no temperature differences exist which might drive heat 
across the vertical gap. On the outer edges of P', Pi' and P2 heat 
will be lost which is produced by P'. In order to reduce this as much 
as possible, the whole arrangement may be surrounded by a loose 
insulating material, such as diatomaceous earth. 

Instead of arranging special guard rings, a single heater as large as 
P and P' and test plates as large as Pi and Pi may be used. In this 
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case, however, particular care must be taken to avoid, errors; for 
instance, thermocouple readings close to the edge must not be used. 

Materials which are on the market in cylindrical shape, such as pipe 
insulation (see Figs. II-4, 7, 9, and 11) are generally tested in cylindrical 



Fig. XIII-1. Twin-plate arrangement for determining the thermal conductivity of 
insulating and building materials. 

Cl, C% Cooling plates Pi, P2 Plates of material to be tested 

H Electrical heating plate Pi, P2 Guard ring plates 

Electrical ring heater 

form. A uniformly wound coil with a tube as a core serves as the 
electric heater, and a second tube is placed over this. The second tube 
is covered with the material to be tested. Thermocouples are arranged 
I on the second tube and on the cylindrical outer surface of the material 
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to be tested. For this case, Eq. III-6 for hollow cylinders must be 
applied in order to determine the value of h. Heat losses at the two 
ends of the cylindrical sections can be avoided by guard coils on the 
ends which eliminate heat loss in the axial direction. They correspond 
to the ring heater in the twin-plate device. The end heaters' may be 
dispensed with if long tubes are used and the thermocouples are placed 
at a sufficient distance from the ends. 

XIII-4 Measurement of the Conductivity of Liquids and Gases 

The thermal conductivity of liquids and gases likewise can be deter- 
mined using either plane or cylindrical layers of the fluid. For instance, 
the thermal conductivity of a liquid may be determined by using an 
apparatus consisting of an upper heating disk and a low^er cooling plate 
separated by a space of J^^-in. or less w^hich is filled with the liquid. No 
heat transfer by convection will take place since the density of the fluid 
decreases from the lower to the upper surface. By measuring the heat 
input, the temperatures of the disks and the vertical distance of sepa- 
ration between the plates, the thermal conductivity of the liquid may 
be obtained by means of Eq. II-l. The upward and edge heat losses 
may be avoided By placing a wide Dewar vessel (thermos bottle, see 
Sect. XII-4) over the apparatus. By this method part of the curve in 
Fig. II-2 was determined (Ref. XIII--3). 

The co-axial cylinder method consists of arranging a fine platinum 
wire at the axis o^f a metal tube in w^hich is placed the liquid or gas 
under test. The apparatus is usually placed in a liquid bath in order 
to regulate the surrounding temperature. The space betW'Cen the 
inner surface of the tube and the surface of the wire forms a hollow 
cylinder. The platinum wire serves as an electric heater and as a 
resistance thermometer.* Then, referring to Eq. III-6, the rate of 
heat energy q' and the temperature h are measured electrically. The 
temperature ^2 is practically equal to the temperature of the liquid bath. 
From these and measurements of the wire diameter 2ri and the tube 
width 2 r 2 the thermal conductivity of the liquid or gas may be obtained 
from Eq. III-6. The influence of convection can be minimized by using 
a narrow tube, or in the case of gases, by operating at low pressure 
(see Sect. XII^). 

XIII-5 Measurement of Emissivities 

The accurate measurement of emissivities of materials is rather 
difficult. Only one of the simplest methods, used at room temperature 

* The principle of operation of the electric resistance thermometers is based on the 
well-defined relation between the change of the electric resistance of pure metals and 
the temperature. 
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(Ref. XIII--4), will be briefly described (see Fig. XIII~-2). The hollow 
cylinder a Adth double walls for circulation of cooling water (inlet at h 
outlet with a mercury thermometer and a thermocouple at c) contains 
a total radiation receiver d as shown in Fig. XIII-3. The hollow 
vertical copper cylinder 1 is gold plated on the outside and blackened 
on the inside. The two extremely thin copper disks 2 and 3 are black- 



Fig. XIII- 2 . Device for relative 
measurement of emissivities at ordi- 
nary temperature. 

a Double wall container 
h Cooling water inlet 

c Cooling water outlet with thermometer and 
thermocouple 
d Radiation receiver 
e Thermocouple 

/ Disk of material to be tested, with thermo- 
couple in center. 



Fig. XIII-3. Radiation receiver d of 
the device represented in Fig. XIII-' 2 . 

1 Hollow copper cylinder 
2, 3 Fine copper disks 
4 Bismuth wire 
5, 6 Wires of bismuth-tin alloy 


ened on the side upon which falls the radiation from the outside. These 
disks which receive the radiation through the equal angles (p are heated 
by the incoming radiation. The warming up of the disks is measured 
by a system of couples consisting of a fine bismuth wire 4 and two 
bismuth-tin alloy wires 5 and 6. Wires 4 and 5 and 4 and 6 constitute 
two equal thermocouples having a thermoelectric force of 61 microvolts 
per degree Fahrenheit, or three times as much as an equivalent copper- 
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constantan system of couples. The thermocouples are connected so as 
to oppose each other. In this way, the thermoelectric effects produced 
by the absorption of equal amounts of radiation received by 2 and 3 
cancel; for example, radiation from the copper cylinder 1 to the two 
thermocouples. Only the difference in radiation entering from above 
and below through the equal angles <p produces an effect. 

By referring to Fig. XIII-2, it is ob\dous that one of the two radi- 
ations comes from the disk /, the emissivity of which is desired, and 
the other comes from the inner blackened surface of the container a. 
The surface of container a is maintained at a constant temperature 
which is indicated by the thermocouple e. The radiation of / at a 
temperature which is measured by a thermocouple indicated in the 
figure, is compared with that from a black body at the same temperature 
by replacing / with a black body. If the galvanometer to which the 
thermocouple system d is connected is deflected by an angle 8^ with the 
black body at the opening and hj 8x with the disk/ at the same opening, 
then, aside from some corrections, the emissivity of the material of disk 
/ which was required is ex = B^/dx. 
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CHAPTER XIV 

HEAT TRANSFER IN TEMPERATURE MEASUREMENTS 

XI7-1 Heat Exchange, a Necessary Condition for Temperature 
Measurements 

It is more or less an accepted notion in engineering practice that 
temperature measurements can be made quite accurately. This has 
also been implied in the former chapters of this book. This so-caBed 
trust in temperature measurements is probably based on the high 
standards maintained by the thermometric industry and the avail- 
ability of very accurate industrial thermometric devices. Users of 
precision thermometric instruments often overlook the fact that accu- 
rate results are obtained only through the correct application and use 
of the devices, that is, the purchase, alone, of a good instrument does 
not assure correct temperature measurements. 

Thermometry is dependent upon heat exchange. When a mercury 
thermometer or a thermocouple is brought in contact with a solid or 
fluid the temperature of which is to be measured, the instrument is 
supposed to assume the temperature of the solid or fluid by conduction, 
convection, radiation, or a combination of two or more. A radiation 
pyrometer operates by means of the absorbed energy emitted by a sur- 
face toward which it is directed. It is a fact that conduction, con- 
vection, and radiation occur only when a temperature difference exists. 
Thus, theoretically, there is always a difference between the tempera- 
ture to be measured and that directly indicated by a thermometric 
device. It does not necessarily follow that this temperature difference . 
will be very small. 

XIV-2 Measurement of Surface Temperatures 

Surface temperatures are often measured by use of a thermocouple 
attached to the surface. If the lead wires are arranged so that they go 
directly from the surface in a perpendicular direction, the system is 
comparable to a bar cooled in air, as discussed in Sect. IX-2. Heat 
from the surface whose temperature is desired is conducted to the couple 
contact point, and is thence conducted along the wires and finally dis- 
sipated to the surrounding air. Since a temperature difference must 
exist when heat flows, that part of a body which is in direct contact 
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withNfche thermocouple, will be cooler than neighboring parts of the body, 
and, as a result, the indicated temperature will be too low. The 
error involved will be larger for bodies having small conductivity values, 
since theNheat conducted away by the lead wires will cause a greater 
temperatu^ drop in the body. Thermocouple experiments by H. 
Reiher usin^lates at 96 F in air at 59 F showed that a thermocouple 
read 73 F if thAplate used was cork, and 89 F if made of copper. When 
the joint of the\^ermocouple was arranged at the center of a thin 
copper disk and tms^^^disk in turn attached to the surface, the readings " 
were 90 F for a cork^te and 94 F for a copper plate. The reason for 
the improvement was\hat the heat instead of coming from a tiny 
part of the surface was ta^n from the rather large area under the disk. 
When the lead wires were 1^ on the surface of the disk, a distance of 
4 in., instead of being brou^ directly away from the surface the 
temperature indicated was 96 F\for both the cork and copper plates. 
From these results it is apparen\^^that large temperature measuring 
errors may result if thermocouples ^ not properly placed on the sur- 
face of insulating materials. As a 
general rule, surface thermocouple 
lead wires should always be arranged 
so that they lie close to the surface 
for a distance of a few inches in 
order to reduce the error due to 
conduction of heat along the wires. 

XIV-3 Influence of Conduction and 
Convection on the Measurement 
of the Temperature of Flowing 
Gases 

The important case of a ther- 
mometer immersed in a flowing fluid 
can be studied by means of the 
analysis given in Sect. IX-2. The 
upper part of Fig. XIV-1 shows a 
thermocouple inserted with the junc- 
tion at the bottom of the well W XIV— 1. Temperature distribution 

placed in a pipe through which a along a thermometer well, 

gas flows at constant velocity. 

The temperature of the gas is assumed higher than that of the pipe. 
The length of the well is denoted by L, the cross-sectional area of its 
wall by A, and the circumference by C. The temperature distribution 
of the gas tg across the pipe section is shown in the lower part of the 
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figure It may be assumed that tg is approximately constant without 
introducing a large error. The temperature of the well will increase 
from the pipe wall temperature h to the temperature at the end of 
the well. If the thermocouple wires are fine, heat conduction along 
them may be neglected. Radiation between the well and the pipe 

may likewise be negligible. ■ i j.i. 

It is a matter of simple calculation to determine whether ti, is con- 
siderably lower than the gas temperature tg which is to be measured by 
the thermocouple. If the well is so long that at a: = L no more heat 
will enter, then it is as though the bottom of the well were perfectly 
insulated.* This case is identical with the one dealt with in Ex. IX-1, 
and solved by the formula: 

= ^0 ^ g-mh ffX-13] 


According to the definition, 6 = it — tg and - <o is, both are 
negative. Thus, the error in the reading of the thermocouple in the 

well is 


{tg - to) ■ 


tg - to 

cosh (mL) 


[XIV-I] 


where m = VhC/kA, k = the thermal coiiductivity of the well mate- 
rial, and h can be calculated from Eq. VIII-5. Substituting the con- 
stants formerly given, the equation becomes 

(Nu) = 0.3 [XIV-2] 


where {Nu) = - — and (Re) — , and 

tOg 

D — the outer diameter of the well, 

V = the average velocity, | 
kg = the thermal conductivity, | of the gas. 
vg = the kinematic viscosity J 

Equation XIV-1 shows that the thermometric error increases with 
{t — to). The error becomes larger with increasing k and A and with 
de°creasing L, h, and C, since the hyperbolic cosine of (mL) increases 
with mL. The error may be reduced by any one of the following 

means; i • n • 

1. The application of insulation to the outside of the pipe, wau m 

order to bring the wall temperature to closer to the gas temperature tg. 
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As a rule, the temperatures of hot gases or vapors, for example, super- 
heated steam, should never be measured in uninsulated pipes. 

2. Make the well length L as great as possible. 

3. Make the cross-sectional area A of the wall as small as possible. 

4. Construct the well of material having a low thermal conductivity, k, 

5. Make h as large as possible. This may be accomplished chiefly 
by increasing the velocity. Thus, if a choice is possible the temperature 
of the fluid should be measured where the velocity is high. The so- 
called suction pyrometer (see Sect. XIV~4) is based on this idea. 
From Eqs. XIV-1 and 2 it may also be concluded that an increase in 
the diameter D of the well will improve the results slightly. 

XIV--4 Influence of Convection and Radiation in the Measurement of 
the Temperature of Flowing Gases 

In measuring the temperature of hot gases flowing through ducts the 
influence of radiation from the well to the cold walls must be carefully 
considered. In tests performed at Purdue University (Ref. XIV-1) in 
order to determine the errors to be expected in the measurement of 
flue-gas temperatures, a bare thermocouple indicated a temperature of 
about 775 F instead of the true gas temperature of 850 F. By placing 
the couple in a ?^-in. pipe and drawing a small part of the gas through 
the pipe by means of a steam ejector (suction pyrometer), this error was 
reduced considerably. 

Another method used for reducing the error is to place the thermo- 
couple in a protecting metal shield which does not hinder the gas flow, 
but shields the thermocouple from radiation to the duct w^all. The 
shield receives energy from the thermocouple by radiation and from the 
flowing gas by convection, and it radiates these energy amounts to the 
duct wall. This requires an appreciable temperature difference, and 
therefore the shield temperature is not much below that of the gas. 
Thus the thermocouple must not give up much energy by radiation and 
its temperature remains close to that of the gas. This means the 
error in the temperature reading will be small. 

Consider what happens when the duct wall is not insulated and the 
other precautionary measures are not taken. Consider an element of 
the well W in Fig. XIV-1 close to the inner end where t = constant = ir- 
The radiation of this element to the wall at temperature will be con- 
sidered. Since the temperature of the well element is assumed almost 
constant, heat flow by conduction from the inner end through the wall 
of well W in the direction toward the pipe wall is negligible. Therefore 
the heat delivered by convection from the gas to the well will be given 
up again only by radiation to the wall of the duct and, according to 
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Eq. XII-2, 

PCIV-3] 

In this relation A is the surface of the element under consideration at 
the inner end of the well. The A terms drop out, since one appears on 
each side of the relation. If the duct is wide, then Fa ^ 1. Fe is 
the emissivity factor of the well surface. The coefficient of heat trans- 
fer h may be found from Eq. XIV-2. 

Example XIV-L In a 36-m. diameter circular duct a %-in. diameter 
tube is placed at right angles to the direction of the gas flow. A thermocouple 
housed in the tube indicates a temperature of 800 F. If the duct wall tem- 
perature equals 500 F, the average gas velocity 15 ft/sec, and the specific weight 
of the gas is 0.04 Ib/ft^, calculate the approximate true gas temperature. As- 
sume that the viscosity and thermal conductivity are 0.18(10“^) lb ft“- hr and 
0.02 B hr~^ ft"^ F“^ respectively. The emissivity of the oxidized surface of 
the thermocouple housing is assumed equal to 0.8. 

Solution: 

Tl = 460 -b 800 = 1260 R 
To = 460 -f 500 = 960 R 

fx - 0.18(10-9)36002 slug ft-i hr-' (See Sect. VI-2) 
slug ft-® (See Sect. 1-3) 

32.2 

j- = - = 0.18(10-^)36002 — ft® hr-i 
0.04 

Substituting in Eq 

0.625 
0 . 02 ( 12 ) 

h - 7,44 B hr-i fr^ F-'- 
Substituting this in Eq. XIV-3 

7M(tg - 800) 0.174[12.64 - 9.64]0.8 

and 

tg ■■ 1113 F 

This means that the thermocouple would indicate a temperature approxi- 
mately 300 F too low. 
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XIV-5 Temperature Measurement with Radiation Pyrometers 

Whereas in ordinary thermometiy^ practice radiation is only of interest 
as a source of errors, high-temperature measurement is based almost 
entirely upon radiation. Only the main principles used in radiation 
p^Tometry will be mentioned here. 

The two types of radiation pyrometers are the optical pyrometer 
and the total radiation pyrometer. 

Optical or selective radiation pyrometers are based on the change of 
intensity of radiation vith temperature at a given wavelength. For 
black radiators this is knovm exac% from Planck's law (see Fig. XI-1). 
In this type of instrument the brightness of a glowing thread is com- 
pared wdth the brightness of a surface the temperature of which is to 
be measured. The brightness of the thread is controlled by an electric 
current ’which heats the thread until its brightness is equal to that of 
the surface under examination. Obser'vdng the thread by a telescope 
directed at the surface under consideration enables the obsertmr to 
adjust the current so that the brightness of the two are equal, at which 
point the thread becomes invisible. A red filter is placed at the ocular 
of the pyrometer in order to produce monochromatic light. The device 
is calibrated by comparison with a black body. If the radiating surface 
the temperature of which is desired is not perfectly black, the true 
temperature may be found by use of the emissmty together with the 
lawo of radiation. 

A total radiation device has been discussed previously in the section 
dealing with the experimental determination of emissivities (Sect. 
XIII-5). The total radiation pyrometer is based on the Stefan- 
Boltzmann law and it consists of an optical system w^hich is designed to 
collect radiant energy of all wavelengths and focus the energy on a heat- 
sensitive element such as a thermocouple. The thermoelectric current 
originating from the heating up of the one junction of the thermocouple 
is a measure of the total radiation which in turn is a measure of the 
temperature of the radiating surface. Because this type of pjuometer 
measures the temperature objectively, and not subjectively as the op- 
tical pyrometer, it is possible to indicate and record the reading of 
the instrument at any distance from the observed surface. For non- 
black body radiators the total radiation pyrometer is unsuited. It is 
most frequently used for measuring the temperature of furnaces, which 
may be considered as approximately black radiators. 

For a complete discussion of errors, suitability of various instruments 
for different applications, and other information relative to temperature 
measurements, the reader is referred to the Power Test Codes of the 
American Society of Mechanical Engineers (Ref. XIV~2). 
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CHAPTER XV 

HEAT TRANSFER AND PRESSURE DROP 


XV-1 Practical Importance of the Relation Between Heat Transfer 
and Fluid Friction 

There are various relations between the laws of heat transmission 
and other physical laws. For instance, the ratio of thermal and electri- 
cal conductivity at a given temperature has the same value for all pure 
metals. Another example is the close analogy of mass exchange by 
diffusion of vapors and gases to the convective heat transfer. 

In this chapter one of these analogies, which is most important for 
mechanical and chemical engineering, will be dealt with briefly. It con- 
cerns Osborne Reynolds’ suggestion from the year 1874 (Ref. XY-1) 
that momentum and heat in a fluid are transferred in the same way so 
that in geometrically similar systems, fluid friction and heat transfer 
should be proportional. This analogy is of particular interest for 
engineers because it permits the predicting of heat transfer coefficients 
from friction data which can be obtained with less difficulty from experi- 
ments than heat transfer data. It further sets the conditions of an 
economical optimum in design for the following reasons. 

As has been shown m Chapter VIII, heat transfer generally can be 
improved by increasing the speed of the fluids involved. On the other 
hand the pressure drop of fluids which stream in conducts increases with 
velocity. This means that more power is required for pumping the 
fluids through a heat exchanger. Therefore, the investment costs for 
blowers or pumps and the current costs of pumping power become 
greater, and eventually may offset the economies resulting from the 
reduction of the heating surface which is possible if higher velocities are 
used. Therefore, a compromise has to be found between the require- 
ment of high velocity from the heat transfer view and the necessity of 
not too hi gh velocity in order to keep the costs of overcoming the pressure 
drop within reasonable limits. 

Some excellent presentations and modifications of Reynolds’ analogy 
have been published recently which, however, use more mathematics 
and fluid mechanics than can be required in this text. Therefore, a 
simpler, though theoretically not so strongly founded way will be 
followed which has been indicated by White (Ref. XV-2). 
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XV-2 Basic Idea of Reynolds’ Analogy 

Reynolds' reasoning was, that from a molecular viewpoint fluid friction 
and heat conduction obey the same law of motion. If two neighboring 
layers in a gas moving along a plane wall have different velocity, then 
in the course of molecular motion, molecules of the faster layer pass to 
the slower one and are retarded there, whereas molecules of the slower 
layer come into the faster layer and are accelerated. By this inter- 
change of molecules, the faster layer will be retarded as a whole, and the 
slower one will assume a greater speed. Owing to the mass of the mole- 
cules an exchange of momentum occurs which acts like friction and can 
be expressed by Newton’s equation: 

„ . dv 

= [VI-2] 


where y is the distance from the wall of the layer under consideration. 

This equation was used for the definition of the dynamic viscosity n. 
According to the foregoing remarks may be considered as the mass per 
unit time and area which is interchanged by molecular action between 
layers unit distance apart. 

This can be shown more directly from Eq. VI-2 by comparing the 
latter with Newton's basic law of dynamics which may be expressed in 
the form 

T, dv 

Fs = -m- [XV-1] 


where m is the mass which by the shearing or frictional force Fs is decel- 
erated, and T is the time. 

From Eqs. VI-2 and XV-1 : 


dv , dv 
—m — = juA — 
dr dy 


[XV-2] 


partials being used because now two independent variables (y and r) 
are involved. But according to a mathematical equation about the 
partial differentials of an analytical function z = f{x,y): 


dz dx dy 
dx dy dz 


[XV-3] 


In our case v is a function of y and t. Therefore, corresponding to 
Eq. XV-3 

dv dr dy 
dr dy dv 


[XV^] 
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and by substitution in Eq. XV-2 

m— = M 

or 

or 

m dy 


[XV-5] 


For unit length and time interval, dy/dr = 1, and unit area means 
A — 1. Then /X = m, that is, the viscosity equals the decelerated mass 
as has been asserted. 

If not only velocity differences, but also temperature differences exist 
between the two layers under consideration, the molecules vill not only 
exchange momentum, but vdll likewise convey heat energ}^ from the 
warmer to the colder layer. The heat transportation w^hich coitc- 
sponds to the product p. • dv in Eq. VI-2 obviously should be proportional 
to the interchanged mass ju, the specific heat Cp of the fluid and to the 
temperature difference dt Thus one arrives at 

df 

Q = /'CpA — p;V-6] 


as an equation corresponding to Eq. YI-2. 

However, according to the basic equation of heat conduction 

-kA f pCV-7] 

dy 

The difference in the signs in Eqs. XV-6 and 7 does not matter because 
it is due to a difference in definition only. Equation XV-6 is derived 
from Eq. VI-2 where > 0 means transportation of momentum in the 
direction to the wall, whereas in Eq. XV-7 g > 0 means transportation 
of heat in the direction from the w’^all. 

Thus, comparison of Eqs. XV-6 and 7 shows that Eejmolds' analog}^ 
holds only if pCp = k, or 



The dimensionless group pCp/k occurs in all kinds of convection. It 
has been called Prandtl number and denoted by the symbol (Pr) in 
Sect. VI-7. 

As mentioned in Sect. VIII-2, (Pr) does not change much for gases, 
and with an average value of 0.8, is not far from the value 1. Because 
further (Pr)” usually appears in the equations of heat transfer where 
71 < 1, Reynolds’ analogy will hold for gases with good approximation. 
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It must be kept in mind, however, that Sow in layers, i.e., streamline 
Sow was assumed in the approximate derivation given above. 

XV-3 Generalization of Reynolds’ Concept 

Reynolds' analogy of heat transfer and Suid friction in its original 
form is not only restricted to gases, more exactly to a fluid with (Pr) = 
but also to streamline flow. In order to get rid of the latter restriction 
a fictitious or apparent viscosity fia will be introduced which includes the 
exchange of momentum by mixing of larger aggregates of molecules as 
occurs in turbulence, instead of the merely molecular mixing. In 
making use of fia, one must keep well in mind that this is a magnitude 
which depends upon the motion and so may vary from point to point in 
the stream whereas the molecular or true viscosity ^ is a property of the 
fluid which in a rather wide pressure range depends on the temperature 
only and nothing else. So /Xa may be almost equal to ja at one point 
of the fluid cross section, and 100 times larger than ja at another 
point. 

In a similar way, instead of the thermal conductivity fc, a fictitious or 
apparent conductivity ha may be introduced, which takes care of the 
exchange of heat energy between the mentioned larger aggregates of 
molecules. This magnitude is likewise very different from the true 
thermal conductivity, nor can it be assumed that Ma/^a = ix/h or 
y^aOplha = MCp/fc. Considering for instance an aggregate of molecules 
which moves from one layer to another one and then returns to the 
first one, momentum may be exchanged so rapidly that the exchange 
is completed when the aggregate returns to the original layer, but the 
interchange of heat energy may not be completed in this time interval. 
Then ixa = ix, but ha < h, and by this \iaOplha > iiCpIk. For this and 
other reasons it is not at all certain a 'priori that Reynolds’ analogy can 
be used in turbulent flow. 

A comparison between the experimentally observed distributions of 
temperature and velocity has been considered as a simple test of the 
theory. If the errors due to the made approximations are unimportant, 
then the two distributions in any plane normal to the flow direction 
should be similar. 

The first to perform such experiments was Pannell (Ref. XV-S). 
He measured the velocity and temperature distribution across air which 
streamed in a vertical brass tube, 1.92-in. inside diameter. The velocity 
and temperature in the center line of the tube were Vc = 87.4 ft/sec, 
and tc = 75.9 F respectively, the wall temperature was 4 = 109.4 F. 
The test cross section was far removed from the pipe entrance so that 
the effect of entrance disturbances was greatly diminished. 
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The air velocity and temperature at a radial distance y from the tube 
wall may be called v and t respectively. Introducing further d ^ tw — 
the velocity ratio v/vc can be compared wdth the temperature ratio 6 /6c 
for different values of the distance ratio y/r, where r is the distance of 
the center line from the wall, and Be — tw 



Fig. XV~1. Distribution of velocity and temperature across an air stream in a 
heated tube, according to measurements of J. R. Pannell, in dimensionless rep- 
resentation. 

In this way PannelFs main results are represented in Fig. XV-1 and 
Table XV~1. Similarity ob^dously exists if v/vc = 6/6c for any radial 
distance, and the ratios shown in the figure and the table indeed corre- 
spond to this condition with good approximation. 

TABLE XV-1 

Veloctty and Temperatube Ratio for Turbulent Flow of Air 
IN A Tube, According to Measurements of J. R. Pannell 


y 

r 

V 

Vc 

e 

da 

0 

0 

0 

0.05 

0.66 

0.65 

0.1 

0.73 

0.73 

0.2 

0.80 

0.81 

0.4 

0.89 

0.92 

0.6 

0.96 

0.98 

0.8 

0.99 

1.00 

1.0 

1.00 

1.00 


White (Ref. XV-2) has pointed to the strange fact that this agreement 
is much better than for streamline motion for which Reynolds' analogy 
was derived originally. From this it is seen that Reynolds' theory does 
not hold exactly. The main reasons for the deviations seem to be a 
certain lack in the assumed similarity of heat transfer and interchange 
of momentum in a tube. This, however, cannot be dealt with here in 
detail. 
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XV-4 Equations Based on Re 5 molds’ Analogy 

In order to calculate the heat flow from the wall to the fluid it is 
necessary to know the temperature slope dd/dy in the fluid adjacent to 
the wall. This is very steep for turbulent 
flow as can be seen from Fig. XV-1 where 
the curves begin almost vertically at the left- 
hand side. Fortunately, the slope of the curves 
can be expressed in terms of heat flow and fluid 
friction. 

Figure XV-2 is a schematic representation 
of similar velocity and temperatui^e distribu- 
tions across the fluid stream in a tube. The 
similarity of the two curves can be checked 
by comparing the angles between the horizon- 
tal axis and straight lines from the zero point 
to points of the two curves at the same dis- 
tance y. Let these angles be a and ^ respec- 
tively. It is seen that they increase with 
decreasing v and 6 up to maximum values of 
ao and /So, determined by the tangents of the 
curves at the zero point of the system. 

Obviously, a condition of similarity is, 
that for any distance y between y — 0 and 
y 

tan a tan /3 
tan ao tan /Sq 



Fig. XV~2. Schematic 
representation of velocity 
and temperature distribu- 
tion across a fluid stream 
in a heated tube. (Since 
e = tw — tj the tempera- 
ture t is great where d is 
small, and vice versa.) 


[XV-9] 


The four tangent terms may be expressed by the pertinent geometrical 
and physical data. 

From Fig. XV-2 it is seen that 

tan a = - . pLV-10] 

y 

tan ^ = - [XV-ll] 

y 

According to Eq. VI-2, the frictional force on the wall is 
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According to Eq. XV-7, the rate of heat flow from the wall to the 
fluid is 

q^=-kA ® 

\dy/y^0 

Substituting t ^ - 6, according to the definition of and con- 

sidering that the wall temperature ty, is supposed to be constant, it 
follows that 

qw = kA Q [XY-IS\ 


By substitution from Eqs. XV-10, 11, 12, and 13 in Eq. XY-9 
one obtains 

QwV ^ 
dk fxv 
or 

q%oV k kC'D Cp 

Fyfi IJL fiCp (Pr) 


or 


Qw^iPr) 

Fy^dCp 


[XV-14] 


This, however, has been derived under the assumption of laminar 
flow. For turbulent motion it w^ould hold only if jjLaCpfka - iJiCp/k 
'which is unlikely. Generally, txaCpIka vill be greater than ixCp/k, as 
has been explained in Sect. XV-3. Furthermore, the lack in similarity, 
mentioned at the end of that section, must be considered. 

For these reasons YTiite recommended the use of the following 
modified equation: 


Qw'^miP'F) 

F vfirrfip 


|XV~15] 


where and are mean values, tyn, tw ~ being the mbdng-cup 
temperature as defbied in Sect. VIII-1. 

The exact theory for streamline flow and uniform heating from the 
wall shows that 

<i> = 0.545 [XV-16]* 


lu ^77t* 

* According to Eqs. XV-18 and 19, Fy, - Ap — . E the tube is unifonnly 

(ice) 2 

heated from outside, qw = 2.18(/:/r) as derived first by Callendar (Ref. XV-4). 
Substituting these two expressions in Eq. XV-15, one obtains Eq. XV--16. 
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For turbulent flow, White found the following empirical relation: 




750 VC^ 7.5 
{Re) 


[XV-17] 


This gives rather satisfactory results for Re ^ 6000 and Pr 5 50. 
If the application of Eq. XV-16 is restricted to quite moderate temper- 
ature differences, dm, the calculation of (Re) and (Pr) may be based 
on the average temperature of the flowing fluid. The more compli- 
cated procedure for larger values of dm will not be dealt with in this book. 

The effect of disturbances at the entrance of the tube is not included 
in Eq. XV-17; it will be considered approximately in Sect. XV-5. 

In fluid dynamics it is usual to express the frictional force P™, which 
is exerted on the inner surface of a tube, in terms of the kinetic energy 
per unit volume. This is done by introduction of a dimensionless 
coefficient of friction, often called Fanning’s factor: 


__ FJA FJA 

Vm^ 7 ^ 

^2 g 2 


CXV-18] 


A, in this equation, is the area of the inner surface; p and y are the 
density and specific weight, respectively, as defined in Sect. 1-3. 

The exact theory of streamline flow leads to 

f-W) 

For turbulent flow the empirical equation 

0.08 

/ = CSV-20] 


holds at least up to {Re) = 200,000. 

Solving now Eq. XV-15 for and substituting Fy, from Eq. XV-18 
one obtains 

2g{^r) ^ / 2{Pr) P21V-21]* 

An average film coefficient of heat transfer hm is defined in the usual 
manner by the formula 

= hrnAOm \X.Y-22] 


y 

* According to Sect. IV-1: Cp = pCp = — Cp. 
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From Eqs. XV -21 and 22 : 

hm = 4>f 


7 Cp 

2g{Pry 


4>S 


2(Pr) 


lX\-23] 


The right side of this equation contains only the diameter of the 
tube, the velocity, and some properties of the flowing substance. Thus, 
the coefficient of heat transfer is determined without reference to heat 
flow or temperature. 

In engineering, it is more usual to operate with the pressure drop 
Ap in a pipe than vdth the frictional force exerted on the inner 
surface. Ap may be introduced by means of the following relation 
which is used generally in fluid d^mamics: 


L 2 T 2 


D 2 


D 2< 


i:XV-24] 


L and D in this equation are the length and inner diameter of the pipe, 
respectively. It is seen that expressing lengths in feet and 7 in Ib/cu ft, 
Eq. XV~-24 yields Ap in Ib/sq ft. 

Substitution of / from this equation in Eq. XV--23 leads to 


f. ^ ^ 

4 (Pr) L Vm 


CKV-25] 


for both streamline and turbulent flow. 

According to Eq. XV-19, / is reciprocal to {Re) and by this to Vm for 
streamline flow. Then it follows from Eq. XV-23 that is inde- 
pendent of the velocity. 

For turbulent flow, according to Eqs. XV -20 and 23, h is about pro- 
portional to There is still another influence of velocity, because, 

according to Eq. XV-17, also the factor <j) depends on {Re). But this 
influence is small for great values of (Pe), and for small values of (Pr). 
Water at 100 F, for instance, has (Pr) » 4. With this value, Eq. 
XV-17 yields 4> = 1.83 for {Re) = 40,000, and <^> = 1.77 for {Re) = 
90,000, so that the change in </> is almost negligible. 

Example XV-1. Air at an average absolute pressure of 29.4 Ib/sq in. and 
an average temperature of 70 F is flowing through a pipe line, 1000 ft long, 
3“in. inner diameter, its average velocity being 30 ft/sec. 

Thermocouples attached to the wall show an average temperature of 60 F. 
Assuming the inner surface to be smooth, and neglecting entrance flow dis- 
turbances, find the pressure drop, the heat loss per hour of the air, and the film 
coefficient of heat transfer on the inner pipe w^all. 

Solution: The mean temperature difference is 0m = = 70 — 60 = 

10 F. This is so small that aU properties of the air may be based on the 
average temperature im = 70 F. 



158 HEAT TRANSFER AND PRESSURE DROP [XVnt 

The first thing to be found out is, whether the flow is streamline or turbulent, 
(ije) = 

- 30(3600) = 108,000 ft/hr 
D = = 0.25 ft 


7 = 0.15 Ib/cu ft (at 29.4 Ib/sq in.) 
p = 0.15/32.2 = 0.00466 slug/cu ft 


Because the pressure does not influence much the dynamic viscosity, p can 
be taken from Table VI-2. 

p = 0.105(10-^)36002 = 0.00136 slug ft^i hr“i 


(Re) = 


108,000(0.25) 0.00466 
0.00136 


92,500 


Thus, the flow is turbulent, and in calculating Ap from Eq. XV~24, the fric- 
tion factor / may be taken from Eq. XV-20. 


0.08 



0.0046 


Ap = 4(0.0046) ^ 0-00466(30)^ ^ 1.07 Ib/sq in. 

U,2u 2 

The heat loss may be found either from Eq. Xy-25 or from Eq. XV-21, 
with <j& to be calculated by means of Eq. XV-17. 

iPr) = f 

Cp = 0.24(32.2) = 7.73 B slug-i F-i 

k = 0.0148 B hr-i ft-i P-i (from Fig. II-l) 

00013 ^ 

0.0148 

Cp = (0.15) (0.24) = 0.036 B ff^ F-i (see Sect. IV-1) 

0 -^^^ = 0 82 

1 + 750VQ.71 7.5^0.71 1.029 

92,500 /92,500 


A — LtD 1000 7r0.25 = 785 sq ft 
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Substituting all these terms in Eq. XV-21: 

= 0.82(0.0046) 785(30)10 = 22.5 B/see 

2 ( 0 ./ 1 ) 

or 

- 22.5(3600) = 81,000 B/hr 

The film coefficient of heat transfer is found from Eq. X\^-22. 

hn, °° 10.3 B hr-i F"! 

/So (10) 

This result may be checked without using the pressure drop analog\% by 
means of Eq. VIII-2. Substituting the pertinent values in this equation: 

0 01 4S 

km 0.023(92,500)0-8 (0.71)O-^ = 11.1 B hr“i ft~2 F-^ 

0.25 


The agreement of the results, found in such different ways and using 
partly empirical equations and some simplifications, is satisfactory^ 

It is seen that Eq. Vni~2 leads much easier to the result than the 
analogy method. If, however, Eq. VIII-2 were not known, or in other 
cases where direct heat transfer measurements were missing and \vere 
difficult to perform, the indirect method would be very useful. If, for 
instance, the tubes of a heat exchanger are curved, the influence of 
curvature to heat transfer can be calculated from the results of pres- 
sure drop measurements as \nU be sho^Mi in the next section. 


XV-5 Secondary Influences 

The equations of Sect. XV-4 are based on the assumptions that the 
tubes under consideration are perfectly smooth and straight, that no 
entrance disturbances exist, and that the temperature difference is 
small. In practice, there may be deviations from each of these con- 
ditions. 

Deviations from the first and last of them wiU not be dealt with 
here. However, ordinary^ drawn tubes can be considered as smooth, 
and the equations give good approximations also for moderately rough, 
tubes and for moderate temperature differences. 

The local disturbances of the flow near the entrance of a tube give 
rise to an additional pressure drop which, according to White, can be 
considered by adding 

5(/) = 0.1 7 [XV-BO] 


to the value of f, as found by Eq. XV-20. This is said to hold for 
abrupt entrances. 
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At j? = 10,000, / - 0,008. If LIB = 12.5, Eq. XV-26 yields 
d{f) = 0.008. This means that under the above circumstances the 
entrance loss is just as great as the ordinary loss in an additional length 
8(1) = 12.5 D. 

In calculating the heat transfer the same additions may be provided 
although there is not much experimental evidence about this. 

Concerning the influence of curvature of tubes, Whitens experiments 
on the turbulent flow resistance of closely coiled helixes led to: 


/ = 


0.08 


+ 0.012 



tKV-27] 


where D = the inner pipe diameter, 

C == the mean coil diameter. 

This holds from Re = 15,000 to Re == 100,000. 

The heat transfer in this range may be calculated from Eq. XV--15, 
using Eq. XV~27 instead of Eq. XV-20. 

XV-6 Heat Transfer and Fluid Friction on Plane Surfaces 
In a fluid flowing along a plate, the length L of the plate may be used 
as characteristic length so that 

(Re) = IXV-28] 


Colburn (Ref. XV-fl), in a correlation of experimental data, dis- 
tinguishes between a streamline or viscous region up to (Re) == 20,000 
and a turbulent region above this. 

Using the definition of /, given by Eq. XV-18, he shows that for a 
streamline flow: 


and for turbulent flow: 


/ = 


1.32 

V'(Pe) 

CXV-29] 

0.072 

</{Re) 

[xv-ao] 


These equations correspond to Eqs. XV-19 and 20 for tubes. 

For both, viscous and turbulent range, Colburn further recommends 
the formula 


^ ^ y ^p - C-D 


[XV-31] 


^vhich corresponds to Eq. XY-2Z. 
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In the flow of fluids across tubes, as in tube banks, at the front and 
rear of the obstacle, forces originate which have nothing to do with 
the exchange of momentum. Therefore, the simple relation between 
surface friction and heat transfer fails. Concerning this, the reader is 
referred to the quoted papers of White (Ref. XV-2) and Colburn 
(Ref. XV-5). 


PROBLEMS 

XV-1. Sulfuric acid is pressed through a capillary tube, in. wide and 10 ft long, 
with well-rounded entrance opening. The pressure difference between the inlet and 
outlet is 10 Ib/sq in., the mean temperature of the fluid is 40 F. How much heat 
energy per degree Fahrenheit difference between wall temperature and mean liquid 
temperature will the acid be able to pick up in one hour? 

The density of the acid is 1.84 gram/cu cm, the dynamic \dscosity is 45 centipoise, 
the thermal conductivity is 0.145 B hr”^ ft~^ F~^, and the Prandtl number is 240. 

XV-2. Water at 60 F average temperature flows through a straight piece of drawn 
copper tubing, 1 in. diameter, 20 ft long, and then through a coil of the same pipe, 
having twepty turns of 15-in. mean diameter, which is located in a liquid bath of con- 
stant temperature. The pressure drop in the straight piece of pipe is found to be equal 
to H-in. mercury column, and the wall temperature of the coil is held at the constant 
temperature of 80 F. Calculate the coefficient of heat transfer at the inside surface 
of the coil (a) using the analogy between surface friction and heat transfer and (6) 
without using this analogy. 

X V-3 . Air at standard atmospheric pressure and 50 F average temperature streams 
along a plane plate, 5 ft wide, with an average velocity of 50 ft/sec. A length of 15 
ft of this plate is held at 100 F by a heating device. Neglecting all secondary influ- 
ences, like that of the edges of the plate or of the increase of temperature of the air 
along the plate, calculate how much heat per hour is given up to the streaming air. 
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Abbreviations, British thermal unit, 5 
degree Fahrenheit, 5 
degree centigrade, 6 
Absorption of radiation, 116 
Absorptive power of a surface, 119 
Absorptivity, 120; see also Emmissivity 
definition, 117 
Adams, 49 
Air bubbles, 111 

Air layers, free convection in, 133 
Air spaces, heat transfer through, 133 
in buildings, 134 
Aluminum foil, 134 
Amorphous substance, 7 
Amplitude of oscillation, 51 
Archimedes' law, 69 
Area factor, 126 

for parallel and opposite squares, 12$ 
Asbestos paper insulation, 13, 14 

B 

Baffles, 102 
Biot, 2 

Bismuth versus bismuth-tin thermo- 
couples, 140 

Black body, concept of, 117 
definition, 117 
emitter of energ}^, 120 
measure of absorptivity, 118 
realization of, 117 
Bliss, 113 
Boiler, 97 

Boiling, detour of heat flow in, 112 
of water, formulas under free convec- 
tion conditions, 112 
theory of, 112 
types of, 111 
film. 111 
nuclear, 111 
Boltzmann, 3, 120 
Bonilla, 113 


British thermal unit, 5 
Bubbles, air, 111 
vapor, 112 
stirring effect of, 112 
Building materials, 19 
Buoyancy, 68 


Callend.\e, 155 
Calorie, 6 
Cattle hair, 11 
Celite, 15 
Centipoise, 62 

Coefficient, combined, of convection and 
radiation, 131 

equivalent, of heat radiation, 131 
of free convection, 74 
of friction, 156 

of heat transfer, 3; see also Film co- 
efficient of heat transfer 
overall, 94, 96, 109 
of thermal expansion of ideal gases, 69 
Colburn, 160, 161 
Condensate layer, 109 
Condensation, dropwise, 109 
conditions of existence, 110 
film, 109 

conditions of existence, 110 
mechanism, 110 
theory of, 110 

influence of surface conditions, 110 
mixed, 110 
Condenser, 97 

.Conduction, in steady flow, fundamental 
equation, 2, 45 
definition, 2, 25 

equation of temperature field, 45 
in unsteady state, 36 
equation, 45 
heat stored, 2 

temperature distribution, 38, 44 
through composite cylinder wall, 30 
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Conduction, through composite plane 
wall, 27 

through homogeneous cylinder wall, 29 
through homogeneous plane wall, 25 
transient state, 36 
Conductivity, electrical, 149 
thermal, apparent, 152 
conversion factors, 6 
definition, 2 
fictitious, 152 

influence of temperature, 9 
mean value of, 33 
of alloys, 21 

of amorphous solid substances, 7 

of building materials, 21, 137 

of diatomaceous earth, 17 

of Douglas fir, 19 

of electrical insulators, 7 

of gases, 8, 9, 24, 139 

of insulating materials, 10, 15, 16, 17 

of iron alloys, 22, 23 

of liquids, 7, 24, 139 

of magnesia, 85%, 15 

of metals, 21, 22, 136 

of non-iron alloys, 22 

of refractory materials, 17, 18, 19 

of solids, 24 

of water, 9, 10 

of wood, 19, 20 

related to electric conductivity, 7 
units, 6, 26 

Configuration factor, 126 
Conservation of energy, 37 
Convection, artificial, see Convection, 
forced 

free, 68, 72, 73 
equation of, 67, 69 
in boiling, 112 
on horizontal cylinders, 69 
on horizontal pipes, 74, 75 
on horizontal plates, 73 
on vertical surfaces, 67, 73 
peculiarities, 72 
variables for, 67 
forced, 83 

dimensionless groups, 79 
equation of, 67 
importance, 78 
in boiling, 112 
peculiarities, 78 


Convection, influence of, in temperature 
measurements, 145 
natural, see Convection, free 
Conversion of constant factors in equa- 
tions to different units, 113 
Cooling of viscous oils, 81 
Cork, 11 

Crystalline substance, 7 
Curvature, influence of, to heat transfer, 
159 

Cylinder method, co-axial, of measuring 
thermal conductivity, 138, 139 
of measuring thermal conductivity of 
insulating materials, 138 

D 

Degree centigrade, 6 
Degree Fahrenheit, 5 
Density, apparent, 7; see also Weight, 
apparent specific 
of emission, 119 
of water, 4 

product of specific heat and, 40 
Dewar vessels, 134 
Diffusivity, thermal, 40 
Dimensional analysis, 66 
advantage of, 69, 70 
limitations of, 70, 76 
simplifying assumptions, 70 
Dimensional homogeneity, 66' 
Dimensional soundness of equations, 113 
Dimensionless groups, 61, 70 
denotation, 70 

Dimensionless magnitude, 63 
Dimensions, check of physical, 5 
equations of, 68 
Drops, formation of, 109 

E 

Earth, diatomaceous, 14, 16 
Energy, conservation of, 37 
storage of, 2, 39 
Electric coil, 55 
average temperature of, 57, 58 
durability of, 58 

maximum temperature of, 57, 58 
temperature distribution in, 57 
Electric resistance of metal, temperature 
coefficient, 57 

Electro-magnetic waves, 116, 117 
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Emission, 119 
Emissivity, 120 
definition, 120 
factor, 126 

of aluminum, influence of air, 134 
of metallic surfaces, 122 
of non-metallic surfaces, 122 
Equations, for forced convection, 67 
for free convection, 67 
of gases on a '^’ertical wall, 69 
of dimensions, 68 

Errors, in temperature measurements, 147 
thermometric, 144, 145 
Evaporation on a hot plate, 111 
Experimental determination, see Meas- 
urement 

Exponents, of basic units in dimensional 
analysis, 66 
use of, 73 


Fanning’s factor, 156 
Film boiling, 111 

Film coefficient of heat transfer, average, 
156 

calculated from pressure drop, 157 
definition, 3 

determination without heat flow or 
temperature, 157 

for condensation on horizontal tubes 
and vertical walls, 110 
for dropwise condensation, 109 
for film condensation, 109 
for forced convection, 79 
for free convection, 68, 73 
in boiling at atmospheric pressure, 112, 
113 

in boiling- at different pressures, 113 
quantities influencing, 60 
Film condensation, 109 
conditions of existence, 110 
theory of, 110 
Fins on a plane plate, 93 
Flow, across tubes, 161 
laniinar, 61 
of fluids, 61 
of gases, 145 
streamlined, 61, 156 
transition, -from laminar to ttirbul^t, 
65 ■ • . .r 


Flow, turbulent, 61 

empirical equation, 156 
velocity and temperature ratio for, 
153 

viscous, 61 

Flow resistance, of coiled helixes, 160 
influence of curvature of tubes on, 
160 

Fluid, see Flow 
Fluid friction, 61 

analogy to heat transfer, 149, 150 
on plane surfaces, 160 
Fourier's equation of heat conduction, 2 
Frequency, 51 
Friction, see Fluid friction 
Fritz, 113 

G 

Gas, diatomic, 74, 75 
flowing, temperature measurement of, 
145 
I ideal, 69 

thermal expansion coefficient of, 69 
Gas cells. 111 
Gas film theory, 76, 83 
Gas friction, 82 
Gas radiation, bands of, 119 
Gauss's error integral, 46 
Glass vessel, evacuated, 134 
Glass wool, 11, 12 
Glassy substance, 7 
Grashof number, 70 
Groeber, 105 
Guard coils, 139 
Guard rings, 137 

H 

Heat, specific, at constant pressure re- 
lated to unit volume, 41 
density and, 40 

transmission of, molecular theory, 6 
transportation of momentum and, 151 
units of, 40 

Heat conduction, at variable conductiv- 
ity, 32 • . 

definition, 1 

steady state, definition, 2, 25 
equations for, 45 ; 
through thermocouple wires, 144 
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Heat conduction, unsteady state, 2, 36 
equation for, 45 

temperature distribution, 38, 44 
Heat conductivity, see Conductivity, 
thermal 

Heat convection, definition, 1, 60 
forced, 78 
free, 68, 72 

in vertical air layers, 133 
Heat energy, 2 
of condensation, 109 
Heat exchange, net, 125 
Heat exchangers, classes, 97 
counter-flow, 97 

effect of conduction and convection in 
85, 93 

economic optimum in design of, 149 
parallel flow, 97 
shell and tube, 102 
Heat flow, definition, 2 
per linear foot, 29 
rate of, 2 

through insulation at unknown surface 
temperature, 132 

Heat loss, of cylindric apparatus in the 
axial direction, 139 
Heat radiation, see also Radiation 
at low temperatures, 122 
between an enclosed body and the en- 
closure, 125 

between a small enclosure and the en- 
closure, 126 

between equal parallel and opposite 
squares, 127 

between parallel black planes, 123 
between parallel planes of different 
emissivity, 124 
by aluminum foil, 134 
by carbon dioxide, 119 
by metallic surfaces, 122 
by non-metallic surfaces, 122 
by oxidized surfaces, 123 
by steam, 119 
definition, 2 
equilibrium, 119 

equivalent coefficient of heat transfer, 
131 

from sun to earth, 116 

fundamental equation, 3 

general equation of net interchange, 127 


Heat radiation, heat exchange by, 123 
influence of, in temperature measure 
ments, 145 
net exchange, 125 

superimposed on conduction and cor 
vection, 130 

transfer mechanism, 116 
velocity, 116 
Heat resistance, 29 

Heat sources within a heat conductii 
body, 55 

Heat transfer, convective, analogy ; 
fluid friction, 149 
molecular viewpoint, 150 
and fluid friction on plane surfaces, 
160 

from horizontal tubes in still air, 131 
from protruding rod, 85 
influence of curvature on, 159 
through air spaces, 133 
Heat transmission, Ixjtwer^n fluids 
through a plane wall, 93 
between fluids through a cylindrirt 
wall, 95 

laws and relation to other physica' 
laws, 149 
molcc\ilar, 6 

momentum exchange and, 151 
Heat wave, 50 
Heated enclosure, 119 
Heating, of viscous oils, 81 
Heating and cooling cycle, 36 
Heilman, 132 
Hottel, 127 


Insinger, 113 

Insulation, for building purposes, 1 1 
for heating and process work, 13 
high-temperature, 14 
in high-ternperature pipe service, 16 
in the power generation field, 14 
low-temperature, 11 
materials of, 10 

apparent specific weight of, 16, 17 
brick of natural Celite, 15 
metallic, 13, 134 
Insulation efficiency, 133 
Intensity of radiation, 118 
Inverse square law, 127 
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Jakob, 111, 113 
Joule, 58 

K 

Kilo-calorie, definition, 6 
Kirchhoff, 117 
Kirchhoff's law, 118, 120 


Lanainar flow, 61 
Langmuir, 76, 83 
Layer of condensate, 109 
Length, characteristic, 70 
Links, 113 

Log mean temperature difference, 98, 
100 

correction factors for, 101 
Lorenz, 69 

M 

Magnesia, 85%, 14, 16 
Mass velocity, 64 
McAdams, 132 
McMillan, 132 
Mean temperature difference, 98 
Mean value of thermal conductivity, 33 
Measurement, of emissivities, 139, 140 
of temperatures, of flowing gases, 143, 
145 

of furnaces, 147 
of surfaces, 142 

of thermal conductivity, of building 
and insulating materials, 137 
of gases, 139 
of liquids, 139 
of metals, 136 
Metal foils, 13 
Mixing movements, 82 
Mixing-cup temperature, 78 
Momentum, transportation of heat and, 
151 

N 

Newton’s equation of convective heat 
transfer, 3 

Newton’s equation of fluid friction, 150 
Nuclei, 111 
Nusselt, 67, 110 
Nusselt number, 70 


0 

Ohm’s law, 58 

Oils, viscous, cooling and heating of, 81 
Opaque body, 117 
emissivity of, 120 
Oscillation, amplitude of, 51 
Overall coefficient of heat transfer, 94, 
96, 109 


Pannell, 152, 153 
Perry, 113 
Phase, 109 
change of, 109 
Pipes, see Tubes 
Planck’s law, 120 
Poise, 62 
Poiseuillb, 62 
Pores, filled with water, 7 
insulating effect of, 7 
Pound force, 4 
Pound mass, 4, 64 
Poundal, 64 
Prandtl, 82, 83 
Prandtl number, for gases, 80 
physical property, 70, 151 
Pressure drop, 157 
in a tube, 65 
in streamline flow, 62 
near entrance of tube, 159 
Pyrometers, optical, 147 
selective, radiation, 147 
suction, 145 
total radiation, 147 

R 

Radiant heat, see Heat radiation 
Radiant interchange, net, 124 
Radiation, see also Heat radiation 
exchange of, 123 
general laws of, 116 
intensity, 118, 122 ■ 
monochromatic, 118, 120 
reflection of long-wave and short-wave, 
122 

total, 120 
wavelength of, 119 
Radiation pyrometers, 147 
Radiation receiver, 140 
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Radiator, gray, 121 
non-black body, 147 
Rate of heat flow, 2 
Reflectivity, 117 

Refractory materials, 17; see also Insula- 
tion 

Rbiher, 143 

Resistance thermometer, 139 
Reynolds, 61 

analogy between heat transfer and 
pressure drop, 149 
basic idea, 150 

modification for streamline flow, 155 
modification for turbulent flow, 156 
restricting condition for, 151 
restrictions of, 152, 153 
Reynolds number, 61, 65 
critical, 65 
definition, 70 
Rice, 77, 83 
Rock wool, 11, 12 
Rod, end uninsulated, 91 
protruding, 86 

S 

Shields, protecting, in thermometry of 
flowing gases, 145 
Sieder, 80 

Similarity, principle of, 61, 70 
Slag wool, 11 
Slug, 4 

Specific heat; see Heat 
Specific weight; see Weight 
Steady state; see Heat conduction 
Steam bubbles, see Vapor bubbles 
Stefan, 120 

Stefan-Boltzmann^s, constant, 3 
law, 3, 120 
generalization, 121 
Stirring effect of vapor bubbles, 112 
Storage of heat energy, 2, 39 
Streamline flow; see Flow 
Substances, amorphous solid, 7 
crystalline, 7 
glassy solid, 7 
Sunlight, 122 

Surface, absorptive power of, 119 
black, 3; see also Black body 
emissive power of, 110 
projections of, 85 


Surface, reflecting, 134 
Surface film theory, 82 
of forced heat convection, 83 
of free heat convection, 76 
Surface temperature, 72 
amplitude, 51 
periodic change, 50 
sudden change, of cylinder, 49 
of infinitely thick plane wall, 46 
of plane wall of finite thickness, 4 
of sphere, 49 


Tate, 80 

Temperature, above the liquid level o 
boiling liquid, 112 

and velocity, ratio for turbulent flow 
153 

comparison of distribution of, 152 
distribution of, across an air stream 
153 

coefficient of electric resistance, 57 
deviation from average, 51 
difference, definition, 73, 78 
in heat exchangers, 101 
distribution, 91 
along a rod, 86 
along a thermometer well, 143 
in unsteady state, 38, 44 
inside vapor bubbles, 112 
mixing-cup, 78 
of a boiling liquid, 112 
of a flowing liquid, 72 
of surface, see Surface temperature 
sudden exposure of wall, 105, 106 
Temperature field, 45 
Temperature gradient, 37 
Temperature measurements, 142, 147 
errors in, 144, 145, 147: , 

Temperature radiation, 118 
Thermal conductivity, see Conductivity 
Thermal resistance, 29 
Thermocouples, arrangement of, 143 
bismuth versus bismuth-tin, 140 
heat losses by, 143 
Thermos bottles, 134 
Time lag, 50 ■ , ■ ■ . ' 

Transient state, see Heat conduction, 
unsteady state 
Transmissivity, 117 
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Transportation of momentum and heat, 
151 

Cube bank, 161 

heating of fluids flowing normal to, 
82 

Cubes, bare horizontal, 132 
co-axial, 95 
concentric, 95 

flow disturbance near entrance of, 159 
drawn, 159 

heat transfer for streamline flow inside, 
80 

heat transfer for turbulent flow inside, 
79 

heating of fluids flowing normal to, 81 
insulated horizontal, 132 
Turbulent flow, see Flow 
Twin-plate arrangement, 137, 138 

U 

Units, 3, 4 

basic, exponents of, in dimensional 
analysis, 66 
consistency of, 3 

Unsteady state,* see Heat conduction 


Vapor bubbles, 112 
Vapor layer, 111 

Vaporization, cells of, 111; see also Boil- 
ing 

Velocity and temperature, comparison 
of distributions of, 152 


Velocity and temperature, distribution 
of, across an air stream, 153 
ratio of, for turbulent flow, 153 
ViDMAR, 57 
Viscosity, absolute, 63 
apparent, 152 
dynamic, 63 
definition, 61 

influence of temperature, 81 
molecular interpretation, 150 
of air, 63 
of water, 62 
physical dimensions, 62 
units, 62, 64 
fictitious, 152 
kinematic, 63 
Visibility, relative, 63 

W 

Wall, conduction, see Conduction 
Warming, of a plane plate, 36 
Wavelength, of radiation, 119 
White, 149, 153, 156, 160, 161 
Williamson, 49 

Wires, heat convection for normal flow 
to, 81 

Weight, apparent specific, of building 
materials, 21 

of insulating materials, 16, 17 
of refractory materials, 21 
of wood, 20 
specific, 4 
of water, 4 



